Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,810,377 T→G 8.5% intergenic (+96/‑73) eno → / → ftsB enolase/cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,810,3770TG 8.5% 538.9 / 11.3 192intergenic (+96/‑73)eno/ftsBenolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (79/89);  new base G (15/6);  total (97/95)
Fisher's exact test for biased strand distribution p-value = 3.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

AGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGTAA  >  NC_002947/1810229‑1810526
                                                                                                                                                    |                                                                                                                                                     
agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCTc                                                                                                                                                      >  1:4286864/1‑150 (MQ=255)
 gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCg                                                                                                                                                     >  2:408362/1‑150 (MQ=255)
 gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCg                                                                                                                                                     <  1:1583097/150‑1 (MQ=255)
 gCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCg                                                                                                                                                     <  1:3641109/150‑1 (MQ=255)
  cAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCgg                                                                                                                                                    >  1:1504687/1‑150 (MQ=255)
   aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                                                                   <  1:3690829/150‑1 (MQ=255)
   aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCGCGGc                                                                                                                                                   >  1:64704/1‑150 (MQ=255)
     cTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCa                                                                                                                                                 <  1:1063784/150‑1 (MQ=255)
     cTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCa                                                                                                                                                 <  1:3437354/150‑1 (MQ=255)
       gggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAgg                                                                                                                                               <  1:4132913/150‑1 (MQ=255)
        ggCGCCACAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGc                                                                                                                                              <  2:2817641/150‑1 (MQ=255)
        ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGc                                                                                                                                              >  1:2207639/1‑150 (MQ=255)
         gcgcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTTCTGTACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                                                                             <  1:2656726/150‑1 (MQ=255)
           gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTc                                                                                                                                           <  2:4155958/150‑1 (MQ=255)
           gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTc                                                                                                                                           <  1:3940465/150‑1 (MQ=255)
           gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTc                                                                                                                                           <  1:1212620/150‑1 (MQ=255)
            ccAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCg                                                                                                                                          >  2:2559801/1‑150 (MQ=255)
             cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt                                                                                                                                         >  1:1270651/1‑150 (MQ=255)
              aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTg                                                                                                                                        >  2:3710774/1‑150 (MQ=255)
              aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTCGtc                                                                                                                                        <  2:2318083/150‑2 (MQ=255)
                aGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGcg                                                                                                                                      >  1:1773979/1‑149 (MQ=255)
                 gCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTa                                                                                                                                     <  2:1907220/150‑1 (MQ=255)
                   ggTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTAtt                                                                                                                                   >  2:1783435/1‑150 (MQ=255)
                    gTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGCGGCCCGGCTTCGGGGTATTc                                                                                                                                  >  1:188155/1‑150 (MQ=255)
                      ttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCGGGCCGCGGCCCGGCGTCGGGCTATTCGa                                                                                                                                >  1:32432/1‑150 (MQ=255)
                          cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGcc                                                                                                                            >  2:1139147/1‑150 (MQ=255)
                                   gCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                                                                   >  2:3508808/1‑150 (MQ=255)
                                    ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAATCCAGGCCTCGGCCAGGCTTAGTGCTATTAGAGGCCCCCAACTaaa                                                                                                                  >  1:3072935/1‑149 (MQ=255)
                                    ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAAc                                                                                                                  <  2:3850055/150‑1 (MQ=255)
                                     cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                                                                       >  1:4511922/1‑112 (MQ=255)
                                     cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                                                                       <  2:4511922/112‑1 (MQ=255)
                                     cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCGAGGCCAGGCGTCGTTCTAATACATGCCCCCAACTaaaa                                                                                                                 >  1:46591/1‑148 (MQ=255)
                                      gAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                                >  2:463902/1‑150 (MQ=255)
                                      gAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                                <  2:1265785/150‑1 (MQ=255)
                                            cgcgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgc                                                                                                          <  1:1152632/150‑1 (MQ=255)
                                            cgcgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAcgc                                                                                                          <  1:4239740/150‑1 (MQ=255)
                                              cgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                                                                        <  2:4490539/150‑1 (MQ=255)
                                                 cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAg                                                                                                     <  2:3121727/150‑1 (MQ=255)
                                                 cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAg                                                                                                     <  2:2499155/150‑1 (MQ=255)
                                                   aaGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGcc                                                                                                   >  2:1298750/1‑150 (MQ=255)
                                                      cAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttt                                                                                                >  2:627062/1‑150 (MQ=255)
                                                      cAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttt                                                                                                <  1:2544402/150‑1 (MQ=255)
                                                       aaGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCtttt                                                                                               <  1:631924/150‑1 (MQ=255)
                                                              ggTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGTCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                                                                        <  2:3875212/150‑1 (MQ=255)
                                                              ggTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                                                                        <  2:3024587/150‑1 (MQ=255)
                                                               gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAAGCGGCAGCCTTTTTAAACTgg                                                                                       >  1:3469097/1‑150 (MQ=255)
                                                               gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTGTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTgg                                                                                       >  2:2599382/1‑150 (MQ=255)
                                                                   aaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATaa                                                                                   <  1:895357/150‑1 (MQ=255)
                                                                   aaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACTGGATaa                                                                                   >  2:2001811/1‑150 (MQ=255)
                                                                        cagcagCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACAGGATAACGTGa                                                                              >  2:1531619/1‑150 (MQ=255)
                                                                          gcagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATg                                                                            <  1:2454290/150‑1 (MQ=255)
                                                                           cagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc                                                                           <  2:142626/150‑1 (MQ=255)
                                                                           cagcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgc                                                                           >  2:3802426/1‑150 (MQ=255)
                                                                            agcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATgcg                                                                          >  1:310411/1‑150 (MQ=255)
                                                                                ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCCGGGTTCGTGCGATTCGAGGCCCCCAAGTAAAAACTACGCCGCACCCTTTTTAACCTGGGTAACGTGATGCGCAGt                                                                      >  1:2348108/1‑150 (MQ=255)
                                                                                ggAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCACGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCACCCGTTTTAACCTGGGTAACGTGATGCGCAGt                                                                      >  1:572426/1‑150 (MQ=255)
                                                                                    ttGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCt                                                                  >  2:1870255/1‑150 (MQ=255)
                                                                                      gTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                                                                <  1:2328921/150‑1 (MQ=255)
                                                                                      gTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACCAATCCTGGCCTCGGCCAGGCTTCGTGCTAATCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGAGCAGTACCTAt                                                                >  2:2070993/1‑150 (MQ=255)
                                                                                       tAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCCAAGTAACAAGTACGCGGCCGCCTTTTTAACCTGGATAAAGTGATGCGCAGTcccgctt                                                               >  2:2263366/1‑146 (MQ=255)
                                                                                        aGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                              <  2:1237195/150‑1 (MQ=255)
                                                                                        aGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                              <  1:2032235/150‑1 (MQ=255)
                                                                                        aGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                              <  1:3321004/150‑1 (MQ=255)
                                                                                          gAACGTTCGTACTGACCTGTCCTACGTACTAACGGGTGTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgt                                                            <  1:2385627/150‑1 (MQ=255)
                                                                                          gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCACGGCCCGGCGTCGTGGTATTCGAGGCCCCCAAGTAACAAGTACGCGGCAGCCTTTTTTAACTGGATTACGTGGTGCGCACTCCCTATTGgt                                                            >  1:4270278/1‑150 (MQ=255)
                                                                                          gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCTCGGCCAGGCTTCGTGGGATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACTGGATAACGTGATGCGCCGTCCCGATTGgt                                                            >  1:268892/1‑150 (MQ=255)
                                                                                            aCGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgttg                                                          <  2:439098/150‑1 (MQ=255)
                                                                                               ttCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTc                                                       <  1:4391874/150‑1 (MQ=255)
                                                                                                     cTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCGATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTGCTCGtt                                                 >  2:3819361/1‑150 (MQ=255)
                                                                                                      tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCGATTGGTTGTTCCTCGttg                                                >  2:673256/1‑149 (MQ=255)
                                                                                                      tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACTGGATAACGTGATGCGCAGTCCCGATTGGTTGTTCCACGtta                                                >  1:337798/1‑149 (MQ=255)
                                                                                                        aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGGTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAAACTGGATAACGTGATGCGCAGTCCCGATTGGTTGTTACACGTTATg                                              >  1:2887166/1‑150 (MQ=255)
                                                                                                        aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCGTCGTGCGATTCGAGGCCCCGAAGTAACAAGGACGCGGCCGCCGTTTTAACCTGGAGAACGTGATGCGCACTCCCGATTGGTTGTTCCTCGATCTg                                              >  1:2974176/1‑150 (MQ=255)
                                                                                                            ttCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGctct                                          <  2:2186900/150‑1 (MQ=255)
                                                                                                               cTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGc                                       <  1:2794349/150‑1 (MQ=255)
                                                                                                                tACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCt                                      >  2:3666362/1‑150 (MQ=255)
                                                                                                                tACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCt                                      <  2:4286864/150‑1 (MQ=255)
                                                                                                                   gTACTAACGGGTGTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGtt                                   <  2:3094372/150‑1 (MQ=255)
                                                                                                                         aCGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTgg                             <  1:419018/150‑1 (MQ=255)
                                                                                                                          cGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGc                            <  1:199061/150‑1 (MQ=255)
                                                                                                                          cGGGTTTCTGTCGACGAAGCCAGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGc                            <  2:2888562/150‑1 (MQ=255)
                                                                                                                            ggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCt                          <  1:1139147/150‑1 (MQ=255)
                                                                                                                            ggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCt                          >  2:1401674/1‑150 (MQ=255)
                                                                                                                               ttCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCa                       <  2:2018161/150‑1 (MQ=255)
                                                                                                                                    tCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTAcc                  <  2:2923196/150‑1 (MQ=255)
                                                                                                                                      gacgaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGc                <  2:3191151/150‑1 (MQ=255)
                                                                                                                                       acgaAGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGcc               >  2:904061/1‑150 (MQ=255)
                                                                                                                                         gaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtg             >  1:3522015/1‑150 (MQ=255)
                                                                                                                                         gaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtg             >  2:3347074/1‑150 (MQ=255)
                                                                                                                                            gCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGTTGTTGGGTGGCCTGCAGTACCGCCTgtgg          >  1:848127/1‑150 (MQ=255)
                                                                                                                                            gCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGTTGTTGGGTGGCCTGCAGTACCGCCTgtgg          >  1:2842109/1‑150 (MQ=255)
                                                                                                                                            gCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgg          >  2:530728/1‑150 (MQ=255)
                                                                                                                                             ccTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtggg         >  2:2760674/1‑150 (MQ=255)
                                                                                                                                             ccTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtggg         >  1:18620/1‑150 (MQ=255)
                                                                                                                                             ccTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtggg         >  2:840315/1‑150 (MQ=255)
                                                                                                                                              cTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggt        <  1:3710774/150‑1 (MQ=255)
                                                                                                                                               tGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACATGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgggggtg       >  2:759193/1‑150 (MQ=255)
                                                                                                                                                ggCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg      <  1:4095047/150‑1 (MQ=255)
                                                                                                                                                  ccTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTgggtgggt    >  2:3961898/1‑150 (MQ=255)
                                                                                                                                                    tCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACAGCCTGTGGGTGGGTaa  >  2:1985377/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGTAA  >  NC_002947/1810229‑1810526

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