Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,902,787 G→C 18.6% 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,902,7870GC18.6% 438.2 / 101.6 194intergenic (‑182/+52)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (68/85);  new base C (12/24);  total (85/109)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-03

Only 100 of 196 total aligned reads displayed.

GAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAA  >  NC_002947/5902638‑5902936
                                                                                                                                                     |                                                                                                                                                     
 agGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGc                                                                                                                                                      <  2:3014167/150‑1 (MQ=255)
   gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                                    <  2:5490110/150‑1 (MQ=255)
   gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCAc                                                                                                                                                    >  1:2437723/1‑150 (MQ=255)
      tGAAGGGGTTTTGCGCCTGGCTTGAGGGCTATAGTAGTGAGGCGGGATTATTCTTAAAATTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAt                                                                                                                                                 <  2:1594211/150‑1 (MQ=255)
       gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                                <  2:2571112/150‑1 (MQ=255)
         aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAAAAAGGTGGTCGATGCACCATCGc                                                                                                                                              >  1:1628247/1‑150 (MQ=255)
           gggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCCTCGCAt                                                                                                                                            >  1:7275917/1‑150 (MQ=255)
             gTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTc                                                                                                                                          <  1:7156023/150‑1 (MQ=255)
                ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgcg                                                                                                                                       <  1:4495857/150‑1 (MQ=255)
                 tGCGCCTGGCTTGCGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                      <  1:3237718/150‑1 (MQ=255)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAACATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                      <  1:6373628/150‑1 (MQ=255)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                      <  1:2897629/150‑1 (MQ=25)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                      <  1:5348723/150‑1 (MQ=25)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                      <  1:5145820/150‑1 (MQ=25)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                      <  1:5664134/150‑1 (MQ=25)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                      <  1:6844708/150‑1 (MQ=255)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                      <  2:5622860/150‑1 (MQ=255)
                 tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                                      <  1:7426910/150‑1 (MQ=37)
                 tGCGCCTGGCTTGAGGGAAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                      <  1:2471874/150‑1 (MQ=255)
                  gcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                                     <  2:423356/150‑1 (MQ=255)
                      cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTaaa                                                                                                                                 >  2:471562/1‑150 (MQ=255)
                       tGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                                <  1:976594/150‑1 (MQ=255)
                       tGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                                <  2:5969732/150‑1 (MQ=255)
                         gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAccc                                                                                                                              <  1:6365762/150‑1 (MQ=255)
                                  cAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCCCGAATAAGGTGGTCGAATCACCAACGCAATAGCCGGTAAAACCCCTccaaac                                                                                                                     >  1:4718022/1‑146 (MQ=255)
                                   aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAACCCAATCGCGGGTAAACCCGCTCCTACCa                                                                                                                    >  2:1470135/1‑150 (MQ=255)
                                   aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCCAAAAACGTGGGCGATGCACCATAGCAATCGCGGGTAAAACCCCTCCGACCa                                                                                                                    >  2:2455220/1‑150 (MQ=255)
                                     tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                                  <  2:4842263/150‑1 (MQ=255)
                                     gagtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATACCCATTTAAGAGAACCCGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                                  <  2:4842836/149‑1 (MQ=255)
                                        tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGc                                                                                                               >  2:1746875/1‑150 (MQ=255)
                                          gtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGccc                                                                                                             >  2:6352433/1‑149 (MQ=255)
                                               gCGTGATTATTATTAAAAATAATGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTTGGAGGGGGTTTACCCGCTAATACGGTGGTCGTTGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATgg                                                                                                        <  1:5462321/150‑1 (MQ=38)
                                               gCGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTACGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAGGGTGGTCGATGCCCCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATgg                                                                                                        <  2:4720969/150‑1 (MQ=38)
                                               gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                                >  1:2651652/1‑126 (MQ=255)
                                                 ggtattattCTTAAAAATACTTTTTAAGAATTTTTATATAACCATTTAATAGAACACGCCTCTGGTTGGAGGGGTTTTACCCGCTAATAGGGTGGTCGATGCCCCATCGAATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGc                                                                                                      <  2:2678195/147‑1 (MQ=255)
                                                        ttCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATAACGTGGTCGATTCACCAACCCAGTAGCCGGGAAAACCCCGCCAACCCCGGGCCCGATGGGCCGGCAgg                                                                                               >  1:3804856/1‑150 (MQ=255)
                                                        ttCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAACGTGGTCGATGCACCATCGCAATAGCGGGTAAACCCGCTACTAACAGGGGCCAAATGGGCCTTCAgg                                                                                               >  1:2965699/1‑150 (MQ=39)
                                                         tCTTAAAAATACTGTTTCAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGc                                                                                              <  2:1944431/150‑1 (MQ=255)
                                                         tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGc                                                                                              >  2:2029102/1‑150 (MQ=255)
                                                           ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGATGCACCAACGCAGTCGCGGGTAAAACCGCTCCAACCAGGGGCCCTATGGGCCGTCAGGCgt                                                                                            >  1:4158811/1‑150 (MQ=255)
                                                           ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGCATACGGTGGTCGATGCACCATCGCAGTCGCGGGTAAACCCCCTCCGACCCGGGGCCATATGGGCCGTCAGGCgt                                                                                            >  2:3132362/1‑150 (MQ=39)
                                                           ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                            >  2:4164623/1‑150 (MQ=255)
                                                           ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATAGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCGTCAGGCgt                                                                                            >  2:2089450/1‑150 (MQ=255)
                                                             aaaaaTACTGTTTAAGATTTTTTATATAACCATTTAAGAGAACACGCCTCTGGTTGGAGCGGGTTTACCCGCGAATACGGGGGTCGATGCACCATCGCATTGGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgtgt                                                                                          <  2:2007111/150‑1 (MQ=255)
                                                             aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAAAACGGTGGTCGATTCACCAACCCAATAGCGGGTAAACCCCCTCCTACCAGGGGCCATATGGGCCTTAAGGCgggt                                                                                          >  2:3275768/1‑150 (MQ=38)
                                                                  tACTGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGTTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTg                                                                                     <  2:2045899/150‑1 (MQ=255)
                                                                   aCTGTTTAATAATTTTTATATAACCTTTTAAGCCAACCCGCCTCTGGTAGGAGCGGGTTTACCCGCTAATAGGGTGGTCGATGCACCATCGCATTCGGGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGc                                                                                    <  2:1876561/150‑1 (MQ=38)
                                                                   aCTGTTTAATAATTTTTATATAACAATTTAAGACAACCCGCCTCGGGTTGGAGGGGGTTTACCCGCGAATACGGGGGTCGCGGCACCATCGCATTCGCGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGc                                                                                    <  1:471562/150‑1 (MQ=255)
                                                                       tttAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGGTTCTGCGGCt                                                                                >  2:3471123/1‑150 (MQ=255)
                                                                          aaGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                             >  2:6028946/1‑150 (MQ=255)
                                                                          aaGAATTTTTATAAACCCATTTAAGACAACCCGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                             <  2:852880/150‑1 (MQ=38)
                                                                                ttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAGGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCt                                                                       <  2:7336955/150‑1 (MQ=255)
                                                                                  ttATATAACCATTTAAGAGAACACGCCTCTGGTTGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                     <  2:3925841/150‑1 (MQ=255)
                                                                                     tataACCATTTAAGAGAACACGCCTGTGGTCGGAGCGGTTTTACCCGCGAATCGGGTGGTCGATGCACCAGCGCATGCGCGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                                  <  1:1711364/150‑1 (MQ=255)
                                                                                     tataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGCTGCACCATAGCATTGGCGGGTAAACCCGCTCCAACCACGGGGCCTATTGGCCTTCCGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                                  >  1:1159999/1‑150 (MQ=255)
                                                                                     tataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCCCGAATAAGGTGGTCGATTCACCATCGCATTAGCGGGTAAACCCGCTCCGAACAGGGGCCCGAATGGCCATCCTGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                                  >  1:2310966/1‑150 (MQ=255)
                                                                                        aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCAGATGGGCCTGCAGGCGTGTTTCTGCGGCTCGCGGCTCGTGTACCAg                                                               >  2:5277768/1‑150 (MQ=255)
                                                                                        aaCCATTGAAGACAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAGGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                               <  2:3768702/150‑1 (MQ=255)
                                                                                                   agaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAg                                                    <  2:6804490/150‑1 (MQ=255)
                                                                                                        acGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCCTCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACCGCGCCg                                               >  2:5434378/1‑150 (MQ=255)
                                                                                                         cGCCTCTTGTTGGCGGGGTTTTACCCGCTAATAGGGTGGTGGATGCACAATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                              <  1:890037/150‑1 (MQ=255)
                                                                                                            ctctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGTGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                           <  2:6706239/150‑1 (MQ=255)
                                                                                                                 gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAca                                      <  2:1064673/150‑1 (MQ=255)
                                                                                                                  tAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                     <  2:2435608/150‑1 (MQ=255)
                                                                                                                      aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                 <  1:6097106/150‑1 (MQ=255)
                                                                                                                      aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGCTCTTGTACCAGGCCACGTACCGCGCCGGCAGGAACagcagc                                 >  2:2636230/1‑150 (MQ=255)
                                                                                                                       gCGGGTTTACCCGCGAAAACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCCGAGGGGCCTTCAGGCGTGGTTCTGCGGCTCGCGGGTCTTGTACCCGGCCACGTCCAGCGCCGGCAGGCACAGCAGCg                                >  1:3899778/1‑150 (MQ=38)
                                                                                                                          ggTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCa                             <  2:7275917/150‑1 (MQ=255)
                                                                                                                               aCCCGCGAATTGGGGGGTCGCTGCACCTTCGCTTTGGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg                        <  1:1341114/150‑1 (MQ=255)
                                                                                                                               aCCCGCGAATAGGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg                        <  2:1826704/150‑1 (MQ=255)
                                                                                                                                 ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                      >  2:4776911/1‑150 (MQ=25)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                      >  1:1819314/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                      <  1:609056/150‑1 (MQ=255)
                                                                                                                                 ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                      >  1:1510834/1‑150 (MQ=255)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                     >  1:4148427/1‑102 (MQ=17)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                     <  2:4148302/102‑1 (MQ=17)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                     <  2:4148434/102‑1 (MQ=17)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                     >  2:75136/1‑150 (MQ=37)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                     >  1:2137069/1‑150 (MQ=37)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                     >  2:3076854/1‑150 (MQ=25)
                                                                                                                                  cgcgAATGCGATGGTGCATCGACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                     >  1:6776639/1‑150 (MQ=37)
                                                                                                                                  cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                     <  2:1221238/150‑1 (MQ=255)
                                                                                                                                   gcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg                    >  2:2517929/1‑150 (MQ=37)
                                                                                                                                   gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg                    >  1:1094163/1‑150 (MQ=255)
                                                                                                                                      aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa                 >  1:4577552/1‑150 (MQ=255)
                                                                                                                                      aaTACGGTGGTCCATCCACAATCCAATTCGCGGTAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa                 <  1:2208377/150‑1 (MQ=255)
                                                                                                                                         tcGGTGGTCTCTGCACCATCGCATTCGCGTGTAACCCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga              <  1:3377449/149‑1 (MQ=255)
                                                                                                                                         aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga              <  2:3197298/150‑1 (MQ=255)
                                                                                                                                         aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga              <  2:1783753/150‑1 (MQ=255)
                                                                                                                                          cGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat             >  2:1673293/1‑150 (MQ=255)
                                                                                                                                           ggtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatg            >  1:4062881/1‑150 (MQ=255)
                                                                                                                                             tggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatgat          >  1:6649020/1‑150 (MQ=255)
                                                                                                                                               gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCACCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATcc        >  1:1010944/1‑150 (MQ=255)
                                                                                                                                                tCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc       >  2:568423/1‑150 (MQ=255)
                                                                                                                                                tCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc       >  2:5749731/1‑150 (MQ=255)
                                                                                                                                                tCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc       >  2:6401791/1‑150 (MQ=255)
                                                                                                                                                  gATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGc     >  2:2879687/1‑150 (MQ=255)
                                                                                                                                                  gATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGc     >  2:2100607/1‑150 (MQ=255)
                                                                                                                                                   aTGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGcc    <  2:1425691/150‑1 (MQ=255)
                                                                                                                                                   aTGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCAGATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGACCAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGcc    >  2:6279505/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAA  >  NC_002947/5902638‑5902936

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: