Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 65,813 A→C 6.7% intergenic (‑33/+42) polB ← / ← araD DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,8130AC 6.7% 382.4 / 13.5 130intergenic (‑33/+42)polB/araDDNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (61/59);  new base C (6/3);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 131 total aligned reads displayed.

GCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCC  >  NC_000913/65671‑65962
                                                                                                                                              |                                                                                                                                                     
gCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATATAACCTGCCTGAGCCACGCTGCAATTCCATCAAAAAACCAGGCTTGAGTATAGCCTg                                                                                                                                                >  1:1381104/1‑150 (MQ=255)
 cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTgg                                                                                                                                               <  2:1413765/150‑1 (MQ=255)
 cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTgg                                                                                                                                               <  2:665130/150‑1 (MQ=255)
   aGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGTAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGACACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGtt                                                                                                                                             <  2:482225/150‑1 (MQ=255)
   aGCGTAACCTGCAACGGCCCGTTGTCCGTCGCAAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGtt                                                                                                                                             >  2:1446890/1‑150 (MQ=255)
    gCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                            >  2:1225128/1‑150 (MQ=255)
     cGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTc                                                                                                                                           >  1:796763/1‑150 (MQ=255)
      gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttttt                                                                                                                                          <  2:905002/150‑3 (MQ=255)
      gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCg                                                                                                                                          >  2:382553/1‑150 (MQ=255)
        aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCGCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCACCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        <  2:567289/150‑1 (MQ=255)
        aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        <  1:1732010/150‑1 (MQ=255)
        aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        >  1:831151/1‑150 (MQ=255)
        aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        >  2:721021/1‑150 (MQ=255)
        aaCCTGCAACGGCCCGTTGTCCGGCGCCAGCCATAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGttttttg                                                                                                                                        <  1:1106173/150‑5 (MQ=255)
         aCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGttt                                                                                                                                       <  2:1380313/150‑1 (MQ=255)
          ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      <  2:128031/150‑1 (MQ=255)
           cTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGa                                                                                                                                     >  1:396313/1‑150 (MQ=255)
            tGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGAt                                                                                                                                    <  1:951280/150‑1 (MQ=255)
             gCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGAtt                                                                                                                                   <  1:1599495/150‑1 (MQ=255)
                 cGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCt                                                                                                                               <  1:44659/150‑1 (MQ=255)
                   gCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCtgt                                                                                                                             >  1:1628033/1‑150 (MQ=255)
                     ccGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                           <  2:175750/150‑1 (MQ=255)
                      cGTTGTCCGTCGCCAGCCAGAAGGAGACTTGTGTCCCTTGCGGTGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGt                                                                                                                          <  1:1566702/150‑1 (MQ=255)
                          gTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCCGGCGTGAGTATAGCCTGGTGTCGTTTGATTGGCTGTGGTTTta                                                                                                                      >  2:862831/1‑150 (MQ=255)
                             cGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATAc                                                                                                                   >  1:523541/1‑150 (MQ=255)
                               tCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAg                                                                                                                 <  2:355324/150‑1 (MQ=255)
                               tCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAg                                                                                                                 <  1:568909/150‑1 (MQ=255)
                                cGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGt                                                                                                                >  2:1404579/1‑150 (MQ=255)
                                      ccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTAc                                                                                                          <  1:1011619/150‑1 (MQ=255)
                                       cagAATGAGACTTCTGTCGCTTGCGGGGTGTCCCGCCAGTGTCGGGGTAAGATAAACCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACt                                                                                                         <  2:254913/150‑1 (MQ=255)
                                       cagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACt                                                                                                         >  1:1420185/1‑150 (MQ=255)
                                       cagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCGTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACt                                                                                                         <  2:989873/150‑1 (MQ=255)
                                        agAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTg                                                                                                        >  2:349152/1‑150 (MQ=255)
                                         gcaGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCCGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGc                                                                                                       <  1:1369272/148‑1 (MQ=255)
                                           aGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGccc                                                                                                     <  2:831151/150‑1 (MQ=255)
                                           aGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGccc                                                                                                     >  1:1726263/1‑150 (MQ=255)
                                                   tCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATg                                                                                             <  1:129995/150‑1 (MQ=255)
                                                   tCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATg                                                                                             <  1:1466929/150‑1 (MQ=255)
                                                    cTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGc                                                                                            >  2:762101/1‑150 (MQ=255)
                                                     tGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGcc                                                                                           <  1:1453623/150‑1 (MQ=255)
                                                     tGTCCCTTGCGGGGTGTCCCGCCAGTGACGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGcc                                                                                           <  1:1633424/150‑1 (MQ=255)
                                                       tCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCAGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                         >  1:680028/1‑150 (MQ=255)
                                                          cTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTcgc                                                                                      <  1:1446397/150‑1 (MQ=255)
                                                            tGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCgcgc                                                                                    >  1:488979/1‑150 (MQ=255)
                                                               ggggTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCAt                                                                                 >  2:583958/1‑150 (MQ=255)
                                                                 ggTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGc                                                                               >  2:1660584/1‑150 (MQ=255)
                                                                 ggTGGCCCGCCAGTGTCGGGTTAAGATAATACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAAACAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGTCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGc                                                                               <  2:932988/150‑1 (MQ=255)
                                                                   tgtCCCGCCAGTGTTGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATg                                                                                >  2:295015/1‑147 (MQ=255)
                                                                   tgtCCCGCCAGTGTTGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATg                                                                                <  1:295015/147‑1 (MQ=255)
                                                                   taacccGACACTGGCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCtt                                                                             >  2:651426/4‑150 (MQ=255)
                                                                        cGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTAGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCa                                                                        >  2:316106/1‑150 (MQ=255)
                                                                         gCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAg                                                                       >  1:1413273/1‑150 (MQ=255)
                                                                           cAGTGTCGGGGTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGAt                                                                     <  2:1517511/150‑1 (MQ=255)
                                                                             gtgtCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAg                                                                   >  1:1523474/1‑150 (MQ=255)
                                                                               gtCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTTTAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAgtg                                                                 >  1:122437/1‑150 (MQ=255)
                                                                                tCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAgtgt                                                                <  1:895699/150‑1 (MQ=255)
                                                                                 cGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGtt                                                               >  2:1508495/1‑150 (MQ=255)
                                                                                  gggTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGttt                                                              >  2:1604423/1‑150 (MQ=255)
                                                                                    gTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTAt                                                            >  1:1033686/1‑150 (MQ=255)
                                                                                    gTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGGTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTAt                                                            <  2:834977/150‑1 (MQ=255)
                                                                                      tAAGATAAATCCTGCCAGCGCCACGCTGAAAATCCATCAAAACACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATcc                                                          <  2:240747/150‑1 (MQ=255)
                                                                                        aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                        >  1:1309699/1‑150 (MQ=255)
                                                                                         gATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagc                                                       >  1:1215693/1‑150 (MQ=255)
                                                                                         gATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagc                                                       >  1:891219/1‑150 (MQ=255)
                                                                                         gATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagc                                                       <  1:1446890/150‑1 (MQ=255)
                                                                                           tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTCATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagcag                                                     >  2:244017/1‑150 (MQ=255)
                                                                                            aaaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCagcagc                                                    >  2:370606/1‑150 (MQ=255)
                                                                                            aaaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCagcagc                                                    >  2:1103331/1‑150 (MQ=255)
                                                                                             aaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCg                                                   >  2:498899/1‑150 (MQ=255)
                                                                                             aaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCg                                                   >  2:594205/1‑150 (MQ=255)
                                                                                              aaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGt                                                  >  1:1210795/1‑150 (MQ=255)
                                                                                                         gcCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCataa                                       >  2:60426/1‑149 (MQ=255)
                                                                                                          ccACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATc                                      >  1:799234/1‑150 (MQ=255)
                                                                                                           cACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATcc                                     <  2:523541/150‑1 (MQ=255)
                                                                                                            aCGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCg                                    <  1:576039/150‑1 (MQ=255)
                                                                                                             cGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCgg                                   <  1:155154/150‑1 (MQ=255)
                                                                                                             cGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCgg                                   <  2:796763/150‑1 (MQ=255)
                                                                                                              gCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGt                                  >  2:1663894/1‑150 (MQ=255)
                                                                                                                  aaaaTCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACt                              <  1:1508495/150‑1 (MQ=255)
                                                                                                                   aaaTCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTg                             <  2:1184330/150‑1 (MQ=255)
                                                                                                                      tCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTgcgg                          <  1:1718748/150‑1 (MQ=255)
                                                                                                                        cATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTAGGCAGATAGTGTTTTTGCAGCAGCGTGTGCTGCATATCCGGTTACTgcggcg                        >  2:162805/1‑150 (MQ=255)
                                                                                                                        cATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTgcggcg                        >  2:1553337/1‑150 (MQ=255)
                                                                                                                         aTCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGgcgc                       >  1:647071/1‑150 (MQ=255)
                                                                                                                           cAAAAAACCAGGCTTGAGTATAGCCTGGTTTGGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCGATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     <  2:795647/150‑1 (MQ=255)
                                                                                                                           cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCtt                     >  2:403712/1‑149 (MQ=255)
                                                                                                                           cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  1:105670/1‑150 (MQ=255)
                                                                                                                           aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  1:215454/5‑150 (MQ=255)
                                                                                                                           aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  1:823557/5‑150 (MQ=255)
                                                                                                                           aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  1:503243/5‑150 (MQ=255)
                                                                                                                             cgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTTTCAAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc                   >  1:1625610/3‑150 (MQ=255)
                                                                                                                             cgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc                   >  1:1498423/3‑150 (MQ=255)
                                                                                                                             aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc                   >  2:1012858/1‑150 (MQ=255)
                                                                                                                                 aCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTAATACAGTCATTACTGCCCGAAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGAc               <  2:1516681/150‑1 (MQ=255)
                                                                                                                                  ccAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACg              <  1:941599/150‑1 (MQ=255)
                                                                                                                                   cAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACgg             >  1:722176/1‑150 (MQ=255)
                                                                                                                                       cTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGa         <  1:498899/150‑1 (MQ=255)
                                                                                                                                        ttGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGaa        <  2:411152/150‑1 (MQ=255)
                                                                                                                                         tGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAat       >  2:404787/1‑150 (MQ=255)
                                                                                                                                              ataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCACGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGATTATccc  >  2:1556962/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCC  >  NC_000913/65671‑65962

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