Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,371 T→G 17.2% intergenic (+61/‑82) 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,781,3710TG17.2% 353.2 / 9.1 184intergenic (+61/‑82)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/110);  new base G (29/5);  total (68/116)
Fisher's exact test for biased strand distribution p-value = 3.71e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

TCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGC  >  NC_002947/3781222‑3781520
                                                                                                                                                     |                                                                                                                                                     
tCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCCTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCt                                                                                                                                                       <  2:7309133/150‑1 (MQ=255)
  aGCATCTGGGAGGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCCGc                                                                                                                                                     >  2:2027838/1‑150 (MQ=255)
  aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                     <  2:1474710/150‑1 (MQ=255)
  aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                     <  2:6291632/150‑1 (MQ=255)
  aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGATTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                     <  2:4995896/150‑1 (MQ=255)
   gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCAGTGAGCCCTCGCCGTGGGGGGCCGGCg                                                                                                                                                    >  1:3363/1‑150 (MQ=255)
     aTCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCggc                                                                                                                                                  >  1:6778332/1‑148 (MQ=255)
     aTCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCACAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGAc                                                                                                                                                  >  1:1075830/1‑150 (MQ=255)
         gggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCCGCGACaaaa                                                                                                                                              >  2:1538485/1‑150 (MQ=255)
          ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCaaaaa                                                                                                                                             >  1:7062805/1‑150 (MQ=255)
           gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAg                                                                                                                                            >  1:4853140/1‑150 (MQ=255)
            aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                           <  1:6727431/150‑1 (MQ=255)
             aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                          <  2:3128731/150‑1 (MQ=255)
               ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                        <  1:7093057/150‑1 (MQ=255)
               ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                        <  2:6131025/150‑1 (MQ=255)
               ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACCAAAAGgcgc                                                                                                                                        >  2:3827874/1‑150 (MQ=255)
                cGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAAAGCGca                                                                                                                                       >  1:6217799/1‑150 (MQ=255)
                cGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACCAAAAAgcgcc                                                                                                                                       >  1:4436856/1‑149 (MQ=255)
                 gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGACAAAAAAgcgccc                                                                                                                                      >  1:4761838/1‑148 (MQ=255)
                 gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                      >  2:4314039/1‑150 (MQ=255)
                 gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGGCTGCGACCAAAAGGCGcac                                                                                                                                      >  1:2441623/1‑150 (MQ=255)
                      gccgccCCAGAGCCAAGGCTTCTTGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACgg                                                                                                                                 <  1:1856289/150‑1 (MQ=255)
                       ccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGGCGGCGACCAAAAGGCGccccaccgc                                                                                                                                >  1:5857033/1‑147 (MQ=255)
                             cAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGGTTGGCTGATCCGAGTGTTCATTATGGCCGCGGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGt                                                                                                                          <  2:5083199/150‑1 (MQ=255)
                             cAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGaaatagt                                                                                                                          >  1:7301279/1‑145 (MQ=255)
                              agagCCATGGATTCATGGTGGTAGGCACAGTGTTACTCCTGCCCCTTTCGCTTGGCTGATCCGACTTTACAGGAGGGCCGCCGGCATGCGCAGCCTTTGAGCCCGCGCGTTGGTGCGCCTGCCACAAAACGGCGCCCCACGGAACTTGTa                                                                                                                         <  2:4646623/150‑1 (MQ=255)
                              agagCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTa                                                                                                                         <  2:6157943/150‑1 (MQ=255)
                                    aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGGGGCAAAATGGGGCCCCACGGAACTTGTTCAGccc                                                                                                                   >  2:1682182/1‑150 (MQ=255)
                                    aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATGTACAGccc                                                                                                                   >  2:1091762/1‑150 (MQ=255)
                                         ttCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGGCAAAAAGGCGCACCACGGAACGTGTACAGCCCGCAAc                                                                                                              >  1:3031410/1‑150 (MQ=255)
                                            aTGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAAGCGCACCACGCAACAATTACAGCCCGCAAAGGa                                                                                                           >  2:2392180/1‑150 (MQ=255)
                                               cTGGTTGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCTTGTTGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                        <  2:5387268/150‑1 (MQ=255)
                                                 ggTAGGCACCGTGTTCACCCTCCCCAGTTCGCTGGGCTGACCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGACCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                                      <  2:4635608/150‑1 (MQ=255)
                                                 ggTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                                      <  2:7403259/150‑1 (MQ=255)
                                                      gCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACATGTACAGCCCGCAACCGATGCGGCTTaa                                                                                                 >  2:2928276/1‑150 (MQ=255)
                                                       cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCTTGTTGGCCGCCGGATTGCGCCGCCATTGCGCCCCCGCGGGGGTGCGCCTGCTACAAAATGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                <  1:2299921/150‑1 (MQ=255)
                                                        acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACCGAACATGTACAGCCCGCAACCGATGCGGCTTAATg                                                                                               >  2:4814878/1‑150 (MQ=255)
                                                        acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                       <  1:895880/110‑1 (MQ=255)
                                                        acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                       >  2:895880/1‑110 (MQ=255)
                                                         cacTGTTCATCCTGCCCTGTTCGATTGGCTGATCCGACTGTTCATGATGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgc                                                                                              <  2:1823422/150‑1 (MQ=255)
                                                          acTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACCGAACATGTAAAGCCCCCAACCGATTCGGCTTAAAgcg                                                                                             >  2:153898/1‑150 (MQ=255)
                                                           cggTTCATCCGGCCCAGTTCGGTGGGCTGATCCGACTGTTCTTTATGGCCGCGGGCATGCGCTGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                            <  2:6516643/148‑1 (MQ=255)
                                                                cATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCACAAAAGCCCCCCCCCGCACGTGTTAAGCCCGCCACCGATGCGGCGGAAATCGCGccca                                                                                       >  1:2188379/1‑148 (MQ=255)
                                                                 aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAAAGCGCCCCCCCGAACTTTTAAAGCCCGCCACCGGTGCGGCTTAATGCGCGcccaa                                                                                      >  1:5477360/1‑147 (MQ=255)
                                                                   ccTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCCACCGATGCGGCTTAATGCGCGCCCTAtt                                                                                    >  2:4732963/1‑150 (MQ=255)
                                                                     tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGGCTGGGTGAAAAATGGGCGCCACGGGAATTTTAAAGGCCGCAAAGGTTTGGGGTTTTTGGGCGGCATAtttg                                                                                  >  2:971041/1‑148 (MQ=255)
                                                                     tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaa                                                                                                                                              >  2:6605886/1‑90 (MQ=255)
                                                                     tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaa                                                                                                                                              <  1:6605886/90‑1 (MQ=255)
                                                                         cAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCTGCGGCAAAAAGGCGCACCCCCGAACATGTAAAGCCCCCCACCCCTGCGGCGTAATGCGCGCCAGATTAGcgag                                                                              >  1:5255696/1‑147 (MQ=255)
                                                                                 gtGGCTGATCCTACTGTTCATTTTGGCCGCCGGAATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                      <  2:5998767/149‑1 (MQ=255)
                                                                                   ggCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAt                                                                    <  2:5768634/150‑1 (MQ=255)
                                                                                         tCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                              <  2:3611880/150‑1 (MQ=255)
                                                                                           cGAGTGTTCCTGCTGGCCGCCGGCATGCGCAGCCATTGCGCCCCCGCCGTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                            <  2:3890851/150‑1 (MQ=255)
                                                                                             agtGTTCATGATGGCCGCCGGCTTGCGCCGCATTGGCGCCCTCGCCTTGGTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgt                                                          <  1:550509/148‑1 (MQ=255)
                                                                                             aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGAACAAAATGCGCCCCACCGAACGTTTACAGCCCCCAACGGATGCGGCGTAATTCGCGCCAGATTAGGGCGGGGGACGCAAGCGTTcaagt                                                          >  2:2189387/1‑150 (MQ=255)
                                                                                               tGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGCGCCCCCCGGCACTTGTACAGCCCCCCCCCGCTGCGGCGTAAAGAGCCCCCCATTAGCGCGGGTGCCACCTGCTTTaacctgc                                                        >  2:5103865/1‑150 (MQ=255)
                                                                                                 ttCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCCCCCCCGAACGTGTACACCCCCCCACCGCTGCGGCGTAATGCGCGCCAAAATCGCTCGGGGGCCGCATGCGTTaacctgccc                                                      >  2:4174343/1‑148 (MQ=255)
                                                                                                    ttgttgGCCGCGGGCATGCGCCGCATTTGCGCCCTCGCGGGGGGGCGCCGCCCACAAAAAGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                   <  1:1334256/146‑1 (MQ=255)
                                                                                                    atgatgGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAACGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                   <  2:5095761/150‑1 (MQ=255)
                                                                                                    atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAACGCGCCCCCCCGCACTTGTAAACCCCCCACACGCTGCGGCTTAAAGCGCGCCAAATTCGCTCGGGGGGCGCATGCTTTaaactgcacgtg                                                   >  2:4349396/1‑150 (MQ=255)
                                                                                                     tgtgggCCGCGGGCATGCGCACCCATTGCCCCCGCGCGGTGGTGCGCCTGCTAAAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                                  <  1:1589008/146‑1 (MQ=255)
                                                                                                      gagggCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCAcgtgcacgtgca                                                 <  2:314092/147‑1 (MQ=255)
                                                                                                         ggCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGACAAAAAGGCGCACCACGGAAAGTTTAAAGCTAGCCAACGAAGCGGCTTAATTCGCGCCCAATTAGCGCGGGTGACTCAAGCTTTaacgtgcacgtgagttt                                              >  2:3370813/1‑145 (MQ=255)
                                                                                                         ggCCGCCGGCAGGGGCCGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACAt                                              <  2:1435505/150‑1 (MQ=255)
                                                                                                           ccgccgTCTTGCGCAGATTTTGCGCCCGCGCGGTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCt                                            <  1:3364999/150‑1 (MQ=255)
                                                                                                           ccgccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtgaaagcat                                            <  2:5114041/150‑8 (MQ=255)
                                                                                                            cgccgGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCGTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCtt                                           <  2:2722357/150‑1 (MQ=255)
                                                                                                             gccgGCTTGCGCAGCCTTTGAGCCCTCGCGGTGGTGGGCGTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  1:5389779/150‑1 (MQ=255)
                                                                                                             gccgGCTTGCGCAGCCATTGTGCCCTCGCCGTGGTGCGCATGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  2:521521/150‑1 (MQ=255)
                                                                                                             gccgGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  2:7039167/150‑1 (MQ=255)
                                                                                                             gccgGCATGCGCCGCCATTGCGCCCTCGCCGTGGGGCGCGGGTGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  2:2403507/150‑1 (MQ=255)
                                                                                                             gccgGCATGCGCAGCAATTGCGCCCTCGCGGTGGTGCCCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  2:277129/150‑1 (MQ=255)
                                                                                                               cgGCTTGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  1:4818482/150‑1 (MQ=255)
                                                                                                               cgGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  1:1225588/150‑1 (MQ=255)
                                                                                                               cgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  2:1209417/150‑1 (MQ=255)
                                                                                                                 gCATGCGAAGCATTTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                      <  1:1215627/150‑1 (MQ=255)
                                                                                                                     gcgcCGCCCTTGCGCCCTCGCCGTGGTGCGCCTGCGCCAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                  <  2:5702982/150‑1 (MQ=255)
                                                                                                                     gcgcAGCCTTTGCGCCCTCGCCGGGGGGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                  <  2:5702903/150‑1 (MQ=255)
                                                                                                                       gcAGCCTTTGCGCCCTCGGGGTGGTGCGCCTGCGAAAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgc                                <  1:1445657/150‑1 (MQ=255)
                                                                                                                       gcAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgc                                <  2:4958865/150‑1 (MQ=255)
                                                                                                                                     cgcGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                  <  2:3311314/148‑1 (MQ=255)
                                                                                                                                     cgcGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                  <  2:1938871/148‑1 (MQ=255)
                                                                                                                                     cTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                  <  2:4022837/150‑1 (MQ=255)
                                                                                                                                     cTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                  <  2:534899/150‑1 (MQ=255)
                                                                                                                                         ccGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGtt              <  2:6021485/150‑1 (MQ=255)
                                                                                                                                          cGTGGTGGGCCTTCTACAAATAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTc             <  2:3235659/150‑1 (MQ=255)
                                                                                                                                          cGTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCATTGTCTGTTCGCCTGCCTCGCCGTGCTGTTc             <  2:3368125/150‑1 (MQ=255)
                                                                                                                                           gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa            >  1:3000716/1‑150 (MQ=255)
                                                                                                                                             gtggcgcATGCAAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc          <  1:314414/147‑1 (MQ=255)
                                                                                                                                             ggtgCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc          <  2:4393134/150‑1 (MQ=255)
                                                                                                                                             ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc          >  2:4879857/1‑150 (MQ=255)
                                                                                                                                             ggtgCGCCTGCCACAAAATGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc          <  1:3827874/150‑1 (MQ=255)
                                                                                                                                             ggtgCGCCTCCAAAAAAATGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc          <  1:115080/150‑1 (MQ=255)
                                                                                                                                               tgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          >  2:7239932/1‑116 (MQ=255)
                                                                                                                                                gcgcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGAGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctg       >  1:196541/1‑150 (MQ=255)
                                                                                                                                                 cgcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgc      >  2:2736157/1‑150 (MQ=255)
                                                                                                                                                 cgcCGGCGACAAATAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTTCCTCGCCGTGCTGTTCAACctgc      <  2:5490932/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgct     <  1:4283517/150‑1 (MQ=255)
                                                                                                                                                    cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  1:3267752/1‑150 (MQ=255)
                                                                                                                                                    cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  2:6082118/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGC  >  NC_002947/3781222‑3781520

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: