Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,398 T→G A228A (GCT→GCG nhaR → DNA‑binding transcriptional activator NhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,3980TG98.0% 149.4 / ‑4.0 49A228A (GCT→GCGnhaRDNA‑binding transcriptional activator NhaR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (21/27);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

ATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCA  >  NC_000913/19153‑19626
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
atGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTtgat                                                                                                                                                                                                                                 >  1:97852/1‑251 (MQ=255)
     cagttCTGACTGTCCGATAGACGCTACGCAGCAGGAAGGCCTGTTCTCTGTTAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTCCAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGaaa                                                                                                                                                                                                                             <  1:197721/247‑1 (MQ=255)
     cATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGGAGAACGGCGGCTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTATAAAACCTCGGCGAGTTTGATGTAGCCGCGTTGATGaa                                                                                                                                                                                                                              >  2:228781/1‑249 (MQ=255)
             aCTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGACTTGGCGAATGTGGCGTGAGTTTCTGGTGCACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGCCGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGTGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGt                                                                                                                                                                                                                    >  1:30778/1‑251 (MQ=255)
               ggttcgttaggatttaaccaccagGAAGGCCTGTGCGCCGTGAGAATGGGCGAATGTTGCGTGAGTTTCTGGTTTACAAAGCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTGTGATTCCTGGGCTACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGACGATGCCGCGTTGATGAAAGCTTTTGGTGc                                                                                                                                                                                                                  <  2:399162/230‑1 (MQ=255)
                     aTAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGCTGc                                                                                                                                                                                                             >  2:1000/1‑250 (MQ=255)
                       aGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGcaca                                                                                                                                                                                                          <  1:245713/251‑1 (MQ=255)
                              aCGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGc                                                                                                                                                                                                      <  2:97852/248‑1 (MQ=255)
                                            ccgggtctcCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGGGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAAGCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCaa                                                                                                                                                                                     <  2:127579/246‑1 (MQ=255)
                                             cTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGCCGTTCAATGTTAGGGCGCCAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGGTGCACAATGCAATCTTCGTTGCCCCAAc                                                                                                                                                                                    >  2:395795/1‑251 (MQ=255)
                                                    ccGTGAGAATTGGCGAATGGGGCGTTAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTTCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTAAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCTCAACGCTTTAt                                                                                                                                                                             <  2:245521/251‑1 (MQ=255)
                                                        gagaATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCa                                                                                                                                                                          <  1:103196/250‑1 (MQ=255)
                                                        gagaATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCa                                                                                                                                                                          <  1:395795/250‑1 (MQ=255)
                                                         agaATTGGCGAATGTGGGGTGATTTTCTGGTGAACAAAGCCACCACCAGAAAAACCGTGCCCGGCTTGACTGGAAGAAAGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGTTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCaga                                                                                                                                                                        <  2:265392/251‑3 (MQ=255)
                                                            agtGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATTAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGa                                                                                                                                                                     <  1:159956/249‑1 (MQ=255)
                                                             ttGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCAAATGAc                                                                                                                                                                    >  2:144751/1‑251 (MQ=255)
                                                                    aTGTGTCGTGAGTGTCTGGTGTACAAATTCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGTCGTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTGAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTTCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATg                                                                                                                                                             <  2:457313/251‑1 (MQ=255)
                                                                      gtgtcGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGcc                                                                                                                                                           <  1:174950/247‑1 (MQ=255)
                                                                                tttCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGATCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCAAACGCTTTATGCATATGACTTTTATTCCGATAAAACTg                                                                                                                                                 >  2:344642/1‑251 (MQ=255)
                                                                                                 caccaccaGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTcgc                                                                                                                                <  1:195085/251‑1 (MQ=255)
                                                                                                     accaGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCg                                                                                                                            >  1:107463/1‑251 (MQ=255)
                                                                                                       caGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCgag                                                                                                                          <  1:252307/251‑1 (MQ=255)
                                                                                                          aaaaaCCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCTGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAAt                                                                                                                       <  2:106778/251‑1 (MQ=255)
                                                                                                                            gTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCAt                                                                                                     >  1:13394/1‑251 (MQ=255)
                                                                                                                               tGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAACGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATTCCGATAAATCTGTCGTAGAAATTGGTCGCGTCGGGAATGTGTTGGAATCGAACCATGCt                                                                                                  >  2:236620/1‑251 (MQ=255)
                                                                                                                                             cTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCg                                                                                    >  1:204474/1‑251 (MQ=255)
                                                                                                                                                             gACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCcg                                                                    >  1:168413/1‑251 (MQ=255)
                                                                                                                                                              acggtCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATGTTTGCTGAGCGGATGATTCAGCACCcgg                                                                   <  2:295730/247‑1 (MQ=255)
                                                                                                                                                                gTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCcggcg                                                                 >  1:186635/1‑251 (MQ=255)
                                                                                                                                                                              cgcAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATccg                                                   <  1:260111/251‑3 (MQ=255)
                                                                                                                                                                                 aaaTTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACa                                                            >  2:347722/1‑239 (MQ=255)
                                                                                                                                                                                 aaaTTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACa                                                            <  1:347722/239‑1 (MQ=255)
                                                                                                                                                                                          aaCTGGTTTAACACCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGatta                                       <  1:484702/251‑1 (MQ=255)
                                                                                                                                                                                                  tAACCCCCAGGGATTAAACGTAGAATTCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTGCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCTAATCTGCAATACGGATTATTCTgcg                                <  2:30778/250‑1 (MQ=255)
                                                                                                                                                                                                     cTCCCAGGGATTAAAAGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGTTGAGCGGATGATTCAGCACCAGGCGGTACAGCGAATCTGCAATAAGGATTATTCTGCGCttt                            >  1:264875/1‑251 (MQ=255)
                                                                                                                                                                                                        ccAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGCTGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCttttt                          >  2:207084/1‑250 (MQ=255)
                                                                                                                                                                                                        ccAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGAACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACACGGAATCTGCAATACGGGTTATTCTGCGCttttt                          >  2:122051/1‑250 (MQ=255)
                                                                                                                                                                                                         cAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCttttt                          >  2:231869/1‑249 (MQ=255)
                                                                                                                                                                                                                tAAACGTAGAAAGCCTCGGCGAGTGTGATGATGCCGTGTTGATGAAAGCTTTTGGTGCGATGCACGATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTGTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATTCTATTTTTGCTGAGCGGATGATTCATCACCCGGCGGTACAGCGAAGCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGc                  <  2:34833/250‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCggc               <  1:231869/251‑2 (MQ=255)
                                                                                                                                                                                                                    cGTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCggcgc             <  2:80865/251‑4 (MQ=255)
                                                                                                                                                                                                                     gTAGAAATCCTCGGCGAGTTTGATGATGCCGCGTTGATGACAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGAAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACAATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTTCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGcgcg            >  2:468193/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaTCCTCGGCGAGTTTGATGATGCCTCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAAt       <  1:259575/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaTCCTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAAt       <  1:397177/251‑1 (MQ=255)
                                                                                                                                                                                                                             ccTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATa                                              >  1:6406/1‑209 (MQ=255)
                                                                                                                                                                                                                             ccTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATa                                              <  2:6406/209‑1 (MQ=255)
                                                                                                                                                                                                                             ccTCGGCGAGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCg     <  1:1000/250‑1 (MQ=255)
                                                                                                                                                                                                                               tCGGCGCGCTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGCTGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATCTGACTTTTATGCCGATAAAACTGTCGCAGCAGTTGGTCGCGTCGCGAATGTGATGGAATGGCACCATGCTATTTTTGTTGAGCGGATGCTTCCGCACCCGGCGGCAACGCGAACCTGCAATCCGTATTAATCTGCGCtttttcagccagcggcgagttattcgta   >  2:253326/1‑227 (MQ=255)
                                                                                                                                                                                                                               ccGGCGAGGTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCa  <  2:186635/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
ATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGAAAAACCGTTCCCGGCTTGTCTGGAAGAACGGCGACTTTTGATTCCTGGGCGACGTTCAATGTTAGGGCGCAAATTGCTTAACTGGTTTAACTCCCAGGGATTAAACGTAGAAATCCTCGGCGAGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCA  >  NC_000913/19153‑19626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: