Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 373,333 C→G 5.5% P138R (CCA→CGA) ‡ mhpF → acetaldehyde‑CoA dehydrogenase II, NAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913373,3330CG 5.5% 468.2 / 16.1 127P138R (CCA→CGA) ‡mhpFacetaldehyde‑CoA dehydrogenase II, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (71/49);  new base G (3/4);  total (74/53)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 128 total aligned reads displayed.

CGCTTTACGCGAAGCGAAACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAGTGGTGGGCG  >  NC_000913/373184‑373479
                                                                                                                                                     |                                                                                                                                                  
cGCTTTACGCGAAGCGAAACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTcc                                                                                                                                                    >  2:1415281/1‑150 (MQ=255)
        gcgaagcgaaACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAAtggtgg                                                                                                                                            >  2:2291901/1‑150 (MQ=255)
          gaagcgaaACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTggcg                                                                                                                                          >  1:1995621/1‑150 (MQ=255)
           aagcgaaACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTGGAATGGTggcct                                                                                                                                         <  1:2302361/150‑3 (MQ=255)
           aagcgaaACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTGGAATGGTggcct                                                                                                                                         <  1:1630883/150‑3 (MQ=255)
           aagcgaaACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTggcgg                                                                                                                                         <  1:451987/150‑1 (MQ=255)
                   ccGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCa                                                                                                                                 >  1:1817349/1‑150 (MQ=255)
                    cGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCAc                                                                                                                                <  2:556102/150‑1 (MQ=255)
                      gATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCAcgc                                                                                                                              >  2:1434126/1‑150 (MQ=255)
                           tCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGgc                                                                                                                         >  2:368835/1‑150 (MQ=255)
                           tCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGgc                                                                                                                         >  1:1381224/1‑150 (MQ=255)
                           tCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGgc                                                                                                                         >  1:1381450/1‑150 (MQ=255)
                             gCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGgcgc                                                                                                                       <  1:667564/150‑1 (MQ=255)
                             gCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGgcgc                                                                                                                       <  1:2477165/150‑1 (MQ=255)
                               ttAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCgt                                                                                                                     <  2:909316/150‑1 (MQ=255)
                                   ttGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTc                                                                                                                 >  2:1908593/1‑150 (MQ=255)
                                    tGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCa                                                                                                                >  1:1204726/1‑150 (MQ=255)
                                      aCCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCAtt                                                                                                              >  1:1023486/1‑150 (MQ=255)
                                       ccTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTa                                                                                                             >  1:2107973/1‑150 (MQ=255)
                                       ccTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGGTCAGTa                                                                                                             >  1:2119787/1‑150 (MQ=255)
                                          gACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGc                                                                                                          >  1:2492961/1‑150 (MQ=255)
                                          gACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGc                                                                                                          <  1:562384/150‑1 (MQ=255)
                                          gACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGc                                                                                                          <  1:2429931/150‑1 (MQ=255)
                                           aCGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGcc                                                                                                         >  1:274580/1‑150 (MQ=255)
                                              ccTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGaa                                                                                                      >  1:1464553/1‑150 (MQ=255)
                                              ccTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGaa                                                                                                      >  1:994066/1‑150 (MQ=255)
                                                tgctgcCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAAt                                                                                                    >  2:1616435/1‑150 (MQ=255)
                                                 gctgcCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTTATTACGCCGAAAtt                                                                                                   <  2:1309217/150‑1 (MQ=255)
                                                    gcCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCCGTTTCACGCGTGGCGCGTGGTCATTACGCCGACATTATc                                                                                                >  2:2231207/1‑150 (MQ=255)
                                                    gcCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATc                                                                                                >  1:1940008/1‑150 (MQ=255)
                                                    gcCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATc                                                                                                >  2:208384/1‑150 (MQ=255)
                                                     ccATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCg                                                                                               >  2:180395/1‑150 (MQ=255)
                                                     ccATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCg                                                                                               >  2:1458218/1‑150 (MQ=255)
                                                     ccATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCg                                                                                               >  2:2397070/1‑150 (MQ=255)
                                                        tCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCtt                                                                                            >  2:2044377/1‑150 (MQ=255)
                                                        tCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCtt                                                                                            >  2:205362/1‑150 (MQ=255)
                                                            cccTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTAt                                                                                        <  2:1817349/150‑1 (MQ=255)
                                                            cccTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTAt                                                                                        <  1:1388503/150‑1 (MQ=255)
                                                             ccTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATc                                                                                       >  2:361228/1‑150 (MQ=255)
                                                                 aCTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCa                                                                                   <  2:213980/150‑1 (MQ=255)
                                                                   tGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGt                                                                                 >  2:1886326/1‑150 (MQ=255)
                                                                   tGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGt                                                                                 >  2:1886323/1‑150 (MQ=255)
                                                                   tGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGt                                                                                 >  1:725963/1‑150 (MQ=255)
                                                                      gTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTaaa                                                                              >  1:2435159/1‑150 (MQ=255)
                                                                      gTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTaaa                                                                              >  1:2139388/1‑150 (MQ=255)
                                                                      gTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTaaa                                                                              >  2:1811622/1‑150 (MQ=255)
                                                                      gTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTaaa                                                                              >  1:2464772/1‑150 (MQ=255)
                                                                        gCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATc                                                                            <  2:60592/150‑1 (MQ=255)
                                                                         ccGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCt                                                                           >  1:2437215/1‑150 (MQ=255)
                                                                          cGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTg                                                                          >  2:350740/1‑150 (MQ=255)
                                                                           ggtggtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCtac                                                                         >  1:1174378/1‑148 (MQ=255)
                                                                           ggtggtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGc                                                                         <  1:521813/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCCGTAAATCTGCCGGa                                                                     >  1:1875236/1‑150 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  2:1995621/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  1:144557/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  2:1208239/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  2:507898/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  1:1171222/150‑1 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     >  1:683363/1‑150 (MQ=255)
                                                                               gtTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGa                                                                     <  1:498681/150‑1 (MQ=255)
                                                                                    ccTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTgg                                                                <  2:568552/150‑1 (MQ=255)
                                                                                    ccTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTgg                                                                <  2:568436/150‑1 (MQ=255)
                                                                                      tCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCa                                                              >  1:1450304/1‑150 (MQ=255)
                                                                                      tCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGGGTTCAGTACGCCGAAATTCTCGCTTCTCTCGCCAGTAAATCTGCCGGACCTGGcc                                                              >  2:220539/1‑149 (MQ=255)
                                                                                           gCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGt                                                         >  1:2038207/1‑150 (MQ=255)
                                                                                             gAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGc                                                       <  2:1651585/150‑1 (MQ=255)
                                                                                                    gATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATc                                                <  1:102028/150‑1 (MQ=255)
                                                                                                     aTCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCg                                               >  2:2543741/1‑150 (MQ=255)
                                                                                                       cAACTGAACGTCTACAGGGTCACCTGCGGCGGCCAGGCCACCCTTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCTGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGAt                                             <  2:267232/150‑1 (MQ=255)
                                                                                                       cAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGAt                                             <  2:2423647/150‑1 (MQ=255)
                                                                                                         aCTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGa                                           >  2:1621006/1‑150 (MQ=255)
                                                                                                         aCTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGa                                           >  2:565960/1‑150 (MQ=255)
                                                                                                            gAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCACTATCGATGAAtt                                        >  1:1572340/1‑150 (MQ=255)
                                                                                                                gTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCGTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAAGAGCGATGAATTTACg                                    >  1:2071117/1‑150 (MQ=255)
                                                                                                                 tCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACgg                                   >  2:364422/1‑150 (MQ=255)
                                                                                                                 tCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACgg                                   >  1:693551/1‑150 (MQ=255)
                                                                                                                 tCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACgg                                   >  2:1025741/1‑150 (MQ=255)
                                                                                                                  cAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGa                                  <  2:838178/150‑1 (MQ=255)
                                                                                                                        ggTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCACTATCGATGAATTTACGGAAACCAc                            >  2:934708/1‑150 (MQ=255)
                                                                                                                        ggTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATGACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCAc                            >  2:2228853/1‑150 (MQ=255)
                                                                                                                          tCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGACACCACta                          >  2:2416724/1‑149 (MQ=255)
                                                                                                                          tCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACtt                          <  1:1818262/150‑1 (MQ=255)
                                                                                                                          tCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTAAGGAAACCAcat                          >  1:657289/1‑148 (MQ=255)
                                                                                                                            aCCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGGTCGCGTGCCAATATCGATGAATTTACGGAAACCACTTcc                        <  2:21811/150‑1 (MQ=255)
                                                                                                                            aCCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTcc                        <  2:21737/150‑1 (MQ=255)
                                                                                                                              cTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCg                      <  1:1228787/150‑1 (MQ=255)
                                                                                                                                      gccaggccaCCATTCCAATGGTGGCGTAAGTTTCACGAGTGGCGCGTTTTCATTACGCCGAAAATATCGCTTCTATCGCCAGTAAATCTGCAGGACCTGGCACGCGTGCCAAAATAGATGAATTTACGGAAACCACATACCGatcacttt              >  1:1016575/1‑143 (MQ=255)
                                                                                                                                      gccaggccaCCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTg              >  1:1016578/1‑150 (MQ=255)
                                                                                                                                       ccaggccaCCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGa             <  2:479767/150‑1 (MQ=255)
                                                                                                                                        caggccaCCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTACCCGAGCCATTGaa            >  2:1882802/1‑150 (MQ=255)
                                                                                                                                         aggccaCCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAg           >  1:300857/1‑150 (MQ=255)
                                                                                                                                         aggccaCCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAg           <  1:1495685/150‑1 (MQ=255)
                                                                                                                                           gccaCCAGTCCAATGGTGGCGGCAGTGTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAgtg         <  2:2520441/150‑1 (MQ=255)
                                                                                                                                            ccaccaTTGGAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAgtgg        >  1:188837/1‑150 (MQ=255)
                                                                                                                                            ccaccaTTGGAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCAGTGAAgtgg        >  1:1496888/1‑150 (MQ=255)
                                                                                                                                            ccaccaTTGGAATGGTGGCGGCAGTTTCACGCGTGGCGCGGGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAgtgg        >  1:623811/1‑150 (MQ=255)
                                                                                                                                            ccaccaTTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAgtgg        <  2:2232070/150‑1 (MQ=255)
                                                                                                                                              accaTTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAgtggtg      <  1:364422/150‑1 (MQ=255)
                                                                                                                                                 aTTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAGTGGTGggc   >  1:410980/1‑150 (MQ=255)
                                                                                                                                                  ttCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCGATTGAAGTGGTggggg  >  2:1012153/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCTTTACGCGAAGCGAAACCGGATATTCGCTTAATTGACCTGACGCCTGCTGCCATCGGCCCTTACTGCGTGCCGGTGGTTAACCTCGAGGCGAACGTCGATCAACTGAACGTCAACATGGTCACCTGCGGCGGCCAGGCCACCATTCCAATGGTGGCGGCAGTTTCACGCGTGGCGCGTGTTCATTACGCCGAAATTATCGCTTCTATCGCCAGTAAATCTGCCGGACCTGGCACGCGTGCCAATATCGATGAATTTACGGAAACCACTTCCCGAGCCATTGAAGTGGTGGGCG  >  NC_000913/373184‑373479

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