Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,338 T→G 10.5% intergenic (+55/‑88) 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,3380TG10.5% 265.4 / 4.2 116intergenic (+55/‑88)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/66);  new base G (10/5);  total (45/71)
Fisher's exact test for biased strand distribution p-value = 2.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CGTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAC  >  NC_002947/3778189‑3778485
                                                                                                                                                     |                                                                                                                                                   
cgTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCgtggt                                                                                                                                                     <  2:4259419/150‑1 (MQ=255)
 gTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCTCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCgtggtg                                                                                                                                                    <  1:4365845/150‑1 (MQ=255)
    gATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgcgc                                                                                                                                                 >  1:4201052/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)
        aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                             <  2:209648/150‑1 (MQ=255)
         gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCg                                                                                                                                            >  2:1685569/1‑150 (MQ=255)
          cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                           <  1:1528051/150‑1 (MQ=255)
          cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGccc                                                                                                                                           >  1:4332736/1‑148 (MQ=255)
              tGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGGCaaa                                                                                                                                       >  1:2641933/1‑150 (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)
                   aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                  <  2:2566684/150‑1 (MQ=255)
                   aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                  <  1:751610/150‑1 (MQ=255)
                   aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                  <  1:751583/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)
                           cgccgcCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGGGCACCACg                                                                                                                          >  2:5009600/1‑150 (MQ=255)
                            gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAAGCGccccaccc                                                                                                                         >  1:2310132/1‑148 (MQ=255)
                            gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGGCAAAAACGCGccccaccg                                                                                                                         >  1:3740877/1‑148 (MQ=255)
                            gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGCCAAAAAAGCGccccaccg                                                                                                                         >  1:3741035/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)
                                                       ggTAGGAACACGGTTCATCCGGCCCGGTTCGCTTGGCTGATCCGACTGTCCATGATGGCCGCCGGCTGGCGCAGCCTTTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                              <  2:726352/150‑1 (MQ=255)
                                                          aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACATGTAAAGCCCGCCACGGGTGCGGCtt                                                                                           >  2:1649388/1‑150 (MQ=255)
                                                          aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACGGCACTTGTACAGCCCGCAACGGATGCGGCtt                                                                                           >  2:4229500/1‑150 (MQ=255)
                                                           ggCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCGCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTa                                                                                          >  2:1450608/1‑150 (MQ=255)
                                                             cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGGGCACCACGGTACTTGTACATCCCGCAACGGATGCGGGTTAAt                                                                                        >  2:4241719/1‑150 (MQ=255)
                                                             cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                        <  2:3320556/150‑1 (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)
                                                                                       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)
                                                                                               tCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACTTGTACCGCCCGCAACCGATGCGGCTTAAAGCGCGCCAAATTAGCTCGGGTGAAGCATGCTTTc                                                      >  2:4340255/1‑150 (MQ=255)
                                                                                                 ctattGTTCATGTTGGCCGCCGGCATGCGCACCTTTTGAGCCCTCGCGGTGGTGCGCCTGCTAAAAAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                    <  1:4924828/146‑1 (MQ=255)
                                                                                                 cGACTGTCCATGTTGGCCGCGGGAATGGGCTGCCATTGAGCCCCCGCGGTGGTGGGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                    <  2:1087181/150‑1 (MQ=255)
                                                                                                  gACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACATGTACAGCCCCCAACCGAAGCGGCTTAAAGCGCGCCATATTCCCACGGGTGACTCAAGcataaacg                                                   >  2:1236940/1‑143 (MQ=255)
                                                                                                   aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                          >  2:1668957/1‑62 (MQ=255)
                                                                                                   aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                          <  1:1668957/62‑1 (MQ=255)
                                                                                                   aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAc                                                                                                                                          >  2:1668874/1‑62 (MQ=255)
                                                                                                     tGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGTCGCCACAAAGGCGCACCACCGAAAATTTAAAGCCCCCAAAGGATGCGGCTTAAATCGCGCCAAATTAGCGCGGGGGAAACAAGCtttaaactgg                                                >  2:975531/1‑143 (MQ=255)
                                                                                                         ctggtgggcggcCGGCTGGCGCCGCCTTTGCGCCCGCGCGGTGGGGCGCCTGCGACAAAATGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgt                                            <  1:4404159/144‑1 (MQ=255)
                                                                                                         cctgttgGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGGGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgt                                            <  1:3770104/145‑1 (MQ=255)
                                                                                                          attatgCCCGCCGCAACGCCCACCAATTGACCCCACGCCGTGGGCCCCCAGCCAAAAAAACCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                           <  1:34715/147‑1 (MQ=255)
                                                                                                            gatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTAAAGCCCGCAACGGATGCGGCTTAAAGCGCGCCCAATTCGCTCGGGTGACTCATGCGTTCAcgtgcacgtgca                                         >  2:4543695/1‑150 (MQ=255)
                                                                                                             tgggCCGCCGGCTTGCGCTGCATTTGAGCCCTCGCGGTGGGGCGCCTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACgtgcacgtgcac                                        <  1:4871345/148‑1 (MQ=255)
                                                                                                               ggCCGCCGGCATGCGCCGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAGGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACAt                                      <  2:1548377/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)
                                                                                                                   gccgGCATGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                  <  1:3336210/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)
                                                                                                                       gCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                              <  2:4029617/150‑1 (MQ=255)
                                                                                                                       gCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                              <  2:4175346/150‑1 (MQ=255)
                                                                                                                               aGCCTTTGCGCCCTCGCGGTGGGGGGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcct                      <  1:4043189/150‑1 (MQ=255)
                                                                                                                               aGCCATTGAGCCCTCGCCGTGGTGGGCCAGCGACAAAAAGGCGCACCACGGCACGTGTACAGCCCGCAACGGATTCGGCTTAATGCGCGCCATATTCGCTCGGGGGACTCAGTCGTTCACGTGCACGTGCACAACTTTTTAGGTTAgcct                      >  2:2877536/1‑150 (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)
                                                                                                                                            tCGCCGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt         <  2:188391/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)
                                                                                                                                                     |                                                                                                                                                   
CGTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAC  >  NC_002947/3778189‑3778485

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: