Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,144 G→C 5.4% intergenic (+17/‑94) aphA → / → yjbQ acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,1440GC 5.4% 492.4 / 15.6 147intergenic (+17/‑94)aphA/yjbQacid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (70/69);  new base C (4/4);  total (74/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 148 total aligned reads displayed.

AATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGT  >  NC_000913/4269996‑4270293
                                                                                                                                                    |                                                                                                                                                     
 aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                                                                     <  2:2160916/150‑1 (MQ=255)
   gATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                                                                                                                   <  1:1541300/150‑1 (MQ=255)
     tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                                                                                                                 >  1:1991448/1‑150 (MQ=255)
     tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                                                                                                                 <  2:594701/150‑1 (MQ=255)
        tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                                                                                                              >  2:887950/1‑150 (MQ=255)
        tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTcc                                                                                                                                              >  1:1902215/1‑150 (MQ=255)
        tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTcc                                                                                                                                              >  2:1349917/1‑150 (MQ=255)
         aCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTccc                                                                                                                                             >  1:1218769/1‑150 (MQ=255)
          ccgccgCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                                                                            <  2:892902/150‑1 (MQ=255)
            gccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGaa                                                                                                                                          >  1:516620/1‑150 (MQ=255)
             ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                                                                         <  2:1628095/150‑1 (MQ=255)
             ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                                                                         <  1:113388/150‑1 (MQ=255)
              cgcACGCGATGTCGGCGCTCGTGGGATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATg                                                                                                                                        <  2:2115537/150‑1 (MQ=255)
               gcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGc                                                                                                                                       <  2:410639/150‑1 (MQ=255)
               gcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctctg                                                                                                                                       <  2:2215963/150‑6 (MQ=255)
               gcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctctg                                                                                                                                       <  2:487125/150‑6 (MQ=255)
                cACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCt                                                                                                                                      >  2:760515/1‑150 (MQ=255)
                 aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTg                                                                                                                                     <  1:878343/150‑1 (MQ=255)
                  cgcgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGt                                                                                                                                    >  1:690771/1‑150 (MQ=255)
                   gcgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtt                                                                                                                                   >  1:1253019/1‑150 (MQ=255)
                      aTGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGAGCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttt                                                                                                                                <  1:975758/150‑1 (MQ=255)
                       tGTCGGCGCTCGTGGTCTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTGGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                                                                               <  1:1683781/150‑1 (MQ=255)
                       tGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                                                                               >  2:820138/1‑150 (MQ=255)
                        gTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttttt                                                                                                                              <  2:1991448/150‑1 (MQ=255)
                          cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                            >  2:229282/1‑150 (MQ=255)
                          cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                            >  2:631588/1‑150 (MQ=255)
                          cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                            <  1:797574/150‑1 (MQ=255)
                          cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                            >  1:1028640/1‑150 (MQ=255)
                          cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                            >  1:1039931/1‑150 (MQ=255)
                           ggCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                                                                           <  1:1250248/150‑1 (MQ=255)
                           ggCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                                                                           >  2:1185015/1‑150 (MQ=255)
                                tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                      <  1:1034565/150‑1 (MQ=255)
                                tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                      <  1:862933/150‑1 (MQ=255)
                                 cGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                                                                     >  2:17230/1‑150 (MQ=255)
                                   tGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCtt                                                                                                                   <  2:1954593/150‑1 (MQ=255)
                                      tATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                >  2:1123191/1‑150 (MQ=255)
                                      tATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                                                                                <  1:2202835/150‑1 (MQ=255)
                                       aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAAGTCAGAATACTGCCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                                                                               >  2:1284172/1‑150 (MQ=255)
                                        tCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATcc                                                                                                              >  2:809748/1‑150 (MQ=255)
                                         ccGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                                                                             >  1:52489/1‑150 (MQ=255)
                                             aTTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                                                                         <  2:834939/150‑1 (MQ=255)
                                                    gcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCtt                                                                                                  <  1:2112414/150‑1 (MQ=255)
                                                    gcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCtt                                                                                                  <  1:1423728/150‑1 (MQ=255)
                                                     cgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtg                                                                                                 >  1:1767838/1‑150 (MQ=255)
                                                       ccTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgct                                                                                               >  2:1680558/1‑150 (MQ=255)
                                                       ccTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCACGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgct                                                                                               >  2:359948/1‑150 (MQ=255)
                                                         tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                                                                             <  1:889530/150‑1 (MQ=255)
                                                         tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACATTTATCCTTTCGCTGTCTtgctgc                                                                                             >  2:187387/1‑150 (MQ=255)
                                                          ccAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCa                                                                                            >  2:1038281/1‑150 (MQ=255)
                                                             aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc                                                                                         <  2:969212/150‑1 (MQ=255)
                                                             aCTCTACCTACAAACCCATGCCACAAGCGGGAGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAGCCTTTCGCTGTCTTGCTGCAAAc                                                                                         <  1:75468/150‑1 (MQ=255)
                                                              ctctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt                                                                                        <  1:17230/150‑1 (MQ=255)
                                                               tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTg                                                                                       <  2:1785167/150‑1 (MQ=255)
                                                                ctacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATAACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa                                                                                      >  1:2159692/1‑150 (MQ=255)
                                                                     tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                                 <  2:904905/150‑1 (MQ=255)
                                                                     tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                                 >  2:16827/1‑150 (MQ=255)
                                                                     tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                                 <  1:2015377/150‑1 (MQ=255)
                                                                       cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGTGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                                                                               <  2:575386/150‑1 (MQ=255)
                                                                       cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                                                                               >  2:1636941/1‑150 (MQ=255)
                                                                             cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTATTAGAGtttt                                                                         >  1:428597/1‑150 (MQ=255)
                                                                              ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                        <  2:1405331/150‑1 (MQ=255)
                                                                              ttGACACAAGCGTGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                        <  2:453485/150‑1 (MQ=255)
                                                                               tGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTAt                                                                       >  2:1830541/1‑150 (MQ=255)
                                                                                  cacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                                    >  2:904406/1‑150 (MQ=255)
                                                                                   acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg                                                                   <  1:413302/150‑1 (MQ=255)
                                                                                       gCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagc                                                               <  1:1234825/150‑1 (MQ=255)
                                                                                          ggTGCGTGTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGGCGTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcagc                                                            <  1:1121412/150‑1 (MQ=255)
                                                                                           gTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAgcagca                                                           <  2:926369/150‑1 (MQ=255)
                                                                                             gCGTTTGGTGAAGAGTTTATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTATATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGcaca                                                         <  1:134142/150‑1 (MQ=255)
                                                                                              cGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAt                                                        >  2:685722/1‑150 (MQ=255)
                                                                                                   ggTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGt                                                   >  1:2102145/1‑150 (MQ=255)
                                                                                                    gTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCATTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTa                                                  >  2:127946/1‑150 (MQ=255)
                                                                                                      gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATc                                                >  2:205160/1‑150 (MQ=255)
                                                                                                       aaGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCa                                               <  2:1072442/150‑1 (MQ=255)
                                                                                                          agGTGATCGTCAATTCAGAATCCTGACAGAGCGTGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGACAGTTTTATCAAGGAGCAGCACATGTGGTATCaaaa                                            <  1:175109/150‑1 (MQ=255)
                                                                                                           ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTCTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                           <  2:516620/150‑1 (MQ=255)
                                                                                                           ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGATCTCCAGAATGCTGTTTTTTTAATCAAACCTTTATCCTATCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                           >  2:767103/1‑150 (MQ=255)
                                                                                                                       ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCtt                               <  1:1123191/150‑1 (MQ=255)
                                                                                                                       ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCtt                               >  1:1181987/1‑150 (MQ=255)
                                                                                                                         cAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAg                             >  2:236128/1‑150 (MQ=255)
                                                                                                                         cAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAg                             <  2:1723945/150‑1 (MQ=255)
                                                                                                                          aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGAATCAAAAGACGCTCACGCTTAgc                            <  1:902995/150‑1 (MQ=255)
                                                                                                                                cTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCaaa                      >  2:1960616/1‑150 (MQ=255)
                                                                                                                                cTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCaaa                      >  2:2053114/1‑150 (MQ=255)
                                                                                                                                   aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtct                   <  1:631588/150‑1 (MQ=255)
                                                                                                                                   aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtct                   >  1:233800/1‑150 (MQ=255)
                                                                                                                                   aCAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtct                   <  1:229282/150‑1 (MQ=255)
                                                                                                                                    cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtctc                  >  1:638649/1‑150 (MQ=255)
                                                                                                                                    cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtctc                  >  1:1180188/1‑150 (MQ=255)
                                                                                                                                     agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg                 >  1:1023397/1‑150 (MQ=255)
                                                                                                                                      gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt                >  1:1648522/1‑150 (MQ=255)
                                                                                                                                          gggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt            >  2:255002/1‑150 (MQ=255)
                                                                                                                                          gggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGt            <  1:447980/150‑1 (MQ=255)
                                                                                                                                          gggAGAGCATCACGCTCTACCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCgtgtgt            >  1:1542022/1‑150 (MQ=255)
                                                                                                                                           ggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTACGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGtt           >  1:679240/1‑150 (MQ=255)
                                                                                                                                               agCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCAt       >  1:1973510/1‑150 (MQ=255)
                                                                                                                                                   tGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTgg   <  1:359948/150‑1 (MQ=255)
                                                                                                                                                   tGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTgg   <  1:1680558/150‑1 (MQ=255)
                                                                                                                                                   tGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTgg   >  1:1743298/1‑150 (MQ=255)
                                                                                                                                                    gATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTAAGCGCCAAATCTCGTGGGTTTCATCTGGt  >  2:493349/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGT  >  NC_000913/4269996‑4270293

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: