Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,899,760 G→C 13.5% intergenic (‑182/+52) 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,7600GC13.5% 564.9 / 68.4 214intergenic (‑182/+52)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (88/93);  new base C (7/22);  total (99/115)
Fisher's exact test for biased strand distribution p-value = 1.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 214 total aligned reads displayed.

AGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCA  >  NC_002947/5899612‑5899908
                                                                                                                                                     |                                                                                                                                                    
    ttGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCa                                                                                                                                                 <  1:4603395/150‑1 (MQ=255)
      gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                               <  1:574207/150‑1 (MQ=255)
       aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCg                                                                                                                                              >  2:3138440/1‑150 (MQ=255)
         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATcaca                                                                                                                                            >  1:971708/1‑150 (MQ=255)
         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                            >  2:8061502/1‑150 (MQ=255)
         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                            <  1:2892846/150‑1 (MQ=255)
          gggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCCCGAATAAGGTGGTCGATGCACCATCGCAt                                                                                                                                           >  2:3036861/1‑150 (MQ=255)
            gTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGTTTACCCCCCAAAAGGTGGTCGCATCACCAGAGCAtta                                                                                                                                         >  1:3030864/1‑149 (MQ=255)
              tttGCGCCTGTCTTGAGGGCAATAGTAGAGAGGCGTTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCAATTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                       <  2:5111976/150‑1 (MQ=39)
              tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                       <  2:8401084/150‑1 (MQ=255)
               ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                                      <  2:7014522/150‑1 (MQ=37)
               ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                                      <  1:7789583/150‑1 (MQ=37)
               ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgcg                                                                                                                                      <  1:1439329/150‑1 (MQ=255)
                tgtgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                     <  2:3828108/147‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                     <  1:1581254/150‑1 (MQ=25)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                     <  1:2357530/150‑1 (MQ=25)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                     <  2:177000/150‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                     <  1:5136882/150‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                     <  1:4149537/150‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                     <  1:6550190/150‑1 (MQ=255)
                tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                     <  1:5885514/150‑1 (MQ=25)
                   gcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTa                                                                                                                                  >  2:1423786/1‑150 (MQ=255)
                       ggCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGTATTCGCGGGTAAAcc                                                                                                                              <  1:8947558/150‑1 (MQ=255)
                        gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAccc                                                                                                                             <  2:1825920/150‑1 (MQ=255)
                          ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGc                                                                                                                           <  2:8898291/150‑1 (MQ=25)
                             tgggCAATAGTCTTGAGGCGGTATTATTCTTAAAAATACTGTTTAAGATTTTTTATATAACCATTTAAGAGAACACCCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTcc                                                                                                                        <  2:7044180/149‑1 (MQ=18)
                               ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                      >  2:2720996/1‑150 (MQ=255)
                               ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAAAACGCCTCTGGTAGGAGCGGGTTTAACCCCGAATACCGTGGTGGATGGAACATCGCAGTCGCCGGTAAACCCCCTccaa                                                                                                                      >  1:7990189/1‑148 (MQ=255)
                                  aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACACCCGAATACGGTGGTCGATGCACCAACGCATTCGCGGGTAAACCCGCTCCTACCa                                                                                                                   >  1:613616/1‑150 (MQ=255)
                                    tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCCCCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                                 <  2:8214377/150‑1 (MQ=255)
                                    tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                                 <  2:3603024/150‑1 (MQ=255)
                                     agtagtGAGGCGGTATTATTCTTAAAAATACTGTTTAAGAATTTTAATATAACCATTTAAGCGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGTTGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAggg                                                                                                                <  2:7682038/150‑1 (MQ=38)
                                     agtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAAACGCGAATACCGTGGTCGATTCACCATAGCAATCGCGGGTAAAACCCCTCCAACCAggg                                                                                                                >  1:6317822/1‑150 (MQ=255)
                                           gAGGCGGGATTATAAATAAAAATACTCTTAAACAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGCAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGCTAAACCCGCTCCTACCAGGGGCCata                                                                                                          <  2:8020514/150‑1 (MQ=38)
                                              gCGGTATTTTTCTTAAAAATACTGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCTTTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATgg                                                                                                       <  2:6570135/150‑1 (MQ=255)
                                              gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATgg                                                                                                       >  1:7628507/1‑150 (MQ=255)
                                                   atttttCTTAAAATTACGGTTTAAGAATTTTTATATAACCATTTATTTGAACCCGCGTCTTGTTGGAGCGGGTTTACCCGCGAATACGTTGGTCGAGGCACCATCGCATTCGCGTGTAAACCCGCTCCTACCAGGGGCCATATGGGCCtt                                                                                                  <  2:2421254/150‑1 (MQ=25)
                                                    ttattCTTAAAAATACTGGTTAAGAATTTTTATAGAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTGTTCCCCCGCCTAAGGTGGTCGAGGCACCCGAGCCTTCGCGGGGGCACCCGCGCCTACCAGCGGCCCGATggggcggc                                                                                                 >  1:5055548/1‑145 (MQ=255)
                                                    ttataCTAAAAAATAATGTTTAATAATTTTTTTATAACCTTTTAAGATAACACGCCTCTGTTTGGAGCGGTTTTACCCGCTATTACGGTGTTCGATGCACCATCGCATTCGCGGTAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTc                                                                                                 <  2:4629576/150‑1 (MQ=255)
                                                     tattatTAAAAATACTGTTAAAGATTTTTTATATAACCTTTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCa                                                                                                <  2:1128194/150‑1 (MQ=255)
                                                     tattCTTAAAAATACTGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCa                                                                                                <  1:7616149/150‑1 (MQ=255)
                                                        tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACCCGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGc                                                                                             <  2:3405783/150‑1 (MQ=255)
                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAccgtat                                                                                                                                              <  2:6988321/99‑5 (MQ=255)
                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                           >  2:5229843/1‑150 (MQ=255)
                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                           >  2:2701404/1‑150 (MQ=255)
                                                                aTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAGGTGGTCGATGCACCAACGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCCTATGGGCCTTCAGGCGTGTTTCt                                                                                     >  2:8653785/1‑150 (MQ=255)
                                                                 tACTTTTTAACAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTTGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTg                                                                                    <  2:845365/150‑1 (MQ=255)
                                                                    tttttACGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTTGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAACCCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCgg                                                                                 <  2:7660381/148‑1 (MQ=255)
                                                                      tttAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                       <  1:1070061/94‑1 (MQ=14)
                                                                      tttAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCt                                                                               >  2:3572037/1‑150 (MQ=255)
                                                                       ttAAGAATTTTTATATAACCATTTAAGAGAACCCGCCTCGGGTAGGAGCGGGTTTACCCGCGAATCCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTc                                                                              <  2:227951/150‑1 (MQ=255)
                                                                       tgaaGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAGGGTGGTCGATGCCCCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTc                                                                              <  2:4791944/148‑1 (MQ=255)
                                                                        tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAAAAAGGTGGGCGATGCACCATCGCATTAGCGGGTAAACCCGCTCCTACCAGGGGCCATAGGGGCCTTCAGGCGTGTTTCTGCGGCTcg                                                                             >  1:6081098/1‑150 (MQ=39)
                                                                        tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTATGCGGCTcg                                                                             >  1:4484190/1‑150 (MQ=255)
                                                                          ataaTTTTTATATAACAATTTAACACAACACGCCTCTGGTAGGAGCGGTTTTACCCGCCAATACGGTGGTCGATGAACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgcg                                                                           <  2:1589587/148‑1 (MQ=255)
                                                                          aGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgcg                                                                           <  2:514491/150‑1 (MQ=255)
                                                                               ttttATATAACCATTTAAGAGAACGCGCCTCTGGTAGGAGCGGGTTTTACCCCGAAAAACGTGGTCGATGCACCAAAGCCTAAGCCGGTAAAACCCCACCAACCAAGGGCCCAAATGGCCATAAGGCCTGTTTTTGCGGCGCGCgggtat                                                                      >  1:6288230/1‑146 (MQ=255)
                                                                                   atataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                                  <  2:4959498/150‑1 (MQ=255)
                                                                                      taACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCa                                                               <  2:5702608/150‑1 (MQ=255)
                                                                                       aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCCCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                              >  2:200790/1‑150 (MQ=255)
                                                                                               acgagaCCCCGCCTGGGGTAGGCGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAGCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                                      <  2:5831930/148‑1 (MQ=38)
                                                                                                    aaCACGCCTCTGGTAGGAGCGGGTTTACCCGCGCATACGGGGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCAGATGGGCCTTCAGGCGGGTGTCTGCGGCTCGCGGCTCTTGTACCAGGCCACGTACAgcg                                                 >  2:5644890/1‑150 (MQ=38)
                                                                                                    aaCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAGGCACCATCGCCATCGCGGGGAAACCCGCGCCGACCCGGGGCCCGCGGGGCCGTCAGGCGGGTTTCGGCGGCACGCGGATCGTGTGCCCGGCCCCGGACAgcg                                                 >  2:1072483/1‑150 (MQ=255)
                                                                                                       acGCCTCTGGTAGGAGCGGGGTTAACCGCGAATACGGTGGTCGATGCACCATCGCAATCGCGGGTAAACCCCCTCCGACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTTTAACAGGCCACGTACAGCGCCg                                              >  1:2976087/1‑150 (MQ=255)
                                                                                                           ctgtGGTCGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGAAACCCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                          <  1:7356338/147‑1 (MQ=255)
                                                                                                            tctGGTAGGAGCGGGTTTACCCGCGACTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTTCGGCTCGCGGATCTTTTACCAGGCCACGTACAGCGccggccgg                                         >  2:7069500/1‑150 (MQ=255)
                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCCTCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAca                                     <  2:8798773/150‑1 (MQ=255)
                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATCCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAca                                     <  1:1423786/150‑1 (MQ=255)
                                                                                                                 tcggcgCGGGTTTCCCCGCTAATAGGGTGGTCGATGCACCATCGATTTCGCGGTAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                    <  1:7029733/145‑1 (MQ=255)
                                                                                                                 tAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                    <  2:5401090/150‑1 (MQ=255)
                                                                                                                     aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                <  2:7628507/150‑1 (MQ=255)
                                                                                                                     aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                <  2:6401974/150‑1 (MQ=255)
                                                                                                                     aGCGGGTTTACCCGCGAATACGGGGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                <  2:4199532/150‑1 (MQ=255)
                                                                                                                               aCCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg                       >  1:550184/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                                  >  2:8832238/1‑89 (MQ=14)
                                                                                                                                 ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  1:7139244/1‑150 (MQ=25)
                                                                                                                                 ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  2:1902247/1‑150 (MQ=25)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  1:968244/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  1:2112917/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  2:7207346/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  1:6383061/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  2:622822/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                     >  2:6160820/1‑150 (MQ=255)
                                                                                                                                  cgcgAATCCGGTGGGCGCTCCACCAGCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                    <  1:3494680/150‑1 (MQ=255)
                                                                                                                                  cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                    >  2:2069406/1‑150 (MQ=255)
                                                                                                                                  cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                    <  2:6850531/150‑1 (MQ=255)
                                                                                                                                   gcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg                   >  1:331475/1‑150 (MQ=37)
                                                                                                                                   gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg                   >  2:3140336/1‑150 (MQ=255)
                                                                                                                                      aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa                >  1:2647286/1‑150 (MQ=255)
                                                                                                                                        tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACg              >  1:3581040/1‑150 (MQ=255)
                                                                                                                                        tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACg              <  2:294346/150‑1 (MQ=255)
                                                                                                                                         aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga             <  2:7701644/150‑1 (MQ=255)
                                                                                                                                           ggtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCAGATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatg           >  2:3599549/1‑150 (MQ=255)
                                                                                                                                              ggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATc        <  2:5873063/150‑1 (MQ=255)
                                                                                                                                               gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATcc       >  2:7466897/1‑150 (MQ=255)
                                                                                                                                                tCGATGCACCCTCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc      <  2:4502031/150‑1 (MQ=255)
                                                                                                                                                tCGATCAACCATCGATTTCGCGGCTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc      <  1:3138440/150‑1 (MQ=255)
                                                                                                                                                  gATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                    <  2:2477034/134‑1 (MQ=255)
                                                                                                                                                  gATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                    >  1:2477034/1‑134 (MQ=255)
                                                                                                                                                    tGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCa  >  1:8764432/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCA  >  NC_002947/5899612‑5899908

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: