Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,114,537 G→T 5.7% intergenic (‑65/+32) fpr ← / ← glpX ferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,114,5370GT 5.7% 708.5 / 19.4 190intergenic (‑65/+32)fpr/glpXferredoxin‑NADP reductase; flavodoxin reductase/fructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base G (92/87);  new base T (4/7);  total (96/94)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 194 total aligned reads displayed.

GCACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTG  >  NC_000913/4114388‑4114686
                                                                                                                                                     |                                                                                                                                                     
gCACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCg                                                                                                                                                       >  2:874835/1‑150 (MQ=255)
   cGGGGGCGTGATCGGTGAGTCTAAACAGGGCGTCGGTCCAGTTCTGCACTTGAGTGACTGTGGCTGTTACCCTATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAg                                                                                                                                                    <  2:1658956/150‑1 (MQ=255)
          gttaaCGGTGAGACTAAACCGGGCGTCGGTCCAGTTCTGCAATTTTGTGACTTTGAATGTTACCCAATCATAAATTTTTTTCTCCTTTTTTTATTGACTTGCCTTTTCTTCGTTCTCCTAACAAAAAAATTCCAGTCCCTATGGACTGGa                                                                                                                                             <  1:2190413/147‑1 (MQ=255)
          gTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTTGTGACTTTGCCTGTTACCCAATCAGCCTTGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGa                                                                                                                                             <  1:2076213/150‑1 (MQ=255)
           tGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTGAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaa                                                                                                                                            <  2:1591732/150‑1 (MQ=255)
               cGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCTTCGGGACTGGaatttt                                                                                                                                        <  2:518420/150‑5 (MQ=255)
                  tGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCaa                                                                                                                                     >  1:1317496/1‑150 (MQ=255)
                  tGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCaa                                                                                                                                     >  1:183132/1‑150 (MQ=255)
                      aatAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATTTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCTAAAAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGa                                                                                                                                 <  1:967197/148‑1 (MQ=255)
                      aCTAAACCGGGCGTCGGTCCAGTTCTGCACTTTCGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGa                                                                                                                                 <  1:1334033/150‑1 (MQ=255)
                      aCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGa                                                                                                                                 >  2:2572963/1‑150 (MQ=255)
                      aCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGa                                                                                                                                 >  2:2291528/1‑150 (MQ=255)
                       cTAAACCGGGCGTCGGTCCAGTTCTGCACTTTTGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                                <  1:539496/150‑1 (MQ=255)
                       cTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGat                                                                                                                                <  1:1080229/150‑1 (MQ=255)
                         aaaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatca                                                                                                                              <  1:1239182/150‑1 (MQ=255)
                          aaCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGatcaa                                                                                                                             >  1:2429305/1‑150 (MQ=255)
                            cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTAGTTCTCCGAACAAAAAAATTCCAGTACCGAAGGACTGGAAGGCTCAATCGatcaaat                                                                                                                           >  1:560543/1‑150 (MQ=255)
                                   tCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCag                                                                                                                    >  1:1437942/1‑150 (MQ=255)
                                   tCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGAACAAATCAATCag                                                                                                                    >  1:1974143/1‑150 (MQ=255)
                                    cGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCaga                                                                                                                   >  1:742582/1‑150 (MQ=255)
                                     ggTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCagag                                                                                                                  <  2:1058747/150‑1 (MQ=255)
                                       tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                >  2:1625303/1‑150 (MQ=255)
                                       tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGa                                                                                                                >  2:1671791/1‑150 (MQ=255)
                                        ccAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAt                                                                                                               >  1:837624/1‑150 (MQ=255)
                                        ccAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGAt                                                                                                               <  1:642616/150‑1 (MQ=255)
                                                gCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCt                                                                                                       <  2:338605/150‑1 (MQ=255)
                                                 cACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                      >  2:1224645/1‑150 (MQ=255)
                                                 cACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                      >  2:2106206/1‑150 (MQ=255)
                                                 cACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                      <  2:2353528/150‑1 (MQ=255)
                                                 cACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                      <  2:925423/150‑1 (MQ=255)
                                                 cACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTg                                                                                                      >  1:2301395/1‑150 (MQ=255)
                                                  aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGc                                                                                                     >  2:574408/1‑150 (MQ=255)
                                                  aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTCCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGc                                                                                                     >  1:1110531/1‑150 (MQ=255)
                                                   cTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCa                                                                                                    <  1:1857524/150‑1 (MQ=255)
                                                    tttAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCCATCGATCAAAACAATCAGAGGATGTGCACCTGCAt                                                                                                   >  1:1671944/1‑150 (MQ=255)
                                                    tttAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTAGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAt                                                                                                   >  2:1681775/1‑150 (MQ=255)
                                                     ttAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAtt                                                                                                  <  1:1599523/150‑1 (MQ=255)
                                                     ttAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCAtt                                                                                                  <  1:1634758/150‑1 (MQ=255)
                                                        gTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTATTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTcc                                                                                               <  2:523816/150‑1 (MQ=255)
                                                         tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCg                                                                                              <  1:1182851/150‑1 (MQ=255)
                                                         tGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCg                                                                                              <  1:1183028/150‑1 (MQ=255)
                                                          gACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCgg                                                                                             <  2:1251905/150‑1 (MQ=255)
                                                          gACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCgg                                                                                             <  2:2230747/150‑1 (MQ=255)
                                                             tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTc                                                                                          >  2:2439633/1‑150 (MQ=255)
                                                             tttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTc                                                                                          >  2:1258308/1‑150 (MQ=255)
                                                              ttGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCt                                                                                         >  2:171812/1‑150 (MQ=255)
                                                                gCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCttt                                                                                       <  2:2515290/150‑1 (MQ=255)
                                                                  cTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGc                                                                                     >  1:839285/1‑150 (MQ=255)
                                                                        cccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCa                                                                               >  1:874017/1‑150 (MQ=255)
                                                                         ccAATCTGCCTTGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCTTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTTGAAGGCTCAATCGATCAAATCAATCAGATGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                              <  1:177068/150‑1 (MQ=255)
                                                                         ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                              <  1:1066204/150‑1 (MQ=255)
                                                                         ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                              <  2:2406904/150‑1 (MQ=255)
                                                                         ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                              <  1:1672244/150‑1 (MQ=255)
                                                                         ccAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTAACTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAg                                                                              >  1:1991464/1‑150 (MQ=255)
                                                                         ccAAACAGCCATGTTATTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATAGATCAAATCAATCAGAGGATGTGCACCTGAATTTCCGGGTCTATGCGATCCAg                                                                              >  2:320282/1‑150 (MQ=255)
                                                                            aTCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATa                                                                           >  1:472363/1‑150 (MQ=255)
                                                                                  cATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggat                                                                     <  1:1732259/150‑1 (MQ=255)
                                                                                     gTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgga                                                                  <  1:244345/150‑1 (MQ=255)
                                                                                     gTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgga                                                                  <  2:1614725/150‑1 (MQ=255)
                                                                                     gTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGtggatgga                                                                  <  1:421449/150‑1 (MQ=255)
                                                                                         ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                              <  1:2572963/150‑1 (MQ=255)
                                                                                         ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                              <  1:1210116/150‑1 (MQ=255)
                                                                                         ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                              <  1:2573003/150‑1 (MQ=255)
                                                                                         ttCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGa                                                              <  1:2291567/150‑1 (MQ=255)
                                                                                                 ttttGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                      >  2:441854/1‑150 (MQ=255)
                                                                                                 ttttGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                      >  2:2432403/1‑150 (MQ=255)
                                                                                                 ttttGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATgcggcg                                                      >  2:1790911/1‑150 (MQ=255)
                                                                                                   ttgattgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGaa                                                    <  1:376720/150‑1 (MQ=255)
                                                                                                       ttgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGt                                                >  1:2483295/1‑150 (MQ=255)
                                                                                                        tgaCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTg                                               <  1:2106206/150‑1 (MQ=255)
                                                                                                          aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCg                                             <  1:54518/150‑1 (MQ=255)
                                                                                                          aCTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCg                                             >  2:723374/1‑150 (MQ=255)
                                                                                                             tGCCGTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGACGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATCGTGGATGGACTGAATGCGGCGAATGGTGCGTGa                                          <  1:872503/150‑1 (MQ=255)
                                                                                                               ccTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACt                                        >  2:2518991/1‑150 (MQ=255)
                                                                                                                   aTCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                    >  2:421782/1‑150 (MQ=255)
                                                                                                                   aTCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGcc                                    >  2:2166100/1‑150 (MQ=255)
                                                                                                                        cGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGa                               <  1:687700/150‑1 (MQ=255)
                                                                                                                          ttCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATc                             >  2:1364795/1‑150 (MQ=255)
                                                                                                                                gAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCagcagc                       >  1:1416258/1‑150 (MQ=255)
                                                                                                                                     aaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTc                  >  1:2506908/1‑150 (MQ=255)
                                                                                                                                     aaaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTc                  >  2:1339276/1‑150 (MQ=255)
                                                                                                                                      aaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCg                 >  2:751273/1‑150 (MQ=255)
                                                                                                                                      aaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCg                 <  2:1856916/150‑1 (MQ=255)
                                                                                                                                      aaaaaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGGTTCg                 >  2:1689318/1‑150 (MQ=255)
                                                                                                                                         ccttCCAGTCCTTCGGGACTGGAATGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTa              >  1:294135/3‑150 (MQ=255)
                                                                                                                                         aaTTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTa              >  2:374392/1‑150 (MQ=255)
                                                                                                                                           ttCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGt            <  1:351165/150‑1 (MQ=255)
                                                                                                                                               aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg        >  1:1644267/1‑150 (MQ=255)
                                                                                                                                               aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg        >  1:701491/1‑150 (MQ=255)
                                                                                                                                               aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg        >  1:1502985/1‑150 (MQ=255)
                                                                                                                                               aGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTCCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTcgcg        >  2:142102/1‑150 (MQ=255)
                                                                                                                                                 tCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGat      >  2:1136546/1‑150 (MQ=255)
                                                                                                                                                 tCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGat      >  2:1292330/1‑150 (MQ=255)
                                                                                                                                                 tCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGat      >  2:1862377/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGatat    >  2:531171/1‑150 (MQ=255)
                                                                                                                                                    cGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATAtt   >  1:496162/1‑150 (MQ=255)
                                                                                                                                                    cGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATAtt   >  1:109380/1‑150 (MQ=255)
                                                                                                                                                    cGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATAtt   >  2:1981455/1‑150 (MQ=255)
                                                                                                                                                     gAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGGTGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTg  >  2:1846009/1‑150 (MQ=255)
                                                                                                                                                     gAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTg  >  2:592836/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCCAGTCCCGAAGGACTGGAAGGCTCAATCGATCAAATCAATCAGAGGATGTGCACCTGCATTTCCGGGTCTTTGCGATCCAGATAGTGGATGGACTGAATGCGGCGAATGGTGCGTGACTTGCCGCGGATCAGCAGCGTTTCGGTAGTCGCGATATTG  >  NC_000913/4114388‑4114686

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