Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,064 T→A 8.6% 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 8.6% 309.3 / 16.8 93K351* (AAA→TAA) dhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base A (3/5);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

TAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247918‑1248212
                                                                                                                                                  |                                                                                                                                                    
tAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTaa                                                                                                                                                    <  4:599800/149‑1 (MQ=255)
   aaTAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGaa                                                                                                                                                 >  3:300946/1‑149 (MQ=255)
   aaTAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGaa                                                                                                                                                 <  3:202222/149‑1 (MQ=255)
        ttGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGaa                                                                                                                                                 >  4:216043/1‑144 (MQ=255)
        ttGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGaa                                                                                                                                                 <  3:216043/144‑1 (MQ=255)
        ttGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgc                                                                                                                                            <  4:300946/149‑1 (MQ=255)
        ttGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttag                                                                                                                                            <  4:290704/149‑3 (MQ=255)
         tGGCGAGTTAAACTGCGGGAGTTCTTCTCTCGTTTGGGTTAGGTAGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttagc                                                                                                                                           <  3:316777/149‑4 (MQ=255)
          ggCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGACATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgctg                                                                                                                                          <  1:84216/149‑1 (MQ=255)
               gTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAgg                                                                                                                                     >  1:226006/1‑149 (MQ=255)
               gTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAgg                                                                                                                                     <  2:278818/149‑1 (MQ=255)
                  aaaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACtt                                                                                                                                              <  1:130487/137‑1 (MQ=255)
                  aaaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACtt                                                                                                                                              >  2:130487/1‑137 (MQ=255)
                  aaaCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCat                                                                                                                                  >  3:446118/1‑149 (MQ=255)
                    aCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCatat                                                                                                                                >  4:665900/1‑149 (MQ=255)
                          ggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGATCTTGCTGCCAGGCATATTCTGcc                                                                                                                          <  4:421061/149‑1 (MQ=255)
                           gAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGcc                                                                                                                          >  4:216165/1‑148 (MQ=255)
                           gAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGcc                                                                                                                          <  3:216165/148‑1 (MQ=255)
                                cttcatTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCaa                                                                                                                    >  4:519485/1‑149 (MQ=255)
                                 ttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   >  1:28285/1‑149 (MQ=255)
                                    tttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGtt                                                                                                                >  4:541209/1‑149 (MQ=255)
                                        gTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATg                                                                                                            <  3:468377/149‑1 (MQ=255)
                                         tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCa                                                                                                                                       >  1:453604/1‑121 (MQ=255)
                                         tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCa                                                                                                                                       <  2:453604/121‑1 (MQ=255)
                                         tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGc                                                                                                           <  1:457752/149‑1 (MQ=255)
                                         tttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGc                                                                                                           >  3:221413/1‑149 (MQ=255)
                                               tCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAgg                                                                                                     >  4:18438/1‑149 (MQ=255)
                                                cAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGa                                                                                                    <  2:302531/149‑1 (MQ=255)
                                                   gTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGtt                                                                                                 >  1:288604/1‑149 (MQ=255)
                                                   gTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGtt                                                                                                 >  1:142370/1‑149 (MQ=255)
                                                   gTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGtt                                                                                                 <  2:618214/149‑1 (MQ=255)
                                                     ggACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCg                                                                                               >  1:443443/1‑149 (MQ=255)
                                                      gACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGc                                                                                              <  1:548868/149‑1 (MQ=255)
                                                       aCCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCt                                                                                             <  3:476858/149‑1 (MQ=255)
                                                         cAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTg                                                                                           >  3:547790/1‑149 (MQ=255)
                                                                 gATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTACTGGCTGAGTTCTTTAAGAACTTGCTGCCAGGCATATTCAGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgcc                                                                                   >  4:448347/1‑149 (MQ=255)
                                                                  aTGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCca                                                                                  <  3:541209/149‑1 (MQ=255)
                                                                     cAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               >  2:237850/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAGGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  1:636072/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  1:189752/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  3:169718/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  2:661318/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             >  2:73879/1‑84 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             <  1:73879/85‑2 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTc                                                                                                <  1:490159/127‑1 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTc                                                                                                >  2:490159/1‑127 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          <  1:226121/149‑1 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          >  4:498978/1‑149 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          >  2:650724/1‑149 (MQ=255)
                                                                          gaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTcc                                                                          >  4:653057/1‑149 (MQ=255)
                                                                             cGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGAAGTTCCGGa                                                                       <  2:522798/149‑1 (MQ=255)
                                                                              gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      >  1:563080/1‑149 (MQ=255)
                                                                                  acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  1:272032/149‑1 (MQ=255)
                                                                                    aTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      <  2:80755/111‑1 (MQ=255)
                                                                                    aTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      >  1:80755/1‑111 (MQ=255)
                                                                                    aTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTa                                                                >  1:476825/1‑149 (MQ=255)
                                                                                       aaTATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACa                                                             <  2:288604/149‑1 (MQ=255)
                                                                                        atatAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAg                                                            <  2:476825/149‑1 (MQ=255)
                                                                                           taGCGTGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtgt                                                         <  4:446118/149‑1 (MQ=255)
                                                                                              cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATg                                                      <  3:653057/149‑1 (MQ=255)
                                                                                                aGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGt                                                    <  3:57561/149‑1 (MQ=255)
                                                                                                 gCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTg                                                   <  1:650724/149‑1 (MQ=255)
                                                                                                   ttGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAc                                                 <  3:651587/149‑1 (MQ=255)
                                                                                                      tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGt                                                                                                                 <  1:278009/82‑1 (MQ=255)
                                                                                                      tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGt                                                                                                                 >  2:278009/1‑82 (MQ=255)
                                                                                                       aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      <  3:208868/44‑1 (MQ=255)
                                                                                                       aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      >  4:208868/1‑44 (MQ=255)
                                                                                                       aGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGa                                             >  3:228527/1‑149 (MQ=255)
                                                                                                          tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGaaaa                                          <  3:665900/149‑1 (MQ=255)
                                                                                                           atTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      <  2:130168/120‑1 (MQ=255)
                                                                                                           atTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      >  1:130168/1‑120 (MQ=255)
                                                                                                           atTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAg                                         >  3:186252/1‑149 (MQ=255)
                                                                                                            ttCATCAGCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAGTGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGa                                        >  2:372769/1‑149 (MQ=255)
                                                                                                              catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAtt                                      >  2:381968/1‑149 (MQ=255)
                                                                                                                tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTg                                                   >  4:597341/1‑134 (MQ=255)
                                                                                                                tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTg                                                   <  3:597341/134‑1 (MQ=255)
                                                                                                                  atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCgg                                  <  1:45772/149‑1 (MQ=255)
                                                                                                                  atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCgg                                  >  1:631040/1‑149 (MQ=255)
                                                                                                                  atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGATGACCAGAAAAGATTGCgg                                  >  1:52490/1‑149 (MQ=255)
                                                                                                                     cAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCaa                               >  3:688748/1‑149 (MQ=255)
                                                                                                                      aGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAat                              <  1:237850/149‑1 (MQ=255)
                                                                                                                      aGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAat                              <  1:661318/149‑1 (MQ=255)
                                                                                                                      aGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAat                              <  2:142370/149‑1 (MQ=255)
                                                                                                                           cTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGt                         >  2:463787/1‑149 (MQ=255)
                                                                                                                            tGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTc                        >  4:504065/1‑149 (MQ=255)
                                                                                                                                    cTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGc                >  1:202148/1‑149 (MQ=255)
                                                                                                                                       agcaaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTgg             >  4:461204/4‑149 (MQ=255)
                                                                                                                                         tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCt           >  2:151575/1‑149 (MQ=255)
                                                                                                                                         caaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAcc                                                >  3:665098/2‑112 (MQ=255)
                                                                                                                                         caaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAcc                                                <  4:665098/111‑1 (MQ=255)
                                                                                                                                           aGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTc         <  4:221413/149‑1 (MQ=255)
                                                                                                                                           aGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTc         <  2:28285/149‑1 (MQ=255)
                                                                                                                                           aGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTc         <  1:381968/149‑1 (MQ=255)
                                                                                                                                           aGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCTTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTc         <  3:461204/149‑1 (MQ=255)
                                                                                                                                             ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTg       >  2:667287/1‑149 (MQ=255)
                                                                                                                                                tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  1:418040/1‑149 (MQ=255)
                                                                                                                                                  tAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAGTGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTc  <  1:372769/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247918‑1248212

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