Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,064 T→A 5.1% K351* (AAA→TAA)  dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,0640TA 5.1% 724.9 / 21.0 196K351* (AAA→TAA) dhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base T (91/95);  new base A (4/6);  total (95/101)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 208 total aligned reads displayed.

AGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247916‑1248212
                                                                                                                                                    |                                                                                                                                                    
         gTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             <  4:1619056/149‑2 (MQ=255)
         gTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTACGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             <  3:92089/149‑2 (MQ=255)
          ttGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTa                                                                                                                                                     >  1:289982/1‑140 (MQ=255)
           tGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgct                                                                                                                                           <  4:444379/149‑1 (MQ=255)
             gCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgctgc                                                                                                                                         >  3:1654451/1‑149 (MQ=255)
                  ttAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTaa                                                                                                                                                    <  4:510221/133‑1 (MQ=255)
                    aaaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgc                                                                                                                                            <  1:256304/139‑1 (MQ=255)
                     aaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      <  4:454691/128‑1 (MQ=255)
                     aaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      >  3:454691/1‑128 (MQ=255)
                     aaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCata                                                                                                                                 >  3:480803/1‑149 (MQ=255)
                           gggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTc                                                                                                                              >  1:466884/1‑146 (MQ=255)
                           gggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGc                                                                                                                           <  3:904549/149‑1 (MQ=255)
                                  cttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCTTGCaa                                                                                                                    >  4:195596/1‑149 (MQ=255)
                                  cttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCaa                                                                                                                    <  2:751535/149‑1 (MQ=255)
                                   ttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgctg                                                                                                                                          >  2:1351920/1‑126 (MQ=255)
                                   ttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   >  2:1605813/1‑149 (MQ=255)
                                       ttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTc                                                                                                               <  4:1833515/149‑1 (MQ=255)
                                         cGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  2:1659604/149‑1 (MQ=255)
                                         cGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  2:1145452/149‑1 (MQ=255)
                                         cGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  1:1353914/149‑1 (MQ=255)
                                           tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   >  2:868138/1‑141 (MQ=255)
                                           tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   <  1:868138/141‑1 (MQ=255)
                                           tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGc                                                                                                           >  2:874491/1‑149 (MQ=255)
                                                gTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             >  4:908995/1‑142 (MQ=255)
                                                  cAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGa                                                                                                    <  3:153263/149‑1 (MQ=255)
                                                  cAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGa                                                                                                    <  1:1481144/149‑1 (MQ=255)
                                                   aGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAg                                                                                                   <  4:1666972/149‑1 (MQ=255)
                                                    ggTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATg                                                                                                                  <  1:1610606/133‑1 (MQ=255)
                                                        gACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGc                                                                                              <  1:874491/149‑1 (MQ=255)
                                                          ccAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCtt                                                                                            <  2:403689/149‑1 (MQ=255)
                                                           cAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  4:851754/131‑1 (MQ=255)
                                                           cAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTg                                                                                           >  4:1193613/1‑149 (MQ=255)
                                                            aGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGt                                                                                                                 >  4:1414295/1‑126 (MQ=255)
                                                                   gATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTc                                                                                                                              >  2:1387990/1‑106 (MQ=255)
                                                                   gATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgcc                                                                                   >  4:1364603/1‑149 (MQ=255)
                                                                   gATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgcc                                                                                   >  2:269872/1‑149 (MQ=255)
                                                                     tGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcag                                                                                 >  2:1179314/1‑149 (MQ=255)
                                                                       cAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               >  4:249995/1‑149 (MQ=255)
                                                                         gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      <  4:731251/124‑1 (MQ=255)
                                                                         gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGt                                                                                                  >  1:1026513/1‑128 (MQ=255)
                                                                         gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGt                                                                                                  <  2:1026513/128‑1 (MQ=255)
                                                                          aagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      >  1:524412/1‑75 (MQ=255)
                                                                              tCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               <  1:134924/142‑1 (MQ=255)
                                                                              tCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               >  2:1321753/1‑142 (MQ=255)
                                                                                gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      <  2:1793239/149‑1 (MQ=255)
                                                                                 tCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTc                                                                                                >  1:1295864/1‑122 (MQ=255)
                                                                                    acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  4:1727326/149‑1 (MQ=255)
                                                                                    acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  >  1:1249000/1‑149 (MQ=255)
                                                                                    acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  2:1088276/149‑1 (MQ=255)
                                                                                             taGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTc                                                                           <  3:1373946/131‑1 (MQ=255)
                                                                                               gCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAAc                                                              <  2:1711615/142‑1 (MQ=255)
                                                                                               gCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAAc                                                              >  1:1711615/1‑142 (MQ=255)
                                                                                                cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgcc                                                                                      >  3:705123/1‑117 (MQ=255)
                                                                                                cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATg                                                      <  2:1519807/149‑1 (MQ=255)
                                                                                                     ttGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgct                                                                                                                                           <  1:936294/59‑1 (MQ=255)
                                                                                                        tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  >  4:1466446/1‑129 (MQ=255)
                                                                                                        tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  3:1466446/129‑1 (MQ=255)
                                                                                                        tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTaa                                                               <  2:1423136/132‑1 (MQ=255)
                                                                                                         aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      >  1:1077790/1‑92 (MQ=255)
                                                                                                         aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGtt                                                                                                 >  4:1426536/1‑97 (MQ=255)
                                                                                                         aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  >  2:1050734/1‑128 (MQ=255)
                                                                                                         aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  1:1050734/128‑1 (MQ=255)
                                                                                                          gATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGaa                                            <  3:249995/149‑1 (MQ=255)
                                                                                                            tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  2:908871/82‑1 (MQ=255)
                                                                                                            tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             >  1:908871/1‑82 (MQ=255)
                                                                                                             atTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcat                                                                   <  2:473402/123‑1 (MQ=255)
                                                                                                               tcatcaTCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAt                                       >  2:1809193/1‑149 (MQ=255)
                                                                                                               tcatcaTCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAt                                       <  3:1863367/149‑1 (MQ=255)
                                                                                                                catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAtt                                      >  3:19915/1‑149 (MQ=255)
                                                                                                                  tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGc                                    <  3:1755750/149‑1 (MQ=255)
                                                                                                                      ccAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCa                                <  2:1249000/149‑1 (MQ=255)
                                                                                                                          ttGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCg                                   >  1:1547189/1‑142 (MQ=255)
                                                                                                                           tGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATc                           <  3:340218/149‑1 (MQ=255)
                                                                                                                              tGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGt                                                                                                  <  1:1639111/75‑1 (MQ=255)
                                                                                                                              tGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTc                        <  1:1809193/149‑1 (MQ=255)
                                                                                                                              tGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTc                        <  2:24542/149‑1 (MQ=255)
                                                                                                                                gTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      <  2:100194/149‑1 (MQ=255)
                                                                                                                                gTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      >  1:1140359/1‑149 (MQ=255)
                                                                                                                                gTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      >  1:100194/1‑149 (MQ=255)
                                                                                                                                 tATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcat                                                                   <  4:1060971/103‑1 (MQ=255)
                                                                                                                                 tATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcat                                                                   >  3:1060971/1‑103 (MQ=255)
                                                                                                                                       tGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTa                                                        <  2:944501/108‑1 (MQ=255)
                                                                                                                                        ggCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTg              >  1:1218106/1‑149 (MQ=255)
                                                                                                                                          gcaaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGc            >  4:608477/3‑149 (MQ=255)
                                                                                                                                          cTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGc            <  3:83873/149‑1 (MQ=255)
                                                                                                                                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCt           >  3:986050/1‑149 (MQ=255)
                                                                                                                                           tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCt           >  3:1128672/1‑149 (MQ=255)
                                                                                                                                             aGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTGaa                      >  2:545364/1‑136 (MQ=255)
                                                                                                                                              gTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          <  1:202426/83‑1 (MQ=255)
                                                                                                                                              gTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          >  2:202426/1‑83 (MQ=255)
                                                                                                                                              gTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCt        <  1:749902/149‑1 (MQ=255)
                                                                                                                                               ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               <  1:1400283/77‑1 (MQ=255)
                                                                                                                                               ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAc                                                 >  3:1420511/1‑107 (MQ=255)
                                                                                                                                               ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGa                                        <  1:659954/116‑1 (MQ=255)
                                                                                                                                               ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGa                                        >  2:659954/1‑116 (MQ=255)
                                                                                                                                                tCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      >  2:1633677/1‑133 (MQ=255)
                                                                                                                                                  tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGGAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGTTTGCGGCAATATCGTCAAGCCCGCTGGCTACTGCtt    >  3:1259514/1‑149 (MQ=255)
                                                                                                                                                  tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  4:378008/1‑149 (MQ=255)
                                                                                                                                                  tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  1:361683/1‑149 (MQ=255)
                                                                                                                                                    tAAGAACTTGCTGCCAGGTATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGTTGACCAGAAAAGATTGCGGCAATAGCGTCAAGCCCGCTGGCTTCTGCTTTc  <  2:703630/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247916‑1248212

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