Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,227,173 A→G 67.7% V389A (GTC→GCC)  glpK ← Glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,227,1730AG67.7% 89.3 / 56.7 65V389A (GTC→GCC) glpKGlycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (12/9);  new base G (19/25);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

CGCATCTGCGGCCGCTCCACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGG  >  ACB143/1227038‑1227305
                                                                                                                                       |                                                                                                                                    
cgcATCTGCGGCCGCTCCACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGa                                                                                                                                      >  4:350121/1‑136 (MQ=255)
cgcATCTGCGGCCGCTCCACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGa                                                                                                                                      >  6:198332/1‑136 (MQ=255)
                ccACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTgg                                                                                                                        >  2:339575/1‑134 (MQ=255)
                 cACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAg                                                                                                                     <  1:147086/136‑1 (MQ=255)
                 cACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAg                                                                                                                     <  5:99052/136‑1 (MQ=255)
                 cACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAg                                                                                                                     <  1:93860/136‑1 (MQ=255)
                    gCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCg                                                                                                                  <  6:334141/136‑1 (MQ=255)
                      aGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGAt                                                                                                                <  5:380955/136‑1 (MQ=255)
                      aGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGAt                                                                                                                >  6:22564/1‑136 (MQ=255)
                         tGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGa                                                                                                             <  1:9585/136‑1 (MQ=255)
                           ccGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACt                                                                                                           >  5:386120/1‑136 (MQ=255)
                            cGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCTTCGCGGGTCTGGTAGGCGATGGACTc                                                                                                          >  2:131325/1‑136 (MQ=255)
                                 aTGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTg                                                                                                     <  1:236972/136‑1 (MQ=255)
                                 aTGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCGGGTAGGCGATGGACTCCAGTg                                                                                                     >  7:185086/1‑136 (MQ=255)
                                    tCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTgcgg                                                                                                  <  1:314759/136‑1 (MQ=255)
                                    tCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGCCGATGGACTCCAGTgcg                                                                                                   <  5:312190/135‑1 (MQ=255)
                                     cGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTgcggc                                                                                                 <  1:94085/136‑1 (MQ=255)
                                     cGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTgcggc                                                                                                 >  4:156508/1‑136 (MQ=255)
                                       gCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGgcgc                                                                                               >  6:266118/1‑136 (MQ=255)
                                       gCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGgcgc                                                                                               <  5:171843/136‑1 (MQ=255)
                                       gCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGgcgc                                                                                               >  8:2633/1‑136 (MQ=255)
                                         cTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCgg                                                                                             <  8:18342/136‑1 (MQ=255)
                                            gAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGga                                                                                            <  1:330254/134‑1 (MQ=255)
                                            gAACTGCATGAGGAAGATGTTGGCCACGGCACCTCCGTACACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGGCTGGTAGGCGATGGACGCCAGTGCGGCGCGgatg                                                                                          <  8:363751/136‑1 (MQ=255)
                                                     gAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGTATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTc                                                                                  >  7:226460/1‑135 (MQ=255)
                                                     gAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCAGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTc                                                                                  >  7:330351/1‑135 (MQ=255)
                                                     gAGGAAGTTGTTGGCCACGGCACAACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACACTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGTCGCGGATGATGTGGTCt                                                                                 >  2:281535/1‑136 (MQ=255)
                                                         aaGTTGTTGGCCACGGCACCACCGTACACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGAGATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCTGATGATGTGGTCTACCt                                                                             <  1:152383/136‑1 (MQ=255)
                                                             tgttgGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCAc                                                                         >  6:325019/1‑136 (MQ=255)
                                                              gttgGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCAc                                                                         <  7:44255/135‑1 (MQ=255)
                                                              gttgGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCAc                                                                         <  3:340091/135‑1 (MQ=255)
                                                              gttgGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCAc                                                                         <  1:82987/135‑1 (MQ=255)
                                                               ttgGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGc                                                                       >  1:294335/1‑136 (MQ=255)
                                                                  gCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCAc                                                                    <  2:179066/136‑1 (MQ=255)
                                                                         caccacCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTGCCTTCACGCCACGTGTCAg                                                             >  5:330014/1‑136 (MQ=255)
                                                                               cGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGaa                                                       <  5:263863/136‑1 (MQ=255)
                                                                                   cACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCCGGCCGaacaac                                                   >  3:232067/1‑136 (MQ=255)
                                                                                    aCGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACg                                                  >  4:169517/1‑136 (MQ=255)
                                                                                    aCGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACg                                                  >  5:104851/1‑136 (MQ=255)
                                                                                    aCGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCAGGTGTCAGGCCGAACAACg                                                  >  6:322119/1‑136 (MQ=255)
                                                                                    aCGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGGCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACg                                                  <  5:187051/136‑1 (MQ=255)
                                                                                      gcgcAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGca                                                <  5:266118/136‑1 (MQ=255)
                                                                                      gcgcAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGca                                                <  6:290638/136‑1 (MQ=255)
                                                                                      gcgcAGCTCGGACAGGCGCAGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGca                                                <  5:309650/136‑1 (MQ=255)
                                                                                              cGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACg                                                                        <  6:26143/104‑1 (MQ=255)
                                                                                              cGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACg                                                                        >  5:26143/1‑104 (MQ=255)
                                                                                               ggACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCCCGTGTCAGGCCGAACAACGCCCCACGGGCa                                       >  3:52973/1‑136 (MQ=255)
                                                                                               ggACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGc                                        <  7:240303/135‑1 (MQ=255)
                                                                                                gACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCAt                                      <  6:335261/136‑1 (MQ=255)
                                                                                                    ggCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAg                                                                                                                     >  5:100562/1‑53 (MQ=255)
                                                                                                    ggCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAg                                                                                                                     <  6:100562/53‑1 (MQ=255)
                                                                                                     gcgcTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACggg                                 <  8:330351/136‑1 (MQ=255)
                                                                                                     gcgcTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACggg                                 <  3:19407/136‑1 (MQ=255)
                                                                                                     gcgcTGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACggg                                 <  1:99624/136‑1 (MQ=255)
                                                                                                         tGGCCACAGTCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCCCGGGCATACGGGTcc                              >  8:359583/1‑135 (MQ=255)
                                                                                                         tGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCAGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTcc                              >  2:276734/1‑135 (MQ=255)
                                                                                                                  tCCTGCTGCATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACgggg                    >  6:139841/1‑136 (MQ=255)
                                                                                                                    ctgctgCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGc                   <  5:198332/135‑1 (MQ=255)
                                                                                                                    ctgctgCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGc                   <  3:350121/135‑1 (MQ=255)
                                                                                                                       ctgcATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCAccc               <  5:371599/136‑1 (MQ=255)
                                                                                                                        tgcATGGCGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCAcc                >  1:7975/1‑134 (MQ=255)
                                                                                                                             ggCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCg         >  3:15350/1‑136 (MQ=255)
                                                                                                                               cGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCgg        >  4:379463/1‑135 (MQ=255)
                                                                                                                               cGTCGAGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCgg        >  8:298283/1‑135 (MQ=255)
                                                                                                                                 tCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTa      >  5:289475/1‑135 (MQ=255)
                                                                                                                                  cGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTa      >  5:366074/1‑134 (MQ=255)
                                                                                                                                    aGGGCATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAgg  >  3:55018/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                    
CGCATCTGCGGCCGCTCCACGCAGGTGCCGAGGATGTCGGCCTGGAACTGCATGAGGAAGTTGTTGGCCACGGCACCACCGTCCACGCGCAGCTCGGACAGGCGCTGGCCACAGTCCTGCTGCATGGCGTCGAGGACATCGCGGGTCTGGTAGGCGATGGACTCCAGTGCGGCGCGGATGATGTGGTCTACCTTCACGCCACGTGTCAGGCCGAACAACGCACCACGGGCATACGGGTCCCAGTACGGGGCACCCAGGCCGGTAAAGG  >  ACB143/1227038‑1227305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: