Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,338 T→G 16.9% 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,3380TG16.9% 256.7 / 14.3 125intergenic (+55/‑88)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/65);  new base G (19/5);  total (55/70)
Fisher's exact test for biased strand distribution p-value = 1.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCT  >  NC_002947/3778191‑3778487
                                                                                                                                                   |                                                                                                                                                     
tCGATCATCATCTGGGAATCCGGATCGCCGCCCCAGAGCCAAGGGTTCTTGGTGTTAGTCTCTCTGTTCATCCTGCCCAGTTCGCTTGGCTGATGCGACTGTTGATGATGGACGCCGGCATGCGGTTCCATTGAGCCCTATACGTGGTgc                                                                                                                                                     <  1:392704/150‑1 (MQ=255)
tCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgc                                                                                                                                                     <  2:522455/150‑1 (MQ=255)
tCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgc                                                                                                                                                     <  1:2328606/150‑1 (MQ=255)
   aTCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGcc                                                                                                                                                  >  1:3055574/1‑150 (MQ=255)
   aTCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGcc                                                                                                                                                  >  1:702851/1‑150 (MQ=255)
    tCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGccg                                                                                                                                                 >  1:4132661/1‑149 (MQ=255)
     cAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTg                                                                                                                                                <  1:1693769/150‑1 (MQ=255)
       gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCg                                                                                                                                              <  1:3364341/150‑1 (MQ=255)
          tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCggcc                                                                                                                                           >  2:784770/1‑147 (MQ=255)
          tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCa                                                                                                                                           >  1:323206/1‑150 (MQ=255)
          tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCa                                                                                                                                           >  2:4073752/1‑150 (MQ=255)
          tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACa                                                                                                                                           >  1:3694609/1‑150 (MQ=255)
           cTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaa                                                                                                                                          <  1:1442938/150‑1 (MQ=255)
               gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                      <  1:194330/150‑1 (MQ=255)
               gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACAAAAAg                                                                                                                                      >  2:2796138/1‑150 (MQ=255)
                aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                     <  2:1027273/150‑1 (MQ=255)
                 aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                    <  1:3971468/150‑1 (MQ=255)
                  gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                   <  2:3644166/150‑1 (MQ=255)
                  gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                   <  1:4121997/150‑1 (MQ=255)
                    cGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                 <  1:2338291/150‑1 (MQ=255)
                     gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                <  2:647164/150‑1 (MQ=255)
                      cAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcacc                                                                                                                               <  1:336665/150‑1 (MQ=255)
                        gcgcCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcaccac                                                                                                                             <  2:1979496/150‑1 (MQ=255)
                          gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACgg                                                                                                                           <  2:4137192/150‑1 (MQ=255)
                          gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACgg                                                                                                                           <  2:1985228/150‑1 (MQ=255)
                          gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCCAAAAGGCGCCCCACgg                                                                                                                           >  2:3217408/1‑150 (MQ=255)
                            cgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGaa                                                                                                                         >  2:686009/1‑150 (MQ=255)
                               cccAGAGCCAAGGCTTCATGCTGGTTGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGAGTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACtt                                                                                                                      <  2:1462517/150‑1 (MQ=255)
                                ccAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACCAAAATGCGCCCCACCGAACtag                                                                                                                     >  2:3039322/1‑148 (MQ=255)
                                      ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACAAAAAAGCGCCCCACCGACCATGTaaacc                                                                                                               >  1:878062/1‑146 (MQ=255)
                                      ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGACAAAAAAGCGCACCACCGAACATGTaaaac                                                                                                               >  1:4400073/1‑146 (MQ=255)
                                       catggCTTCTGGCTGTTTGGCACACTGTACACCCGGCCCAGTTCGGTTGGCTGACCCGAGTGTTCATGTTGGCCGCCGCCATGCGCAGCAATTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAACGGCGCACCACGGAACTTGTACAGcc                                                                                                              <  2:2003742/147‑1 (MQ=255)
                                       catggCTTCATGGTGGTAGGCACACTTTTCATCCTGCCCAGTTCGCTTGGCTTTTCCGACTGTTCATGTTGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGcc                                                                                                              <  2:4272706/147‑1 (MQ=255)
                                       cAAGGCTTCATGCTGGTTGGAAAACTTTTCATCCTGCCCTTTTCGCTTGGCTTTTCCGAGTGTTCATGTTGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGcc                                                                                                              <  2:3230899/150‑1 (MQ=255)
                                       cAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTAAAGcc                                                                                                              >  2:3681425/1‑150 (MQ=255)
                                       cAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGACAGTGTACAGcc                                                                                                              >  2:1707643/1‑150 (MQ=255)
                                       cAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACGGAACATGTACAGcc                                                                                                              >  2:3681288/1‑150 (MQ=255)
                                       cAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACCAAAAGGCGCCCCCCCGAACATGTACAGcc                                                                                                              >  2:1707661/1‑150 (MQ=255)
                                        aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCATCGACAAAAAGGCGCACCACGGAACATGTACAGccc                                                                                                             >  2:1409973/1‑150 (MQ=255)
                                        aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGCGCCCCACCGAACAAGTAAAGccc                                                                                                             >  1:642819/1‑150 (MQ=255)
                                         aGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCg                                                                                                            <  2:4085356/150‑1 (MQ=255)
                                          ggCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAAGCGCACCACCCAACATGTAAAcccccc                                                                                                           >  1:2900344/1‑148 (MQ=255)
                                                  gCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCAAAAAGGCGCACCACGGAACATTTAAAGCCCGCAACCGaag                                                                                                   >  2:85342/1‑148 (MQ=255)
                                                   cTGGTAGCCACACTGTTCATCCTGCCCAGTTGGGTTGGCTGATCCGACTTTTCTTGATGGCCGCCGGCATGCGCTGCCATTTAGCCCTCGCGGTGGTGCGCCTGCTAAAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                  <  2:1306790/150‑1 (MQ=255)
                                                    tGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGGGCCCAAAAAGCGCCCCACCGCACATGTACAACCCCCAAccgcatcg                                                                                                 >  1:257716/1‑143 (MQ=255)
                                                           cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAATGCGCACCCCCGAACATGTACAGCCCCCAACGGATTCGGCTTGAt                                                                                          >  2:536278/1‑150 (MQ=255)
                                                           cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGCCAAAAAGGCGCACCACCCCACATGTAAAGCCCGCAACCGATGCGGCGTAAt                                                                                          >  2:2787428/1‑150 (MQ=255)
                                                            acccggTTCATCCGGCCCTGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCGCGCCGTGGTGCGCCTGCGACAAACAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATg                                                                                         <  2:2347755/145‑1 (MQ=255)
                                                              acTGTTCATCCTGCCCTGTTCGCTTGGCTGATCCGACTGTTCATGTGGGCCGCCGGCATGCGCAGCCTTTGAGCCGGCGCCGGGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcg                                                                                       <  2:3855273/150‑1 (MQ=255)
                                                                 gTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAAgcgcgc                                                                                    >  2:2541952/1‑150 (MQ=255)
                                                                 gTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGAGGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgcgc                                                                                    <  2:1733299/150‑1 (MQ=255)
                                                                 gTTCATCCGGCCCTTTTCGCTTGGCTTATCCGACTTTTCATGTTGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCTACAAAAAGGCGAACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgcgc                                                                                    <  1:2588738/150‑1 (MQ=255)
                                                                  ttCATCCTGCCCTGTTCGGTTGGCTGAGCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCCACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGcc                                                                                   <  1:3358370/150‑1 (MQ=255)
                                                                  ttCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGCCAAAAAGGCGCACCCCCGCACGTGTACAGCCCCCAACCGATGCGGCGTAATGCGCGcc                                                                                   >  2:3395232/1‑150 (MQ=255)
                                                                  ttCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGGCGGCGGCCAAAAAGCGCACCACCGGACGTGTGAACCCCCCCACCGAGGCGGCGGAAGGCGCGcc                                                                                   >  1:3440111/1‑150 (MQ=255)
                                                                      tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCCGCAAAAAAGCGCACCACCCAACCTGTAAAGCCCCCCACCGAAGGGGCGTTATGGGCGcccaat                                                                               >  1:1352498/1‑146 (MQ=255)
                                                                       ccTGCCCAGTTCGCTTGGCTTATCCGACTGTTCATGAGGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATAtt                                                                              <  2:3260711/150‑1 (MQ=255)
                                                                          gCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCCCACCACCCAACATGTAAAACCCCCAACCGAAGCGGCGTAATTCGCGCCAAATTAGc                                                                           >  2:3200009/1‑150 (MQ=255)
                                                                               gTTCGCTTGGCTGATCCGACTGTTCAGGATGGCCGCCCGCAGGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCggg                                                                      <  1:1370586/150‑1 (MQ=255)
                                                                                   gCTTGGCTGATCCGACTGTTCATGCTGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGAc                                                                  <  2:4160951/150‑1 (MQ=255)
                                                                                   gCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCAAAAATGCGCCCCACCGAACGTGTTAAGCCCGCCACCGCAGCGGCTTAATTCGCGCCCAAATAGCGCGGGTggc                                                                  >  1:3265108/1‑148 (MQ=255)
                                                                                    gttGGCTGATCCGACTGTTCATTATGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACt                                                                 <  2:3088945/149‑1 (MQ=255)
                                                                                    gttGGCTGACCGGCTTGTTCCTGTGGGCCGCCGGCATGCGCAGCATTTGCCCCCCCGCGGGGGGGCGCCTGCCAAAAACAGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACt                                                                 <  1:268803/149‑1 (MQ=255)
                                                                                     ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAAGGCGCCCCAACGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCCTATTCGCTCGGGTGACTc                                                                >  2:330435/1‑150 (MQ=255)
                                                                                     ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACCGCACGTGGACAGCCCGCAACGCCTGCGGCGTAATGCGCGCCAAAATCGCGCGGGTGAcgc                                                                >  2:4342383/1‑148 (MQ=255)
                                                                                     ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAACATGTAAACCCCCCAACCGCGGCCGCGTAAATCGCCCCAAAATCGCGCCGGTggaaa                                                                >  2:3910109/1‑146 (MQ=255)
                                                                                       ggCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATGTAAAGCCCGCAAAGGATTCGGGGTAATGCGCCCCCAATTGGCTCGGGTGAcaaaa                                                              >  2:3501630/1‑146 (MQ=255)
                                                                                        gCTGATCGGATGTTTCATGTGGGCGGCGGGCATGCGCAGCCTTTGCGCCCGCGCGGTGGGGGGCCCGCCAAAAAACGCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATg                                                             <  1:3194478/150‑1 (MQ=255)
                                                                                         cTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGCCCAAAAGGCCCCCCACCGACCCTGTTCACCCCGCAAACGAAGCGGCGTAAAGAGCGCCAAATTATCGCGGGTGaaacaatc                                                            >  1:4422411/1‑143 (MQ=255)
                                                                                             tCCGTCTGTTCCTGTGGCCCGCCGGCTGGCGCCGCCTTTGCCCCCGCGCGGTGGTGCGCCTGCCAAAAAAAGCCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                        <  1:1387079/150‑1 (MQ=255)
                                                                                             tCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                        <  2:1284071/150‑1 (MQ=255)
                                                                                              ccGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGCACGTGTACAGCCCCCAACGCATGCGGCGTAATCCGCGCCCGATTAGCGCGGGTGACTCAGGCTGTca                                                       >  2:3180368/1‑150 (MQ=255)
                                                                                               cGACTTTTCATGTTGGCCGCCGGCAGGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCATGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                      <  1:2573940/150‑1 (MQ=255)
                                                                                                        atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGGACTTTTAAAGCCCGCAACGGATGCGGGTTAAAGCGCCCCAAATTCGCTCGGGTGACAGATGCTTTcacgtgcacgtg                                             >  2:2526254/1‑150 (MQ=255)
                                                                                                        atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCTGGGGCCAAAAGGCGCCCCCCCGCACTTTTAAACCCCCCCACGCATTCGGCTTAAATCGCGCCCAAATAGCTCCGGTGAATCATGCATTgaagggcacggg                                             >  2:3791358/1‑148 (MQ=255)
                                                                                                         tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCCCCACGGAACTTGTAAAGCCCCCAACGGGTGCGGCTTAATGCGCGCCAAATTAGCTCGGGTGACTCATGCGTTCacgtgcacgtgc                                            >  2:2246923/1‑150 (MQ=255)
                                                                                                              gccgccGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCTGCGCCAAAAAGGCGCCCCACGGCACATGTAAAGCCCCCAACGGCTGCGGCGTAATGCGCGCCAAAATAGCACGGGGGTCGCATGCGTTAACGTGCACCGGcacaac                                       >  2:2100570/1‑148 (MQ=255)
                                                                                                                cgccgGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCGCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCtt                                     <  2:1129759/150‑1 (MQ=255)
                                                                                                                 gcgggCATGCGCAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:4466031/147‑1 (MQ=255)
                                                                                                                 gcgggCAGGCGCCGCATTTGAGCCCGCGCGGTGGTGCGCCTGCCAAAAAATGCCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:2063360/147‑1 (MQ=255)
                                                                                                                 gccgGCTTGCGCAGCCTTTTCGCCCCCGCCTTGGTGCGCCTGGGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:1622067/150‑1 (MQ=255)
                                                                                                                 gccgGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:4079991/150‑1 (MQ=255)
                                                                                                                 gccgGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  2:3097968/150‑1 (MQ=255)
                                                                                                                   cgGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                  <  2:2058682/150‑1 (MQ=255)
                                                                                                                    gcctggcgcCGCCATTGCGCCCGGGCGGGGGGGCGCATGCCAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTc                                 <  1:1629740/145‑1 (MQ=255)
                                                                                                                     gaaTGCGCAGCCTTTGCGCCCTCGCGTTGGTGCGCCTGCAAAAAAAGGCCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                <  1:1564591/148‑1 (MQ=255)
                                                                                                                        tGCGCAGCATTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGtt                             <  1:2201280/150‑1 (MQ=255)
                                                                                                                         gcgcTGCCAGTGAGCCATCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                            <  2:546763/150‑1 (MQ=255)
                                                                                                                                cATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgcc                     <  2:460141/150‑1 (MQ=255)
                                                                                                                                  ttGCGCCCGGGCGGTGGGGCGCCTGCCAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTc                   <  1:169513/150‑1 (MQ=255)
                                                                                                                                   tGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCg                  <  1:3714892/150‑1 (MQ=255)
                                                                                                                                    gcgCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGc                 <  2:889898/148‑1 (MQ=255)
                                                                                                                                       cccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt              <  2:323206/150‑1 (MQ=255)
                                                                                                                                       cccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt              <  2:3694609/150‑1 (MQ=255)
                                                                                                                                         cTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc            <  2:4223629/150‑1 (MQ=255)
                                                                                                                                          tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt           <  2:1107190/150‑1 (MQ=255)
                                                                                                                                            ggggtggtgCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGt         <  1:3310330/147‑1 (MQ=255)
                                                                                                                                             ccGTGGTGCGCGTGCTAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGtt        <  2:391608/150‑1 (MQ=255)
                                                                                                                                                 ggtgCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc    <  2:960706/150‑1 (MQ=255)
                                                                                                                                                   tgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACct  >  1:4173345/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCT  >  NC_002947/3778191‑3778487

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: