Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 65,811 G→C 6.9% intergenic (‑31/+44) 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,8110GC 6.9% 432.5 / 17.8 131intergenic (‑31/+44)polB/araDDNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (62/60);  new base C (6/3);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 133 total aligned reads displayed.

CTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAAT  >  NC_000913/65663‑65957
                                                                                                                                                    |                                                                                                                                                  
ctTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGt                                                                                                                                                   <  2:1457945/150‑1 (MQ=255)
        gCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATATAACCTGCCTGAGCCACGCTGCAATTCCATCAAAAAACCAGGCTTGAGTATAGCCTg                                                                                                                                           >  1:1381104/1‑150 (MQ=255)
         cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTgg                                                                                                                                          <  2:665130/150‑1 (MQ=255)
         cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTgg                                                                                                                                          <  2:1413765/150‑1 (MQ=255)
            gCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                       >  2:1225128/1‑150 (MQ=255)
              gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttttt                                                                                                                                     <  2:905002/150‑3 (MQ=255)
              gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCg                                                                                                                                     >  2:382553/1‑150 (MQ=255)
                aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                   <  1:1732010/150‑1 (MQ=255)
                aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                   >  2:721021/1‑150 (MQ=255)
                aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                   >  1:831151/1‑150 (MQ=255)
                aaCCTGCAACGGCCCGTTGTCCGGCGCCAGCCATAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGttttttg                                                                                                                                   <  1:1106173/150‑5 (MQ=255)
                 aCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGttt                                                                                                                                  <  2:1380313/150‑1 (MQ=255)
                 aCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGGGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGttt                                                                                                                                  <  2:441776/150‑1 (MQ=255)
                  ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                 <  2:128031/150‑1 (MQ=255)
                   cTGCAACGGTCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAg                                                                                                                                                    >  2:1646951/1‑130 (MQ=255)
                   cTGCAACGGTCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAg                                                                                                                                                    <  1:1646951/130‑1 (MQ=255)
                   cTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGa                                                                                                                                >  1:396313/1‑150 (MQ=255)
                    tGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGAt                                                                                                                               <  1:951280/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)
                               gTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGtt                                                                                                                    >  2:389905/1‑150 (MQ=255)
                                 tGTCCGTCGCCAGCCAGAATGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGtttt                                                                                                                  <  2:1496078/150‑1 (MQ=255)
                                  gTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCCGGCGTGAGTATAGCCTGGTGTCGTTTGATTGGCTGTGGTTTta                                                                                                                 >  2:862831/1‑150 (MQ=255)
                                   tCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTcat                                                                                                                >  2:1688662/1‑150 (MQ=255)
                                       tCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAg                                                                                                            <  1:568909/150‑1 (MQ=255)
                                        cGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGt                                                                                                           >  2:1404579/1‑150 (MQ=255)
                                        cGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGt                                                                                                           >  1:1475412/1‑150 (MQ=255)
                                             gccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTa                                                                                                      >  1:1624193/1‑150 (MQ=255)
                                              ccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTAc                                                                                                     <  2:530416/150‑1 (MQ=255)
                                              ccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTAc                                                                                                     <  1:1011619/150‑1 (MQ=255)
                                              ccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTAc                                                                                                     >  1:1526796/1‑150 (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)
                                                           tCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATg                                                                                        <  1:129995/150‑1 (MQ=255)
                                                            cTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGc                                                                                       >  2:762101/1‑150 (MQ=255)
                                                             tGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGcc                                                                                      <  1:1453623/150‑1 (MQ=255)
                                                               tCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCAGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                    >  1:680028/1‑150 (MQ=255)
                                                                  cTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTcgc                                                                                 <  1:1446397/150‑1 (MQ=255)
                                                                       ggggTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCAt                                                                            >  2:583958/1‑150 (MQ=255)
                                                                         ggTGGCCCGCCAGTGTCGGGTTAAGATAATACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAAACAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGTCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGc                                                                          <  2:932988/150‑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)
                                                                                    aGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATa                                                               >  1:242731/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)
                                                                                           ggTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTa                                                        >  2:24736/1‑150 (MQ=255)
                                                                                           ggTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTa                                                        >  1:577928/1‑150 (MQ=255)
                                                                                            gTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTAt                                                       >  1:1033686/1‑150 (MQ=255)
                                                                                            gTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGGTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTAt                                                       <  2:834977/150‑1 (MQ=255)
                                                                                                aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                   >  1:1309699/1‑150 (MQ=255)
                                                                                                aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                   >  1:452386/1‑150 (MQ=255)
                                                                                                aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCGGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                   <  2:1526796/150‑1 (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:1103331/1‑150 (MQ=255)
                                                                                                    aaaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCagcagc                                               >  2:370606/1‑150 (MQ=255)
                                                                                                     aaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCg                                              >  2:594205/1‑150 (MQ=255)
                                                                                                     aaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCg                                              <  1:1706061/150‑1 (MQ=255)
                                                                                                     aaaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCg                                              >  2:498899/1‑150 (MQ=255)
                                                                                                      aaCCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGt                                             >  1:1210795/1‑150 (MQ=255)
                                                                                                      aaCCTGCCTGCGCCACGCTGAAAATCCATAAAAAAACCAGTCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGt                                             <  1:735788/150‑1 (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                              <  2:796763/150‑1 (MQ=255)
                                                                                                                     cGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCgg                              <  1:155154/150‑1 (MQ=255)
                                                                                                                      gCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGt                             >  2:1663894/1‑150 (MQ=255)
                                                                                                                      gCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGt                             <  1:515555/150‑1 (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)
                                                                                                                                   cAAAAAACCAGGCTTGAGTATAGCCTGGTTTGGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCGATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                <  2:795647/150‑1 (MQ=255)
                                                                                                                                   cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                >  1:105670/1‑150 (MQ=255)
                                                                                                                                   aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                >  1:503243/5‑150 (MQ=255)
                                                                                                                                   aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                >  1:215454/5‑150 (MQ=255)
                                                                                                                                   aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                >  1:823557/5‑150 (MQ=255)
                                                                                                                                     cgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTTTCAAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc              >  1:1625610/3‑150 (MQ=255)
                                                                                                                                     cgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc              >  1:1498423/3‑150 (MQ=255)
                                                                                                                                     aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc              <  2:1627311/150‑1 (MQ=255)
                                                                                                                                     aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc              >  2:1012858/1‑150 (MQ=255)
                                                                                                                                          ccAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACg         <  1:941599/150‑1 (MQ=255)
                                                                                                                                           cAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACgg        >  1:722176/1‑150 (MQ=255)
                                                                                                                                            aGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGc       <  2:349648/150‑1 (MQ=255)
                                                                                                                                               cTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGa    <  1:498899/150‑1 (MQ=255)
                                                                                                                                                ttGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGaa   <  2:411152/150‑1 (MQ=255)
                                                                                                                                                 tGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAat  >  2:404787/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAAT  >  NC_000913/65663‑65957

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: