Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,681 C→G 15.6% V419V (GTG→GTC rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,6810CG15.6% 521.5 / 80.8 192V419V (GTG→GTCrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base C (89/73);  new base G (13/17);  total (102/90)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 192 total aligned reads displayed.

TACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCG  >  NC_002947/471533‑471828
                                                                                                                                                    |                                                                                                                                                   
 aCGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCAc                                                                                                                                                   >  2:5279963/1‑150 (MQ=255)
    gTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGc                                                                                                                                                >  1:2620376/1‑150 (MQ=255)
        aGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                                                                                                            >  2:1962717/1‑150 (MQ=255)
          cTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcat                                                                                                                                          <  1:4206843/150‑1 (MQ=255)
             gttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAg                                                                                                                                       >  2:1742702/1‑150 (MQ=255)
              ttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                                                                                                            >  2:475055/1‑144 (MQ=255)
               tgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                                                                                                     >  2:731387/1‑150 (MQ=255)
                  tCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGa                                                                                                                                  >  1:571170/1‑150 (MQ=255)
                    gTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATg                                                                                                                                >  1:1546532/1‑150 (MQ=255)
                     tttcGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGAtat                                                                                                                               >  2:1655452/3‑148 (MQ=255)
                     tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTTTTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                                                                                               <  2:4604970/150‑1 (MQ=255)
                      ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                              >  2:875315/1‑150 (MQ=255)
                      ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                              <  2:270461/150‑1 (MQ=255)
                      ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                              <  1:826929/150‑1 (MQ=255)
                      ctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGtt                                                                                                                              >  2:2329630/1‑150 (MQ=255)
                         gATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                                                                                                           >  1:4122234/1‑150 (MQ=255)
                             aTGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTcc                                                                                                                       <  2:4281284/150‑1 (MQ=255)
                              tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                                                                                                      <  1:2943244/150‑1 (MQ=255)
                                tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                                                                                                    >  1:389662/1‑150 (MQ=255)
                                  cACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGAt                                                                                                                  >  2:4430935/1‑150 (MQ=255)
                                   aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                                                                                 <  1:3642127/150‑1 (MQ=255)
                                   aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                                                                                 >  1:2378438/1‑150 (MQ=255)
                                   aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATc                                                                                                                 >  2:2406207/1‑150 (MQ=255)
                                    ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                                                                                >  2:2519252/1‑150 (MQ=255)
                                    ccGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCa                                                                                                                <  1:3856609/150‑1 (MQ=255)
                                     cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                                                                                               >  1:3184156/1‑150 (MQ=255)
                                        aaTGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTc                                                                                                            <  2:4324989/150‑1 (MQ=255)
                                          tacGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGa                                                                                                          <  2:2120921/148‑1 (MQ=255)
                                          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                          <  1:4029793/150‑1 (MQ=255)
                                          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                          <  2:1080631/150‑1 (MQ=255)
                                          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                          >  1:4551397/1‑150 (MQ=255)
                                          tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                                                                                                          <  2:2475398/150‑1 (MQ=255)
                                           gCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAg                                                                                                         <  2:2635398/150‑1 (MQ=255)
                                             ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                                                                                                       >  2:1775890/1‑150 (MQ=255)
                                              gATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGaa                                                                                                      <  2:942900/150‑1 (MQ=255)
                                              gATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGaa                                                                                                      <  1:731387/150‑1 (MQ=255)
                                               aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                                                                                     <  1:4400983/150‑1 (MQ=255)
                                               aTGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAAc                                                                                                     >  1:4708495/1‑150 (MQ=255)
                                                    gcgtgcCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                                <  1:2157765/150‑1 (MQ=255)
                                                         ccTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCAc                                                                                           >  2:5398720/1‑150 (MQ=255)
                                                             gTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                                                                                                            <  2:1239666/97‑1 (MQ=255)
                                                             gTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGtt                                                                                       >  1:1098123/1‑150 (MQ=255)
                                                             atCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGtt                                                                                       <  1:3353532/149‑1 (MQ=255)
                                                                ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTTGCCGCGTCGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAGGTTGCCTTCCTGGATCAGGTCGTGGACCTGCCGGCCTCGGTTGgt                                                                                    >  2:1414089/1‑150 (MQ=255)
                                                                ggCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgt                                                                                    >  1:4448453/1‑150 (MQ=255)
                                                                ggCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGgt                                                                                    >  2:5216818/1‑150 (MQ=255)
                                                                 gCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAGTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtg                                                                                   >  2:3497642/1‑150 (MQ=255)
                                                                   gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                 <  1:3467782/150‑1 (MQ=255)
                                                                   gATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTa                                                                                 >  2:5169306/1‑150 (MQ=255)
                                                                    aTCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTAc                                                                                >  1:5137688/1‑150 (MQ=255)
                                                                          cGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTg                                                                          <  1:2537292/150‑1 (MQ=255)
                                                                               gATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGAt                                                                     >  1:55200/1‑150 (MQ=255)
                                                                                   gCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGag                                                                 <  2:920114/150‑2 (MQ=255)
                                                                                    ccTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGaaa                                                                >  2:4247056/1‑150 (MQ=255)
                                                                                    ccTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGaga                                                                <  1:1433117/150‑3 (MQ=255)
                                                                                       gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                             >  1:2820951/1‑150 (MQ=255)
                                                                                       gACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCa                                                             >  1:4739836/1‑150 (MQ=255)
                                                                                       gACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGAtaa                                                             <  1:2661586/150‑3 (MQ=255)
                                                                                          ggATCCACCTGGTTGAATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                          <  1:1884146/150‑1 (MQ=255)
                                                                                          ggATCCACCAGGTGGCATAGGTCGAGAACTTGGAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCa                                                          <  2:3393797/150‑1 (MQ=255)
                                                                                           gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                         >  1:4513983/1‑150 (MQ=255)
                                                                                           gATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAg                                                         >  1:468180/1‑150 (MQ=255)
                                                                                             tCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGAc                                                       <  1:3564483/150‑1 (MQ=255)
                                                                                               caccaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGc                                                     <  2:392192/150‑1 (MQ=255)
                                                                                                accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGACATCACCAGACGCa                                                    >  2:3993499/1‑150 (MQ=255)
                                                                                                accaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCa                                                    >  2:3580962/1‑150 (MQ=255)
                                                                                                 ccaTGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAg                                                   <  2:1905996/150‑1 (MQ=255)
                                                                                                 ccaGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAg                                                   <  2:2183214/150‑1 (MQ=255)
                                                                                                   aGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTATTTCTTGGCGATGGAAATCACCAGACGCAGGt                                                 >  2:1285654/1‑150 (MQ=255)
                                                                                                     gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                               >  2:1328735/1‑150 (MQ=255)
                                                                                                     gTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTc                                               >  2:4080081/1‑150 (MQ=255)
                                                                                                         cATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCt                                           <  2:2850456/150‑1 (MQ=255)
                                                                                                           tAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCg                                         >  1:4807603/1‑150 (MQ=255)
                                                                                                               tCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCa                                     >  2:298188/1‑150 (MQ=255)
                                                                                                                cGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCAt                                    <  1:45901/150‑1 (MQ=255)
                                                                                                                  agaACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAttt                                  >  2:3872103/1‑150 (MQ=255)
                                                                                                                       ttGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCAtttctttc                             >  2:3076891/1‑150 (MQ=255)
                                                                                                                         gTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTtctttctt                           <  2:324741/150‑1 (MQ=255)
                                                                                                                         gTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTtccttctt                           <  1:1655452/150‑1 (MQ=255)
                                                                                                                            gCCGCGACGGTATTCGACCGTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                        <  2:5049857/150‑1 (MQ=255)
                                                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                        <  1:1221674/150‑1 (MQ=255)
                                                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                        <  1:226292/150‑1 (MQ=255)
                                                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                        <  1:5200110/150‑1 (MQ=255)
                                                                                                                            gCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGgc                        <  1:2472355/150‑1 (MQ=255)
                                                                                                                                 gACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGGTGTACTTCTTGGCAATCGAGATAACGAGACGAAGGTTGGCTTCCACCATTTCCTTCTTGgcgatgg                   >  2:1884534/1‑146 (MQ=255)
                                                                                                                                     gTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGc               >  1:1512802/1‑150 (MQ=255)
                                                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCt             >  2:953899/1‑150 (MQ=255)
                                                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCt             >  2:226037/1‑150 (MQ=255)
                                                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGACTTCCTGGATCAGGTCGAGGAACTGCAGGCCa                                                                                            <  1:1332995/71‑1 (MQ=255)
                                                                                                                                       aTTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGACTTCCTGGATCAGGTCGAGGAACTGCAGGCCa                                                                                            >  2:1332995/1‑71 (MQ=255)
                                                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                                <  1:4652760/66‑1 (MQ=255)
                                                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAg                                                                                                >  2:4652760/1‑66 (MQ=255)
                                                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt            >  2:3122064/1‑150 (MQ=255)
                                                                                                                                        ttCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCtt            <  1:1998446/150‑1 (MQ=255)
                                                                                                                                         tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg           >  1:4404624/1‑150 (MQ=255)
                                                                                                                                         tCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg           <  2:2620376/150‑1 (MQ=255)
                                                                                                                                           gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGACGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg           >  1:2315952/1‑148 (MQ=255)
                                                                                                                                           gAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGACGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTg           <  2:2315952/148‑1 (MQ=255)
                                                                                                                                            aaCTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCt        >  1:3797189/1‑150 (MQ=255)
                                                                                                                                                 gTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCAcc   >  2:2662692/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTACTTCTTGGCGATGGAAATCACCAGACGCAGGTTCGCCTCGACCATTTCTTTCTTGGCGCGGCGGGCCTTGGCTTCACCG  >  NC_002947/471533‑471828

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