Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,371 T→G 14.3% 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,3710TG14.3% 342.0 / 4.1 181intergenic (+61/‑82)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/92);  new base G (26/4);  total (84/97)
Fisher's exact test for biased strand distribution p-value = 1.41e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

AGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGC  >  NC_002947/3781224‑3781520
                                                                                                                                                   |                                                                                                                                                     
aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                     <  2:2428033/150‑1 (MQ=255)
 gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCg                                                                                                                                                    >  1:553502/1‑150 (MQ=255)
 gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCg                                                                                                                                                    >  1:645606/1‑150 (MQ=255)
 gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCg                                                                                                                                                    >  2:2429374/1‑150 (MQ=255)
    tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGaga                                                                                                                                                 >  1:2395618/1‑148 (MQ=255)
        ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCaaaaa                                                                                                                                             >  1:84077/1‑150 (MQ=255)
           aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                          <  1:6778617/150‑1 (MQ=255)
            gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                         <  2:4705786/150‑1 (MQ=255)
             ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                        <  2:1804321/150‑1 (MQ=255)
               gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                      >  2:6297695/1‑150 (MQ=255)
               gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAAGgcgccc                                                                                                                                      >  1:2975187/1‑148 (MQ=255)
               gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGGCAAAAAGGCGcac                                                                                                                                      >  2:684482/1‑150 (MQ=255)
                 aGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcacca                                                                                                                                    <  1:4152624/150‑1 (MQ=255)
                    gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACCAAAAAGCGcaccaccg                                                                                                                                 >  2:1186635/1‑148 (MQ=255)
                    gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAACGCGcccccacg                                                                                                                                 >  1:2575766/1‑146 (MQ=255)
                       gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGacc                                                                                                                              >  2:5012301/1‑148 (MQ=255)
                       gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGCCAAAAAGGCGCACCACCGAAc                                                                                                                              >  2:3318446/1‑150 (MQ=255)
                       gccCCAGAGCCAAGGCATCATGCAGGTAGGCACACTTTACATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATAGAGCCCTCGCCGTGGGGCGCCAGCGAAAAAAAGGCGCACCACGGAAc                                                                                                                              >  2:3318514/1‑150 (MQ=255)
                          ccAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACAAAAAAGCGCACCACCGAACGTg                                                                                                                           >  1:4817363/1‑150 (MQ=255)
                           cAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAcaagt                                                                                                                          >  2:3664763/1‑146 (MQ=255)
                                ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCCCGGAACATGTACAGc                                                                                                                     >  1:2805365/1‑150 (MQ=255)
                                ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCACCGAACTTGTACAGc                                                                                                                     >  1:6778474/1‑150 (MQ=255)
                                  aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGccc                                                                                                                   >  2:2768242/1‑150 (MQ=255)
                                       ttCATGCTGGTTGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTTATCCGAGTGTTCATGTTGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAAc                                                                                                              <  2:1123866/150‑1 (MQ=255)
                                        tCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGGCAAAAAAGCGCACCACGCAACGAGTACACCCCGCAACg                                                                                                             >  2:1302005/1‑150 (MQ=255)
                                             cTGGTAGGCACACTTTACATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCCCAGCCATTGAGCCCACGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                        <  2:5651183/150‑1 (MQ=255)
                                                    gCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGCGCACCCCGGCCCTTGTAAAGCCCGCCACGGGTGCGGcggga                                                                                                 >  2:4037812/1‑146 (MQ=255)
                                                     cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGGACATTTAAAGCCCGCAAACGGTTCGGCGTAAt                                                                                                >  2:1189352/1‑150 (MQ=255)
                                                     cacaCTGTTCATCCGGCCCTGTTCGGTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCTGGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                <  2:5919181/150‑1 (MQ=255)
                                                     cacaCTGTTCAGCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCAGGCTGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                <  2:2242925/150‑1 (MQ=255)
                                                      acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATg                                                                                               >  2:6768032/1‑150 (MQ=255)
                                                      acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGGCGGCGGCCAAAACGCGCCCCCCCGGACGGGGACAGCCCCCCACCGCTGCGGcggaaag                                                                                               >  1:6529562/1‑144 (MQ=255)
                                                        acTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGGCAAAAATGCGCAACACGGAACGTGTAAAACCCCCAAAGGAATCGGGTTAAttag                                                                                             >  1:3285582/1‑147 (MQ=255)
                                                          tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAAGCGCACCACCGAACATGTACAGCCCCCAACCGATGCGGCTTAAAtcgcg                                                                                           >  2:958855/1‑150 (MQ=255)
                                                          tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAAATGTACAGCCCGCAACGGATGCGGCTTAATgcgcg                                                                                           >  2:3866792/1‑150 (MQ=255)
                                                           gTTCAGCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgcgc                                                                                          <  2:6731676/150‑1 (MQ=255)
                                                               acccggCCCTGTTCGCTGGGCTGAGCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCAATTGAGCCCGCGCGGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCata                                                                                      <  2:6396377/145‑1 (MQ=255)
                                                               aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGGCAAAAAGGCGCACCACGGGACTTGTACAGCCCGCAACCGATGCGGCTTAATGCGCGCCaat                                                                                      >  2:554356/1‑148 (MQ=255)
                                                                tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACCAAAAGGCGCACCACGGAACATGTACAGCCCGCCACCGATGCGGCGTAATGCGCGCCatat                                                                                     >  1:6170723/1‑150 (MQ=255)
                                                                  cTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGGCAAAAAGGCGCACCACCGCACTTGTAAAGCCCCCCACCGATGCGGCGTAATGCGCGCCCAATTc                                                                                   >  2:365816/1‑150 (MQ=255)
                                                                   tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCCTATTCg                                                                                  >  2:6535870/1‑150 (MQ=255)
                                                                   tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCCCCCAAAACGCGCCCCCCCGCCCCGGTACAGCCCCCCACCGCGGCGGCGTAAGGCGCGCCACAGTCg                                                                                  >  2:4099340/1‑150 (MQ=255)
                                                                      ccttttCGCTTGGCTTATCCGAGTGTTCCTGTGGGCCGCCGGCAGGCGCAGCCTTTGCGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                               <  1:4657699/146‑1 (MQ=255)
                                                                      cccgTTCGCTTGGCTGAGCCTACTTTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                               <  2:4181602/147‑1 (MQ=255)
                                                                      ccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                               <  2:6285884/150‑1 (MQ=255)
                                                                       cAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAAGCCCCCCACCGCACATGTACAGCCCGCCACCGATGCGGCGGAAAGCGCGCCCAATTAGCTCg                                                                              >  1:3873134/1‑150 (MQ=255)
                                                                        ttttCGCTTGGATTATCCTACTTTTCTTGTTGGCCGCGGGCTTGCGCAGCATTTGAGCCCTCGCGGTGGTGCGCCTGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCgg                                                                             <  1:3246548/148‑1 (MQ=255)
                                                                        aGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGACAAAAAAGCGCCCCACCGAACTTGTACAGCCCCCAACGGAGGCGGCGTAATGCGCGCCCTATTAGGTCgg                                                                             >  1:2699179/1‑150 (MQ=255)
                                                                         tttCGCTTGGCTGATCCGACTGTCATTGATGGCCGCGGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCggg                                                                            <  2:2938275/149‑1 (MQ=255)
                                                                          ttCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAAGGGGCCCCAACGAACGAGTACAGCCCGCAACCGATGCGGCTTAATTCGCGCCCTATTCGCACGGGt                                                                           >  2:722941/1‑150 (MQ=255)
                                                                           tCGGTTGGCTGACCCTACTGTTCCTGTGGGCCGCCGGCATGCGCAGCATTTGCGCCCCCGCGGTGGGGGGCCTGCCACAAAACGCCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTg                                                                          <  1:5332953/150‑1 (MQ=255)
                                                                               ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAAAGGATGCGGCTTAATGCGCGCCAAATTCGCTCGGGTGACTc                                                                      >  2:526823/1‑150 (MQ=255)
                                                                               tgggCTGTTCCGCCTGTTCCTGAGGGCCGCCGGCTTGCGCAGCCTTTGCGCCCTCGCGGTGGGGCGCCTGCGCCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                      <  2:6428565/148‑1 (MQ=255)
                                                                                 gggtgttCCGAGTGTTCCTGCGGGCCGCCGGCATGCGCCGCAATTGCCCCCGGGCGGTGGTGCGCCCGCCAAAAAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAt                                                                    <  1:2540838/144‑1 (MQ=255)
                                                                                 ggCTGACCCTCCTGTTCTTTTGGGCCGCCGGCAGGCGCCGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAt                                                                    <  2:2850703/150‑1 (MQ=255)
                                                                                   cTGATCCGACTGTTCATGTTGGCCGCCGGCTTGCGCAGCCATTGCGCCCGCGCGGTGGGGCGCCTGCGACAAAAAGCCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGc                                                                  <  1:4539690/150‑1 (MQ=255)
                                                                                    ggATCCGACTGTTCATGCGGGCCGCCGGCATGCGCAGCCATTGAGCCCGCGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCt                                                                 <  2:6036892/149‑1 (MQ=255)
                                                                                         cgccTGTTCATGTTGGCCGCCGGCAGGCCCAGCCATTGAGCCCGCGCCGGGGGGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                            <  2:325056/147‑1 (MQ=255)
                                                                                         cGACTGTTCTTGTGGGCCGCGGGCCGGCGCCGCCTTTGCCCCCCCGCGGGGGGGCGCCGGCCACAAACCGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                            <  1:3159216/150‑1 (MQ=255)
                                                                                               ttCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAAAATGAAAAGACCGCAACGGAAGCCCCTTAAAGCGCGCCAAATTAGCTCGGGTGACCAATGCATAcaccagcac                                                      >  2:6263151/1‑150 (MQ=255)
                                                                                               ttCATGATGGCCGCCGGAATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcac                                                      <  2:1696620/150‑1 (MQ=255)
                                                                                                tCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCACCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacg                                                     >  1:1830847/1‑150 (MQ=255)
                                                                                                 cATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCCAAAAGGCGCCCCCCCGGCCTAGTAAAGCCCCCAACGGCTGCGGCGTAATGAGCGCCAAAAAAGCTCGGGGGGCGCATGCGTTgacgggcaccg                                                    >  2:2391283/1‑148 (MQ=255)
                                                                                                  atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                   <  2:1638456/150‑1 (MQ=255)
                                                                                                  atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCCCCACCGAACTTTTACAGCCCCCAACCGATGCGGCTTAATGCGCGCCAAATTCGCTCGGGTGTCTCATTCTTTaacgtgcacgtg                                                   >  1:4506306/1‑150 (MQ=255)
                                                                                                   tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTTCGCCTTCGGCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCCGCTTTATGCGCGCCATATTCGCTCGGGTGACTCATGCGGTCacgtgcacgtgc                                                  >  2:5701305/1‑150 (MQ=255)
                                                                                                     atgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACgtgcacgtgcac                                                <  2:6652914/150‑1 (MQ=255)
                                                                                                         ccgccgGCTTGCGCAGCATTTGAGCCCTGGCGGTGTTGCGCCTGCTAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCt                                            <  1:1531081/150‑1 (MQ=255)
                                                                                                          cgccgGCATGCGCAGCAATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCtt                                           <  2:1384629/150‑1 (MQ=255)
                                                                                                           gcgggCATGCGCAGCCTTTGAGCCCTCGCGGGGGTGCGCCTGCGAAAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                          <  1:6665950/147‑1 (MQ=255)
                                                                                                           gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCCACAAAAAGGCGCCCCACCGAACATTTAAAGCCCGCCACGGGTGCGGCGTAATTCGCGCCCAATTAGCGCGGGGGTCTCAAGCGTTCAAGGGCACGGGCACAACttt                                          >  2:2365827/1‑150 (MQ=255)
                                                                                                             cggctggcgcAGCCTTTGACCCCTCGCGGGGGTGCGCCTGCGAAAAACCGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  1:2029030/144‑1 (MQ=255)
                                                                                                             cggctggcgcAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  2:5211215/144‑1 (MQ=255)
                                                                                                             cgGCATGCGCCGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                        <  1:5624061/150‑1 (MQ=255)
                                                                                                             cgGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGGGCACCACCGGACATGTACAGCCCCCCACGGCTGCGGCTTAATGCCCGCCATATTCGCTCGGGTGCATCCGGCGTTCACGTGCACGTGCACATCTTTGt                                        >  2:4394348/1‑150 (MQ=255)
                                                                                                               gcttGCGCTGCATTTGAGCCCTGGCGGTGGGGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                      <  1:2065832/147‑1 (MQ=255)
                                                                                                                   gcgcCGCCATTGAGCCCGCGCCGTGGTGCGCCTGCTAAAAAAAGGCGAACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                  <  2:1840756/150‑1 (MQ=255)
                                                                                                                   gcgcAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                  <  2:255003/150‑1 (MQ=255)
                                                                                                                   gcgcAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAACGTGTACAGCCCCCAACCGATTCGGCTTAAAGCGCGCCAAATTCGCTCGGGTGGCTCATGCTTTAACGTGCACGTGCACAAATTTTTCTGTTc                                  >  2:2161838/1‑150 (MQ=255)
                                                                                                                       agcatttGCGCCCTCGCGGTGGTGCGCCTGCCACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcct                              <  1:3664763/145‑1 (MQ=255)
                                                                                                                          ctttGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCGCGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgcc                           <  1:5391071/148‑1 (MQ=255)
                                                                                                                               cgCCCTCGCGGTGGTGCGCCTGCTAAAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                      <  1:1796599/149‑1 (MQ=255)
                                                                                                                               aGCCCTCGCCGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                      <  2:4610695/150‑1 (MQ=255)
                                                                                                                                 cccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt                    <  2:3931456/150‑1 (MQ=255)
                                                                                                                                 cccTCGCCGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt                    <  1:5546765/150‑1 (MQ=255)
                                                                                                                                  ccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTg                   <  2:3065100/150‑1 (MQ=255)
                                                                                                                                    tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt                 <  2:2631755/150‑1 (MQ=255)
                                                                                                                                      gCCGTGGTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGt               <  1:949263/150‑1 (MQ=255)
                                                                                                                                         gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa            <  1:5061231/150‑1 (MQ=255)
                                                                                                                                          tggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCaa           <  2:2716931/150‑1 (MQ=255)
                                                                                                                                              gcgcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctg       >  2:833114/1‑150 (MQ=255)
                                                                                                                                                 ccTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCGTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgctg    >  2:2347275/1‑150 (MQ=255)
                                                                                                                                                  cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   <  2:6309929/150‑1 (MQ=255)
                                                                                                                                                  cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  2:6107627/1‑150 (MQ=255)
                                                                                                                                                  cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  2:5027782/1‑150 (MQ=255)
                                                                                                                                                  cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  1:4066798/1‑150 (MQ=255)
                                                                                                                                                  cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg   >  2:3817131/1‑150 (MQ=255)
                                                                                                                                                   tGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGc  >  2:6265515/1‑150 (MQ=255)
                                                                                                                                                   tGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGc  >  2:5158543/1‑150 (MQ=255)
                                                                                                                                                   tGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGGTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGc  >  2:864625/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGC  >  NC_002947/3781224‑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: