Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,340 C→G 9.9% intergenic (+57/‑86) PP_3338 → / → PP_3339 cytochrome bd‑type quinol oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,778,3400CG 9.9% 301.6 / 3.7 117intergenic (+57/‑86)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/71);  new base G (11/2);  total (44/73)
Fisher's exact test for biased strand distribution p-value = 3.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGC  >  NC_002947/3778193‑3778489
                                                                                                                                                   |                                                                                                                                                      
gATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgcgc                                                                                                                                                      >  1:4201052/1‑150 (MQ=255)
  tCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATTCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGACGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCt                                                                                                                                                    >  2:789781/1‑150 (MQ=255)
   cAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTg                                                                                                                                                   <  1:828962/150‑1 (MQ=255)
   cAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTg                                                                                                                                                   <  2:596870/150‑1 (MQ=255)
    aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                  <  2:1529980/150‑1 (MQ=255)
     gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCg                                                                                                                                                 >  2:1685569/1‑150 (MQ=255)
      cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                                <  1:1528051/150‑1 (MQ=255)
      cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                                <  1:2741053/150‑1 (MQ=255)
      cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGccc                                                                                                                                                >  1:4332736/1‑148 (MQ=255)
           gggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaa                                                                                                                                           <  2:292669/150‑1 (MQ=255)
           gggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaa                                                                                                                                           >  2:2898448/1‑150 (MQ=255)
            ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGACcaaaa                                                                                                                                          >  1:4149462/1‑150 (MQ=255)
            ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaaa                                                                                                                                          <  2:53256/150‑1 (MQ=255)
             gAAGCCGCAGTGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCCAAAAg                                                                                                                                         >  1:2578853/1‑150 (MQ=255)
             gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                         <  1:4177054/150‑1 (MQ=255)
              aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGGGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                        <  1:4403948/150‑1 (MQ=255)
               aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                      <  1:751583/150‑1 (MQ=255)
               aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                      <  1:751610/150‑1 (MQ=255)
               aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                      <  2:2924179/150‑1 (MQ=255)
               aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAGgc                                                                                                                                      >  1:1152961/1‑150 (MQ=255)
                gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                     <  1:3630090/150‑1 (MQ=255)
                gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                     <  2:3769050/150‑1 (MQ=255)
                 ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGAAAAAAATgcgc                                                                                                                                    >  2:917577/1‑150 (MQ=255)
                 ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                    <  2:2659168/150‑1 (MQ=255)
                   gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGACAAAAAGgcgccc                                                                                                                                  >  2:2155185/1‑148 (MQ=255)
                   gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACAAAAAAGCGcac                                                                                                                                  >  1:3219785/1‑150 (MQ=255)
                   gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCACTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                  >  2:4611488/1‑150 (MQ=255)
                      gcgcCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCTCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcaccac                                                                                                                               <  2:2679778/150‑1 (MQ=255)
                        gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACgg                                                                                                                             >  1:3741083/1‑150 (MQ=255)
                        gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAAGCGccccaccc                                                                                                                             >  1:2310132/1‑148 (MQ=255)
                        gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGGCAAAAACGCGccccaccg                                                                                                                             >  1:3740877/1‑148 (MQ=255)
                             cccAGAGCCATGGCTTCTTGGTGGTTGGCACACTGTTACTCGTGCCCTGTTGGCTTGGCTGTGCCGAGTGTTCTTGTTGGCCGCGGGCTTGCGCAGCCTTTGAGCCCGCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACtt                                                                                                                        <  2:3556534/150‑1 (MQ=255)
                                   gCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCCACAAAGCGCCCCCCCGAACATggaaag                                                                                                                  >  1:233591/1‑145 (MQ=255)
                                    ccAAGGCTTCTTGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTTATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGc                                                                                                                 <  2:2078122/150‑1 (MQ=255)
                                         gCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACGTTTAAAGCCCCCa                                                                                                            >  2:1680492/1‑150 (MQ=255)
                                                gCTGGTAGGCACCCGGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGTGGGCCGCCGGCTTGCGCAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATg                                                                                                     <  2:3943457/150‑1 (MQ=255)
                                                  tGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGGCGGCGGCAAAAAATCGCCCCACCCAACATGTAAAGCCCCCAAcccaagcg                                                                                                   >  1:4818526/1‑143 (MQ=255)
                                                  gggTAGGAACAGTGTTCATCCTGCCCTGTGCGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCTTTGAGCCCCCGCCGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCg                                                                                                   <  2:3907240/149‑1 (MQ=255)
                                                      aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACGGCACTTGTACAGCCCGCAACGGATGCGGCtt                                                                                               >  2:4229500/1‑150 (MQ=255)
                                                       ggCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCGCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTa                                                                                              >  2:1450608/1‑150 (MQ=255)
                                                         cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGGGCACCACGGTACTTGTACATCCCGCAACGGATGCGGGTTAAt                                                                                            >  2:4241719/1‑150 (MQ=255)
                                                             cTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTAAAGCCCGCAACGGATGCGGCGTAATgcgc                                                                                        >  2:1133497/1‑150 (MQ=255)
                                                             cTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGACCAAAAGGCGCCCCACCGAACGTGTTCAGCCCCCCACCGCTGCGGCTGGATgcgc                                                                                        >  2:1133458/1‑150 (MQ=255)
                                                                  ctcccTGCCCCGTTCGCTTGGGTTATCCGACTGTTCTTTTTTGCCGCCGGCTTGGGCAGCAATTGAGCCGTCGGCTTTTTGCGCCTGGTACAAATTGCGGCACCACGGAACTTGTACAGCCCGCAACGGATGCTGCTTAATGCGCGCCat                                                                                   <  2:2608356/147‑1 (MQ=255)
                                                                  cATCCGGCCCTGTGCGCTGGGCTGATCCGACGGTTCATGTTGGCCGCGGGCATGCGCCGCCATTGAGCCCTCGCGGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCat                                                                                   <  2:2608155/150‑1 (MQ=255)
                                                                    tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTTCGCCGGCGACAAAAAGGCCCACCACGGGAAGATTTCATCCCCCAACGGATTCGGCTTAATGCGCGCCagat                                                                                 >  2:2398167/1‑150 (MQ=255)
                                                                     ccTGCCCTGTTCGCTTGGCTGATCCGCCTGTTCATGCTGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATAtt                                                                                <  2:2190844/150‑1 (MQ=255)
                                                                      cTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTc                                                                               <  2:2881914/150‑1 (MQ=255)
                                                                       tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGCCAAAAAAGCGCCCCACCGCACATGTACAGCCCCCCACGGCAGGCGCGTAATGCGCGCCCGATTCg                                                                              >  2:852961/1‑150 (MQ=255)
                                                                       tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCCAAAAGGCGCACCCCCGAACGAGTACAGCCCCCAACCGCTGCGGCGTAATGAGCGCCAAAttag                                                                              >  1:799188/1‑148 (MQ=255)
                                                                       tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACCAAAAGGCGCACCACCGCACATGTAAAGCCCCCCACCGCTGCCGCGTAATTCGCGcccaaatcg                                                                              >  2:701130/1‑143 (MQ=255)
                                                                        gCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGACCGTGTACAGCCCGCAACGGATGCGGCGTAATGCGCGCCCTATTCGc                                                                             >  2:2876836/1‑150 (MQ=255)
                                                                          ccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                           <  2:4201052/150‑1 (MQ=255)
                                                                                cGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAAGCGCACCACCGAACATTTAAAGCCCCCCACCGGTGGGGCTTAAATCGCGCCCTAATAGCACGGGTgt                                                                     >  2:1356616/1‑149 (MQ=255)
                                                                                   ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGGCCAAAAGGCCCCCCCACGCAAGTGTACAGCCCCCAACCGCGGCGGCTTAATGCGCGCCAAATTAGCACGGGTGAcac                                                                  >  2:1881554/1‑148 (MQ=255)
                                                                                   ttGGCTGATCCGACTGTTCAGGTGGGCCGCCGGCAGGCGCAGCCTTTGCGCCCCCGCCGTGGTGCGCCGGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                  <  2:97062/150‑1 (MQ=255)
                                                                                       cTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGGACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAAGc                                                              >  2:5090032/1‑150 (MQ=255)
                                                                                           tCCGAGTGTTCATGATGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                          <  1:3813762/150‑1 (MQ=255)
                                                                                           tCCGACTTTTCTTGTTGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGGGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                          <  1:1234166/150‑1 (MQ=255)
                                                                                             ctattGTTCATGTTGGCCGCCGGCATGCGCACCTTTTGAGCCCTCGCGGTGGTGCGCCTGCTAAAAAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                        <  1:4924828/146‑1 (MQ=255)
                                                                                             cGACTGTCCATGTTGGCCGCGGGAATGGGCTGCCATTGAGCCCCCGCGGTGGTGGGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                        <  2:1087181/150‑1 (MQ=255)
                                                                                              gACTGTTCATGATGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacg                                                       <  2:2080295/150‑1 (MQ=255)
                                                                                              gACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCCCGGAACGTGTAAAGCCCGCAACGGATGCGGGTTAATTCGCGCCATATTCGCGCGGGTGACTCATGCTTTcacg                                                       >  2:5046123/1‑150 (MQ=255)
                                                                                               aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                               >  2:1668874/1‑62 (MQ=255)
                                                                                               aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                               <  1:1668957/62‑1 (MQ=255)
                                                                                               aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                               >  2:1668957/1‑62 (MQ=255)
                                                                                                 tGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGTCGCCACAAAGGCGCACCACCGAAAATTTAAAGCCCCCAAAGGATGCGGCTTAAATCGCGCCAAATTAGCGCGGGGGAAACAAGCtttaaactgg                                                    >  2:975531/1‑143 (MQ=255)
                                                                                                   ttCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCAACACCGAACGTGTAAAGCCCCCAACCGATTCGGCGTAAATCGCGCCCAATTAGCACGGGTGACACAAGCGGTcaagagcac                                                  >  2:2334101/1‑150 (MQ=255)
                                                                                                     cctgttgGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGGGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgt                                                <  1:3770104/145‑1 (MQ=255)
                                                                                                      attatgCCCGCCGCAACGCCCACCAATTGACCCCACGCCGTGGGCCCCCAGCCAAAAAAACCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                               <  1:34715/147‑1 (MQ=255)
                                                                                                      atgatgGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCGGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                               <  2:1161107/150‑1 (MQ=255)
                                                                                                        gatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTAAAGCCCGCAACGGATGCGGCTTAAAGCGCGCCCAATTCGCTCGGGTGACTCATGCGTTCAcgtgcacgtgca                                             >  2:4543695/1‑150 (MQ=255)
                                                                                                         tgggCCGCCGGCTTGCGCTGCATTTGAGCCCTCGCGGTGGGGCGCCTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACgtgcacgtgcac                                            <  1:4871345/148‑1 (MQ=255)
                                                                                                           ggCCGCCGGCATGCGCCGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAGGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACAt                                          <  2:1548377/150‑1 (MQ=255)
                                                                                                           ggCCGCCGGCAGGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACAt                                          <  2:1548454/150‑1 (MQ=255)
                                                                                                            gccgccGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCCCCCCGGAACTTGTAAAGCCCCCAACCGCTGCGGGTTAAAGCGCGCCAAATTCGCTCGGGGGACTCAAGCTTTCACCTGCACGGGCCCAta                                         >  2:1964015/1‑149 (MQ=255)
                                                                                                               gccgGCTTGCGCAGCCTTTGAGCCCTGGCGGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                      <  1:4543695/150‑1 (MQ=255)
                                                                                                               gccgGCATGCGCTGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                      <  2:2067152/150‑1 (MQ=255)
                                                                                                               gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                      <  2:4534374/150‑1 (MQ=255)
                                                                                                               gccgGCATGCGCAGCATTTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                      <  2:454541/150‑1 (MQ=255)
                                                                                                                 gggCTTGCGCCGCATTTGCGCCCTCGGGGTGGTGCGCCTGCCAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                    <  1:3228061/149‑1 (MQ=255)
                                                                                                                 gggCTTGCGCAGCCATTGAGCCCTGGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                    <  1:346352/149‑1 (MQ=255)
                                                                                                                 cgGCATGCGGAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGGGCACCACGGAAGTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCAGATCTTTGt                                    <  2:2009687/150‑1 (MQ=255)
                                                                                                                  gtcttGCGCCCCCTTTGAGCCCTCGGCGTGGGGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTc                                   <  1:2763803/146‑1 (MQ=255)
                                                                                                                  ggCTTGCGCAGCAATTGCGCCCTCGCGGTGGTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTc                                   <  1:346047/150‑1 (MQ=255)
                                                                                                                   gCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                  <  2:4175346/150‑1 (MQ=255)
                                                                                                                           aGCCTTTGCGCCCTCGCGGTGGGGGGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcct                          <  1:4043189/150‑1 (MQ=255)
                                                                                                                             cctttGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgc                        <  2:2310132/147‑1 (MQ=255)
                                                                                                                             ccATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgc                        <  2:3307725/150‑1 (MQ=255)
                                                                                                                              cATTGCGCCCCCGCCGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgcc                       <  2:432046/150‑1 (MQ=255)
                                                                                                                               aTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgcct                      <  2:3892907/150‑1 (MQ=255)
                                                                                                                                  gAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGc                   <  1:2898448/150‑1 (MQ=255)
                                                                                                                                       cTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc              <  1:3743290/150‑1 (MQ=255)
                                                                                                                                       cTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCGCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTgn              <  1:1450608/150‑2 (MQ=255)
                                                                                                                                        tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt             <  2:4997986/150‑1 (MQ=255)
                                                                                                                                        gcGCCGGTGTGCGCAGGCCACCAAATGGGCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt             <  1:789781/149‑1 (MQ=255)
                                                                                                                                        ccgcgggggtgCGCCTGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt             <  1:944937/143‑1 (MQ=255)
                                                                                                                                               ggtgCGCCTGCAACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc      <  1:1472509/150‑1 (MQ=255)
                                                                                                                                                 ggcgcCTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACct    <  1:2437965/149‑1 (MQ=255)
                                                                                                                                                   cgcCTGCCACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgc  <  1:4958171/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGC  >  NC_002947/3778193‑3778489

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