Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,742,811 C→G 26.3% intergenic (+80/+672) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,742,8110CG26.3% 164.7 / 74.4 103intergenic (+80/+672)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (44/32);  new base G (14/13);  total (58/45)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

Only 100 of 110 total aligned reads displayed.

GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCT  >  NC_002947/4742663‑4742957
                                                                                                                                                    |                                                                                                                                                   
 gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTCGCGGCCTTGccc                                                                                                                                                   >  2:14227/1‑150 (MQ=255)
 gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTAGCGGCCGTTccc                                                                                                                                                   >  2:14218/1‑150 (MQ=255)
 gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGAGTAGCGGCCGTGccc                                                                                                                                                   >  2:14212/1‑150 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGACAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCCTTTTTTATTGGGCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  2:3529864/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  1:999681/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  2:3548264/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  2:2746063/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  2:1430733/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2536980/150‑1 (MQ=255)
      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2962530/150‑1 (MQ=255)
       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGACAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                                                             <  2:3029451/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAACGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAtga                                                                                                                                            <  2:132235/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:424062/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:4957254/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:231997/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:3737901/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:1661378/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:4522963/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  2:3411301/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  2:2378896/150‑1 (MQ=255)
        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:4710297/150‑1 (MQ=255)
         cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                           <  1:1600835/150‑1 (MQ=255)
           ccAGCTGTACACCGGCGAGCAGACGCGCGACATCGTTGCCCTGAAAGACCGCAAGTCAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                         <  2:5121366/150‑1 (MQ=255)
             aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                                                       <  2:2454484/150‑1 (MQ=255)
               cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGc                                                                                                                                     >  2:4630259/1‑150 (MQ=255)
                      ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACATCGCGGCCTTGCCCGCCAAGAGGCCGGCCCAGaca                                                                                                                              >  1:4039167/1‑150 (MQ=255)
                             gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                                                                     <  1:1768806/120‑1 (MQ=255)
                             gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                                                                     >  2:1768806/1‑120 (MQ=255)
                             gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCAAAGACACCACaaa                                                                                                                       >  1:5114630/1‑150 (MQ=255)
                               aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                     >  2:3250467/1‑150 (MQ=255)
                                 aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGTACAGACACCACAAAACCt                                                                                                                   <  2:4579727/150‑1 (MQ=255)
                                   gcgcGACATCGTTGCCCTGAACGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTATTCGCGGCCGTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTcc                                                                                                                 >  2:3021437/1‑149 (MQ=255)
                                   gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCGGCCGTGGCCGCGGAGGGGCCGGCGCAGCAACCACACAAcccac                                                                                                                 >  1:1811244/1‑147 (MQ=255)
                                             gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCGGCCTTTCCCGCCCAGAGGCCGGCACAGACACCACAAAACAACAGtttttatac                                                                                                      >  1:2636533/1‑146 (MQ=255)
                                             gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCAGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGttttttgtc                                                                                                      >  2:2732394/1‑146 (MQ=255)
                                              ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            >  1:2830228/1‑112 (MQ=255)
                                              ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  2:2830228/112‑1 (MQ=255)
                                                  cgtGAATGACCGCATGTAAGCTGCTGTAGCCGAACGCAAGATGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                                                                 <  2:2955069/148‑1 (MQ=255)
                                                    tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCTCTTAGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACACTTTTTCAccgccccca                                                                                               >  1:2427192/1‑147 (MQ=255)
                                                    tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccgca                                                                                               >  2:375967/1‑149 (MQ=255)
                                                     gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCCAGGCCGCGAaga                                                                                                                                            >  1:4973472/1‑105 (MQ=38)
                                                     gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:4973472/105‑1 (MQ=38)
                                                     gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                                                                              >  2:3528981/1‑150 (MQ=255)
                                                      aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTcc                                                                                             >  2:4121512/1‑149 (MQ=255)
                                                           ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAg                                                                                        >  2:1130867/1‑150 (MQ=255)
                                                                     gCTGCAGTAGCCGAACGCAAGACGGCGGCCCTCGGCCAGCATTTTTTATTTGTCCTGCGCCGGCCTCTCCGCGGCCTTGCCCCCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                              <  2:4857569/150‑1 (MQ=255)
                                                                     gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                              <  2:2038184/150‑1 (MQ=255)
                                                                                cGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTAGCGGGCGTGCCCGCGCACAGGGCCGCCCAGACACCACCAAAACACAGTTTTTGTCCGCCCCCGCACGCAGCCCAGTACGGGGATTGAACgg                                                                   >  1:718759/1‑150 (MQ=255)
                                                                                 taaCGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                                  <  2:4336313/149‑1 (MQ=255)
                                                                                        cgccggcggcCCCGGGGCAGCCTTTTTTTTTTGGCGGGGGCGGGCCTCTGCGCGGCCGGGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                           <  1:4052856/146‑1 (MQ=255)
                                                                                        agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCGGTGCCGGCCTCTTCGGGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                           <  1:4228777/150‑1 (MQ=255)
                                                                                        acccggcggcCCCCGGCCAGCTTTTTTTATTTGCCCGGGGCCGGCCCCTGCGCGGCCTTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                           <  1:4052878/146‑1 (MQ=255)
                                                                                          aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACCCCGCCCGTTCAGGGTATTGTACGGCGCATccacc                                                         >  1:347068/1‑150 (MQ=255)
                                                                                           agGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            >  2:111800/1‑67 (MQ=255)
                                                                                           agGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:111800/67‑1 (MQ=255)
                                                                                               tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGTACGGCGCCTCCACCACGaa                                                    >  2:1492363/1‑150 (MQ=255)
                                                                                                     cGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                              <  2:4929254/150‑1 (MQ=255)
                                                                                                           gCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACAACACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCCTCCACCACGAACTTGTTGGcccc                                        >  2:2751171/1‑148 (MQ=255)
                                                                                                              ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACACACACCACAAAACCTCACTTTTTCTCCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCAACCACCACGAACGTGTTGGCCACCAt                                     >  2:1480513/1‑150 (MQ=255)
                                                                                                              ttttttATTTGTCCTGTGCCGGCCTCGTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATACGCCCGCGCAAGCAGCCCGTTCAGGGGATTGAACCGCGCATACACCACGAACTTTTTGGccaccag                                     >  2:148571/1‑149 (MQ=255)
                                                                                                              ttttttATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACATGTTGGccaccaa                                     >  2:1573103/1‑149 (MQ=255)
                                                                                                                    aTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                               <  1:2915555/150‑1 (MQ=255)
                                                                                                                    aTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                               >  2:4238772/1‑150 (MQ=255)
                                                                                                                       ggTCCTGGGCCGGCCTCTACGCGGCCTGGCCCGCCAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCAc                            <  1:4121512/149‑1 (MQ=255)
                                                                                                                          ccTGTGCCGGACTCTTCGCGGCCTTGCCCGGGTTGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCTCCGGCACGCt                         <  2:1856613/150‑1 (MQ=255)
                                                                                                                            tgtgCCGGCCTCTTCGCGGCCGTGCCCGCGCAGAGGCCGGCACCGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCATCCCCCACGAACTTGTTGGCCACCCACGCCGGCACGCTac                       >  2:4826392/1‑150 (MQ=255)
                                                                                                                            gggcccGCCCTCTCCGGGGCCTGCCCCGCCACGCGCCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                       <  1:14212/146‑1 (MQ=255)
                                                                                                                                ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTATCCCCCCGCTCACGCAGCCCTTTCAGGGTAATGAAAGGCGCATCCACCACGAACATGTTGGCCACCAAAGCCGGCACGCAaccacc                   >  2:4833365/1‑150 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCAAGGCCGCGAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACATGTTGGCCACCATCGCCGGCACGCTACCACCt                  >  1:3833058/1‑150 (MQ=39)
                                                                                                                                  ggCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg                 >  2:1585289/1‑150 (MQ=255)
                                                                                                                                  ggCCTCTGCGCGGCGTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg                 <  1:4869228/150‑1 (MQ=255)
                                                                                                                                     ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt              >  2:2279864/1‑150 (MQ=255)
                                                                                                                                     ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGGATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt              >  2:4800895/1‑150 (MQ=39)
                                                                                                                                     ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt              >  2:787782/1‑150 (MQ=255)
                                                                                                                                     ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt              >  1:1880942/1‑150 (MQ=255)
                                                                                                                                     ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt              >  1:1881059/1‑150 (MQ=255)
                                                                                                                                      tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             >  2:712574/1‑150 (MQ=39)
                                                                                                                                      tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             >  1:4200876/1‑150 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCGGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             <  2:3471245/150‑1 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                >  2:790249/1‑83 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                <  1:790249/83‑1 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             >  2:2923779/1‑150 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             >  1:2462245/1‑150 (MQ=255)
                                                                                                                                      tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTGTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc             >  2:3567617/1‑150 (MQ=255)
                                                                                                                                       ctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg            >  1:1683863/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  1:2108332/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                >  2:4870392/1‑81 (MQ=37)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                <  1:4870392/81‑1 (MQ=37)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  1:4510248/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  2:649366/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGACACGCTACCACCTGGCTCgg           >  1:4779025/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  1:3415042/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  2:2890580/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  2:3498694/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg           >  1:3415080/1‑150 (MQ=255)
                                                                                                                                         tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt          >  1:2650180/1‑150 (MQ=255)
                                                                                                                                           gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt        >  2:158388/1‑150 (MQ=255)
                                                                                                                                              gcggtccccGCGACGCCGCCGGCACAGACACCACAAAACCTCAGATTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGTTATTGAACGGCGCATACACCACGAACTTGTTTGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGc      <  1:2732394/144‑1 (MQ=255)
                                                                                                                                                 ttGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCt  >  2:5077933/1‑150 (MQ=255)
                                                                                                                                                 ttGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGGGCACCt  >  2:3010025/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCT  >  NC_002947/4742663‑4742957

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: