Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,899,757 G→C 7.7% intergenic (‑179/+55) PP_5172 ← / ← PP_5173 hypothetical protein/RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,899,7570GC 7.7% 463.5 / 24.3 143intergenic (‑179/+55)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (65/65);  new base C (4/7);  total (69/74)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 143 total aligned reads displayed.

CAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCC  >  NC_002947/5899609‑5899903
                                                                                                                                                    |                                                                                                                                                   
cAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGa                                                                                                                                                    <  1:3433180/150‑1 (MQ=255)
 agagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAt                                                                                                                                                   <  2:1462798/150‑1 (MQ=255)
  gagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATg                                                                                                                                                  <  1:1318668/150‑1 (MQ=255)
  gagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATg                                                                                                                                                  <  1:3776126/150‑1 (MQ=255)
  gagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATg                                                                                                                                                  <  1:3476443/150‑1 (MQ=255)
   agGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGc                                                                                                                                                 <  2:4105023/150‑1 (MQ=255)
    ggTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCa                                                                                                                                                <  2:3587426/150‑1 (MQ=255)
    ggTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCa                                                                                                                                                <  2:3587400/150‑1 (MQ=255)
        tGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGTTTAACGCGGATACGGTGGTGGCTGCACCaa                                                                                                                                            >  1:2268852/1‑149 (MQ=255)
          aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAtag                                                                                                                                          >  1:2521090/1‑148 (MQ=255)
          aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAAAACGGTGGTCGATGCACCAACg                                                                                                                                          >  1:2642031/1‑150 (MQ=255)
            ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                        >  1:724371/1‑150 (MQ=255)
            ccggTTTTCCGCCTGGCTTGACCCAAATATTAGAGAGGCGGTATTATACTTAAAAATACTGTTTAACAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                        <  1:3491722/148‑1 (MQ=255)
                 tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                   <  1:3647458/150‑1 (MQ=37)
                 tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   <  2:120130/150‑1 (MQ=255)
                 tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   <  2:2033380/150‑1 (MQ=255)
                 tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   <  2:3395870/150‑1 (MQ=255)
                 tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   <  1:2932741/150‑1 (MQ=255)
                  ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                                  <  1:3486963/150‑1 (MQ=37)
                  ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGATGGTGCATCGACCACCGTATTcgcg                                                                                                                                  <  2:610230/150‑1 (MQ=37)
                  ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                                  <  1:3311276/150‑1 (MQ=25)
                   tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                 <  1:3272683/150‑1 (MQ=25)
                   tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  1:203878/150‑1 (MQ=255)
                   tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  2:4103664/150‑1 (MQ=255)
                   tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  1:4356339/150‑1 (MQ=255)
                   tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  1:3107138/150‑1 (MQ=255)
                   tGCGCCTGGCTGGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  2:2120419/150‑1 (MQ=255)
                    gcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGTTACCCGCGAATACGGGGGTCGCTGCACCATCGCATTCGCggg                                                                                                                                >  2:2851780/1‑150 (MQ=255)
                      gcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTa                                                                                                                              <  1:4024821/150‑1 (MQ=255)
                      gcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTAGCGGGTa                                                                                                                              >  1:2574426/1‑150 (MQ=255)
                        cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTaaa                                                                                                                            <  1:24872/150‑1 (MQ=255)
                           gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATATGGTGGTCGATGCACCATCGCATTCGCCGGTAAAccc                                                                                                                         >  1:107310/1‑150 (MQ=255)
                            cTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 >  2:1713720/1‑141 (MQ=255)
                            cTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                 <  1:1713720/141‑1 (MQ=255)
                             ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAAAACCGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                       >  2:550739/1‑150 (MQ=255)
                                  ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                  >  2:1472755/1‑150 (MQ=255)
                                    cAATAGTAGTGAGGCGTGTTTATTATTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAATACAACACGCCTCTGGTTGGAGCGGGTTTACCCGAGAATAGGGTGGTCGATGCACCAACGCATTCGCGGGTAAACCCGCTCCTAcc                                                                                                                <  2:2470568/150‑1 (MQ=38)
                                     aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCa                                                                                                               >  1:1852309/1‑150 (MQ=255)
                                       tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                            >  2:4305762/1‑150 (MQ=255)
                                         gtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgggg                                                                                                          >  2:2152108/1‑150 (MQ=255)
                                          tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTACACCCGCTCCTACCAGGGGc                                                                                                         >  2:3281464/1‑150 (MQ=255)
                                             tgcggCGGGATTATTCTTAAAAATACTGTTTAAGATTTTTTATATAACCATTTAAGAGAACAGGCCTCTGGTCGGAGCGGGTTTACCCGCGAATCGGGTGGTCGATGCACCATCGCATTCGCGGGAAAACCCGCTCCTACCAGGGGCCat                                                                                                      <  2:2102865/147‑1 (MQ=38)
                                             tGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGGTTCCCGCGCAAAAGGTGGGCGATTCACCATAGCATTCGCGGGTAAACCCGCGCCGACCACGGGGccg                                                                                                       >  1:616753/1‑149 (MQ=255)
                                                          ttCTTAAAAATACTGTTTAAGATTTTTTATATAACCTTTTAAGAGACCACGCCTGTGGTTGGCGGGGGTTTACCCGCGAATAGGGTGGTCGACGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAgg                                                                                         <  2:1225489/150‑1 (MQ=38)
                                                             ttAAAATTACTGTTTAAGAATTTTTATATAACAATTTAAGAGAACACGCATCTTGTAGTAGCTGGTTTAAACGATAATAAGGTGGTAGATGCACCATCGAATTCGCGGGTAAAACCGCTCGTACAAGGGGGCATATTGGGCTTAAGGCgt                                                                                      >  2:3888604/1‑150 (MQ=255)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTAACCCGAAAAAGGGGGTCGATTCACCAACGCAGTCGCGGGGAAACCCCCACCTACCCGGGGCCCTATGGGGCGTCAGGCgg                                                                                      >  1:2404783/1‑149 (MQ=255)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                      <  2:3124750/150‑1 (MQ=255)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCGTCAGGCgt                                                                                      >  1:1862343/1‑150 (MQ=255)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                      <  2:2770154/150‑1 (MQ=255)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAAAACGGTGGTCGATGCACCCTCGCATTCGCGGGTAAACCCGCACCGACCAGGGGCCAGATGGGCCGTCACGCgt                                                                                      >  2:2454633/1‑150 (MQ=38)
                                                             ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACCGTGGTCGATGCACCAACGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                      >  1:1862340/1‑150 (MQ=255)
                                                                     aCTGTTTAAGAATTTTTATATAACCATTTACCAGACCACGCCTCTGGTCGGCGCGGGTTTACCCGCGAATAGGGTGGTCGATGCACAATCGCATTCGCGGTTAACCCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGc                                                                              <  1:2179242/150‑1 (MQ=38)
                                                                     aCTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTCGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGc                                                                              <  1:330183/150‑1 (MQ=255)
                                                                     aCTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGc                                                                              <  2:2995885/150‑1 (MQ=255)
                                                                            aatattttttATATAACAATTTAAGACACCACGCCTCTGGTAGGCGCGGTTTTACCCGCGAATAGGGTGGTCGATGCACCATCGCATTCGGGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                       <  1:291630/145‑1 (MQ=38)
                                                                               aaTTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCAACCAGGGGCCCGCTGGGCCGTCAGGCGTGTTTCTGCGGCGCGCGGa                                                                    >  2:2803179/1‑150 (MQ=38)
                                                                                 tttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACCGTGGTCGATGCCCCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCGTCAGGCGGGTTTCTgcggctcgcggctc                                                                  >  2:513029/1‑150 (MQ=38)
                                                                                 tttttATATAACATTTTAAGAGACCCCGCCTCTGGTAGGAGCGGTTTTACCCGCGAATACGGTGGTCGATGTACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATc                                                                  <  1:1403478/150‑1 (MQ=255)
                                                                                     tatataACCATTTCAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTACACCCGCTCCTACCAGGGGCCAGATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                              <  2:4502601/150‑1 (MQ=255)
                                                                                         taACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATGGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCa                                                          <  2:2341409/150‑1 (MQ=255)
                                                                                         taACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACCGTGGGCGATGCACCATCGCATTAGCGGGTAAAACCGCTCCTACCAGGGGCCATAAGGGCCGTCCGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCa                                                          >  1:4054938/1‑150 (MQ=38)
                                                                                          cacattttAACAGACCACGCCTGTGGTTGGCGGGGGTTTACCCGCTATTACGGTGGTCTAGCCACCATCGCATGCGGGGCAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                         <  1:550739/145‑1 (MQ=255)
                                                                                                 tAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGt                                                  >  1:4102521/1‑150 (MQ=255)
                                                                                                        acccGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgcgc                                           <  2:2268852/147‑1 (MQ=255)
                                                                                                         cacGCCTCTGGTAGGAGCGGGTTTACCCGCGCATACGGTGGTCGGTGCACCATCGCAGTCGCGGGTAAACCCGCTCCTACCAGGGGCCCTATGGGCCGTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCCGGCCACGTACAGCGcc                                          >  2:262133/1‑150 (MQ=38)
                                                                                                         cacGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCCTCTGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCCCGTACAGCGcc                                          >  2:262691/1‑150 (MQ=255)
                                                                                                              ctctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAGCGCCTTCGCGGGTAAACCCGCTCCGACCAGGGGCCATATGGGCCTTCAGGCGGGGTTCTGCGGCTCGCGGGTCTTGGACCCGGCCCCGTACAGCGCCGGCAg                                     >  2:4233601/1‑150 (MQ=37)
                                                                                                              ccctGGTAGGAGCGGGTTTACCCGCGAATGCGATGGGGCATCGACCACCGGATTCGCGGGTAAACCCGCTCCTCCCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACCGCGCCGGCAg                                     >  2:4233669/3‑150 (MQ=18)
                                                                                                              ccctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   >  2:1932444/3‑57 (MQ=255)
                                                                                                              ccctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                   <  1:1932444/55‑1 (MQ=255)
                                                                                                                 tGGTCGGTGGGGGTTTACCGGCGATTAGGGTGGACTATGCACATTCGAATCGGGGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                  <  1:3036890/150‑1 (MQ=255)
                                                                                                                 tGGTAGGAGCGGGTTTACCCCCCCAAAAGGTTGTCGCTGCCCCAACGCAATAGCGGGGAAACCCGCGCCGACCCCGGGCCCCATGGGCCGTCAGGCGGGTTTCAGCGGCTCGCGCAACTTGTACCACGCCCCCAACAGCCCCCGCAGGaa                                  >  2:1917501/1‑150 (MQ=255)
                                                                                                                   gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAca                                <  2:2100378/150‑1 (MQ=255)
                                                                                                                     aGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagc                              >  2:688836/1‑150 (MQ=255)
                                                                                                                        cgCGGGTTTCCCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                           <  2:2016223/149‑1 (MQ=255)
                                                                                                                         gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGTGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCg                          >  1:3237968/1‑150 (MQ=255)
                                                                                                                            ggTTGACCCGCGAATACGGTGGTCGAGGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCa                       <  2:114634/150‑1 (MQ=255)
                                                                                                                             gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAg                      >  2:1703692/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  2:1614810/1‑150 (MQ=25)
                                                                                                                                   ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  2:857671/1‑150 (MQ=25)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGa                                   >  2:2418745/1‑131 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGa                                   <  1:2418745/131‑1 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGACCAGCAGCGTCAGCAAGGt                >  2:3887390/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  1:4090454/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  1:554799/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  2:1034014/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  2:3887373/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  1:1213081/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGGTGGTCGATGCACCATCGAATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                >  2:1033848/1‑150 (MQ=255)
                                                                                                                                    cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc               <  1:1472755/150‑1 (MQ=255)
                                                                                                                                     gcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg              >  1:2754113/1‑150 (MQ=37)
                                                                                                                                     gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg              >  2:768494/1‑150 (MQ=255)
                                                                                                                                     gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg              >  1:3101452/1‑150 (MQ=255)
                                                                                                                                         aTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGTCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc          >  1:3445552/1‑150 (MQ=255)
                                                                                                                                            cGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat       >  2:4419560/1‑150 (MQ=255)
                                                                                                                                              gtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGc                                                                                                <  2:1716590/59‑1 (MQ=255)
                                                                                                                                              gtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGc                                                                                                >  1:1716590/1‑59 (MQ=255)
                                                                                                                                               tggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatgat    >  2:2017636/1‑150 (MQ=255)
                                                                                                                                               tggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatgat    >  2:1588451/1‑150 (MQ=255)
                                                                                                                                                 gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATcc  >  1:3683017/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCC  >  NC_002947/5899609‑5899903

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