Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,703,797 G→A 8.6% intergenic (+82/‑85) hslU → / → PP_5002 protease HslVU ATPase subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,703,7970GA 8.6% 546.7 / 30.1 185intergenic (+82/‑85)hslU/PP_5002protease HslVU ATPase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (78/91);  new base A (4/12);  total (82/103)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 195 total aligned reads displayed.

CGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCACCTAC  >  NC_002947/5703648‑5703944
                                                                                                                                                      |                                                                                                                                                    
 ggcgTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGGCACCCCGCGCCAGCCCGCGCCACCCCACCCCCGTGGGCCGTGGGGGTACGTGCggggg                                                                                                                                                     >  2:5478084/1‑150 (MQ=255)
  gcgTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGGCAACCCCCACCAGCACCCCCCCCCACAACCCCAGGGCCCGTGGTGGTACCTGCGGgggc                                                                                                                                                    >  1:1676519/1‑150 (MQ=255)
    gTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgg                                                                                                                                                  >  2:3801914‑M2/1‑146 (MQ=255)
     tATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcggg                                                                                                                                                 >  1:1239237‑M2/1‑145 (MQ=255)
     tATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcggg                                                                                                                                                 >  2:5402134‑M2/1‑145 (MQ=255)
      aTGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggc                                                                                                                                                >  1:6345917‑M2/1‑144 (MQ=255)
      aTGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggc                                                                                                                                                >  1:159030‑M2/1‑144 (MQ=255)
       tgtgAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggct                                                                                                                                               >  1:4515021‑M2/1‑143 (MQ=255)
         tgAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttg                                                                                                                                             >  1:1044389‑M2/1‑141 (MQ=255)
         tgAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttg                                                                                                                                             >  1:1191925‑M2/1‑141 (MQ=255)
             cAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccg                                                                                                                                         >  1:1998814‑M2/1‑137 (MQ=255)
              aGCCACCTGGGTGAGCTGGCACAGAACTAAGACCTGTCGCGTTACATTGTTTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCGGCACGCGCCACCACAACCCTTGGGCCGGTGGTGATCCGTGCGGGGgcggggttgcccgc                                                                                                                                        >  1:1788073‑M2/1‑136 (MQ=255)
                  aCCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGGAAGCCCCCTCCAGCACCCGCCACCCCAACCCGTTGGGCGGTGGTGATACGGGCGGGGgcgggggtgcccgccaag                                                                                                                                    <  2:472196‑M2/150‑19 (MQ=255)
                  aCCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCCAACCCCCAACACCACACCCCCCCACAAACCCAGGGCCGGGGGTTTTACGTGCGGGGGGGGGCGTTGCCCCCGag                                                                                                                                    >  2:3074619/1‑149 (MQ=255)
                             cTGGCACAGACCGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGaa                                                                                                                        <  2:6691750/150‑1 (MQ=255)
                              tGGCACACCACGACCCCCTGCCGCGTTACATGCTGTAAGCCCGGTTGCGACAGTTTCGCGGGCAAGCCCGCTCCTGCCCGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAt                                                                                                                       <  2:4852921/150‑1 (MQ=255)
                               ggCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCCAGCCCGCGCCAGCCCGCGCCACCACAAACCGTGGGCGGGGGGTGGTCCGTGCGGGGGGGGGCGTGGCCCCGGATAGGCCCCTGAAtc                                                                                                                      >  2:2257635/1‑150 (MQ=255)
                                  acaGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCCAGCCCCCGCCGGCCCCCGCCCCCCCCACCCCGGGGGCGGGGGGGTTCCCGGCGGGGGCGGGCGTGGCCCCGCAGCGGCCCGTGGAtacca                                                                                                                   >  2:3321546/1‑149 (MQ=255)
                                    aGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGGCAGGCCGCGACTGCACCCGCCCCCACAACCCGAGTGGCGGGGGTGTTCAGTGGGGGGGGGGGGGTGGCCCCGCAGAGGCCCGTGAAtaacacc                                                                                                                 >  2:2032298/1‑144 (MQ=255)
                                     gAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccctccc                                                                                                                <  2:2170066/150‑1 (MQ=255)
                                      aaCGAAGACCTGGCGCTTTAATTTCTTTACGACCGGTTTCGACAGTTTCGGGGGCACGCCCGCTCCGGCACGCGCCACCACACCCCTGGGGCCTGTGGTGATCGGTGGGGGGGGGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTcccc                                                                                                               <  1:2128760/150‑1 (MQ=255)
                                       aCGAAGACCTGTCGCGTTACTTTCTGTATGACCGGTTTCGACATTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTccccc                                                                                                              <  2:1778518/150‑1 (MQ=255)
                                         gAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTTTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCAt                                                                                                            <  2:281498/150‑1 (MQ=255)
                                                  tCGCTTTACATTCTTTAAGACGGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGaa                                                                                                   <  2:1162982/150‑1 (MQ=255)
                                                     cgTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCCCCCCCCCGCCCCCACCCCCCCCCCCCACCCCAGGGCCCGGGGGGGTCCCGGCCCGGGCGGGCGGGCCCCCGCACACGCCCCAGGATCCCACCCCCCGAAGCAggccgcg                                                                                                >  2:3257565/1‑145 (MQ=255)
                                                       ttACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTGGGGCCTGTGGTGATCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCtgt                                                                                              <  2:6470325/150‑1 (MQ=255)
                                                         aCATTCTGTATGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCGGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGc                                                                                            <  2:1884121/150‑1 (MQ=255)
                                                           aTTCTTTAAGACCGGTTTCGACAGTTTCGCGGGCACGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCtt                                                                                          <  1:6779650/150‑1 (MQ=255)
                                                             tCTTTAAGACCGTTTTCGACAGTTTCGCGGGCACGCCCGCGCCTGCCCGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGGGGGGGGGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtca                                                                                        <  1:6052728/150‑1 (MQ=255)
                                                                ttttgACGGGTTTCGACTGTTTGGGGGGCCCGCCGGCTCCTGCCCGCGCCACCCCCACGGTTGGGCCTGTGGTGATCGGTGGGGGGGGGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                     <  1:3657746/146‑1 (MQ=255)
                                                                gTAAGACCGGTTTCGCCAGTTGCGGGGGCACGCCGGCTGCTGCACGCGCCCCCACACCCCTTGGGCCTGTGGTGCGCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                     <  2:251/150‑1 (MQ=255)
                                                                gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCACCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGAGCCGTGCGGGGGCGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                     <  1:4449759/150‑1 (MQ=255)
                                                                      ccGGTTTCGACTGTTTCGCGGGCACGCCCGCGCCGGCCCGCCCCCCCCCAACCCTGGGGCCTGTGGTGACCCGGGGGGGGGGGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTccc                                                                               <  1:4950147/150‑1 (MQ=255)
                                                                        ggTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                             <  2:3899641/150‑1 (MQ=255)
                                                                        ggTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                             <  2:4796333/150‑1 (MQ=255)
                                                                        ggTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                             <  2:1593387/150‑1 (MQ=255)
                                                                          tttCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCgaga                                                                           <  2:2427800/150‑1 (MQ=255)
                                                                          tttCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCgaga                                                                           <  2:430428/150‑1 (MQ=255)
                                                                              tactgTTGCGGGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGGGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGAtt                                                                       <  1:5355191/146‑1 (MQ=255)
                                                                              gACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGAtt                                                                       <  2:1944982/150‑1 (MQ=255)
                                                                                  gTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGc                                                                   <  2:1605895/150‑1 (MQ=255)
                                                                                            cAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCg                                                         >  1:6842493/1‑150 (MQ=255)
                                                                                             acgCCCGCTCCGGCACGCGCCACCACACCCCTGGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgc                                                        <  1:5523650/148‑1 (MQ=255)
                                                                                             aaGCCCGCTCCTGCTCGCGCCACCACACCCCTTGGGCCTGTGTTGATCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgc                                                        <  1:3827159/150‑1 (MQ=255)
                                                                                             aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcac                                                         <  1:2824928‑M2/150‑94 (MQ=255)
                                                                                             aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcac                                                         <  1:3043605‑M2/150‑94 (MQ=255)
                                                                                             aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgc                                                        >  2:6613899/1‑150 (MQ=255)
                                                                                             aaGCCCGCCCCCGCACGGATCACCACAGGCCCAAGGGTTGTGGTGGCGCGTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgc                                                        >  1:5921976/1‑150 (MQ=17)
                                                                                                cccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctg                                                     >  1:1441399/1‑150 (MQ=255)
                                                                                                 ccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGTGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctgc                                                    >  1:4839168/1‑150 (MQ=255)
                                                                                                 ccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaagctgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatac                                                     <  1:4343751‑M2/150‑98 (MQ=255)
                                                                                                 ccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctgc                                                    >  1:2019413/1‑150 (MQ=255)
                                                                                                  cGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacg                                                    <  1:4127351‑M2/150‑99 (MQ=255)
                                                                                                   gCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgc                                                   <  2:4626061‑M2/150‑100 (MQ=255)
                                                                                                   gCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctgcct                                                  >  1:6519086/1‑150 (MQ=255)
                                                                                                    cTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgcc                                                  <  2:3320622‑M2/150‑101 (MQ=255)
                                                                                                     tCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaagctgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccg                                                 <  1:2417318‑M2/150‑102 (MQ=255)
                                                                                                     tCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccg                                                 <  1:3911433‑M2/150‑102 (MQ=255)
                                                                                                     tCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcat                                                                                       <  2:1684608/111‑1 (MQ=255)
                                                                                                     tCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccg                                                 <  1:2101494‑M2/150‑102 (MQ=255)
                                                                                                          cACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCa                                           >  1:46610/1‑150 (MQ=255)
                                                                                                           aCGCGCCACCACAACCCTTGGGCCTGTGGTGAGCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCAt                                          <  2:1397311/150‑1 (MQ=255)
                                                                                                             gcgcCACCACACCCCTTGCGCCTGTGGTGAGCCGTGCGGGGGGGGGCTCGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCa                                        <  1:2032298/150‑1 (MQ=255)
                                                                                                             gcgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCa                                        >  1:1370431/1‑150 (MQ=255)
                                                                                                             gcgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCa                                        >  2:5898953/1‑150 (MQ=255)
                                                                                                             gcgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCa                                        >  1:2520775/1‑150 (MQ=255)
                                                                                                             gcgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCa                                        >  1:3096051/1‑150 (MQ=255)
                                                                                                              cgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCaa                                       <  2:4402538/150‑1 (MQ=255)
                                                                                                              cgcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgcca                                                                               >  1:2202704‑M2/1‑40 (MQ=255)
                                                                                                               gcCACCACAACCCTTGGGCCTGTGGTGGCGCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctg                                       <  2:3616826‑M2/150‑112 (MQ=255)
                                                                                                               gcCACCACAACCCTTGGGCCTGTGGTGGCGCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctg                                       <  1:5647779‑M2/150‑112 (MQ=255)
                                                                                                               gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAAc                                      >  2:2883513/1‑150 (MQ=255)
                                                                                                               gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAAc                                      <  1:5898953/150‑1 (MQ=255)
                                                                                                               gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgaccggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctg                                       <  2:6797766‑M2/150‑112 (MQ=255)
                                                                                                                 caccacAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCt                                    >  2:255417/1‑150 (MQ=255)
                                                                                                                  accacAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTg                                   >  2:3499151/1‑150 (MQ=255)
                                                                                                                  accacAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTg                                   >  2:4584895/1‑150 (MQ=255)
                                                                                                                   ccacAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGc                                  >  2:698342/1‑150 (MQ=255)
                                                                                                                    cccaACCCTTGGGCCTTTGGTGATCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGca                                 <  2:775494/148‑1 (MQ=255)
                                                                                                                    cacaCCCCTTGGGCGTGTGGTGATCGGTGGGGGGGGGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGca                                 <  1:207275/150‑1 (MQ=255)
                                                                                                                    cacaACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGca                                 >  2:1035060/1‑150 (MQ=255)
                                                                                                                     acaACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGcac                                >  1:6201430/1‑150 (MQ=255)
                                                                                                                       aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaa                              >  1:491478/1‑150 (MQ=255)
                                                                                                                       aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaa                              >  2:5079942/1‑150 (MQ=255)
                                                                                                                       aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaa                              >  2:6091508/1‑150 (MQ=255)
                                                                                                                         ccccTTGGGCCTGTGGTGATCGGGGGGGGGGGGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaaa                             <  1:5478084/149‑1 (MQ=255)
                                                                                                                          cccTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctggatcggcgct                             <  1:5167016‑M2/150‑122 (MQ=255)
                                                                                                                          cccTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAg                            >  2:5570323/1‑150 (MQ=255)
                                                                                                                           ccTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGc                           >  2:401150/1‑150 (MQ=255)
                                                                                                                             ttGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctggatcggcgcttcg                          <  1:1060270‑M2/150‑125 (MQ=255)
                                                                                                                             gggggCGTGTGGGGATCGGTGGGGGGGGGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCt                         <  1:6399294/148‑1 (MQ=255)
                                                                                                                               gggCGTGTGGTGATCGGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTcc                       <  1:6667237/150‑1 (MQ=255)
                                                                                                                               gggCCTGTGGTGATCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTcc                       <  2:742822/150‑1 (MQ=255)
                                                                                                                                 gCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCaa                     <  2:6186196/150‑1 (MQ=255)
                                                                                                                                     gtggtgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTATGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGAcc                 <  2:5220658/150‑1 (MQ=255)
                                                                                                                                      tggtgATCCGGGGGGGGGCGGGGTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGAccc                <  1:6569129/150‑1 (MQ=255)
                                                                                                                                        gtgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTc              >  1:3737510/1‑150 (MQ=255)
                                                                                                                                           gcgcGTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGc           >  2:6413280/4‑150 (MQ=255)
                                                                                                                                               gTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCa       >  2:4902536/1‑150 (MQ=255)
                                                                                                                                                    ggggCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGc                                                            <  2:2452806/92‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCACCTAC  >  NC_002947/5703648‑5703944

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: