Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,627 A→G S111S (TCT→TCC)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,6270AG98.1% 138.1 / ‑6.5 52S111S (TCT→TCC) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/22);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

ATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTT‑CTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTG  >  NC_000913/7382‑7867
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
atcgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGGCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCAGCACTTT‑CTTTATAAAGCTGTGCGAGGGAACa                                                                                                                                                                                                                                                <  1:76969/248‑1 (MQ=255)
  cgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACAt                                                                                                                                                                                                                                               <  1:563648/249‑1 (MQ=255)
   gCCAGCAGAGCAAAGACAGCGAGAATAATACCTATCACCAGCGGGGGAAAATCAACAGAAAAACTCAGGGCTCGCGCAACGGCGATCGCCTGAATTCCGCTGAAAATTATGACGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATACCCAGCCCGCGCGCCATATACCATGCCGTTCCGACACGCAACTGCCCATTGGCATACCGAT‑TTTTATAAAGCTGTGCGAGGGAACATTc                                                                                                                                                                                                                                             >  2:676572/1‑251 (MQ=255)
     cagcagAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCgg                                                                                                                                                                                                                                           >  1:274353/1‑251 (MQ=255)
           agcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCaaa                                                                                                                                                                                                                                       >  1:283527/4‑249 (MQ=255)
             caaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACg                                                                                                                                                                                                                                     >  1:444787/2‑249 (MQ=255)
              aaaGACAGCGAGAATAATAACTGTCAACAGCGGGGGAAAGTCAAAAGAAAAACTCAGGGATCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCAATTGGCATCACGTT‑ATTTATAAAGCTGTGCGAGGGAACCTTCGGCAAACGAgg                                                                                                                                                                                                                                  >  1:588138/1‑251 (MQ=255)
                    aGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGcc                                                                                                                                                                                                                             >  1:675250/1‑250 (MQ=255)
                        agaATAATACCTGTCAACAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGAGCAACGGCGGTCGCCTGAATTCCGCTGAAAATTATGCAGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCAGCCACGACACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGCACATTCGGAAAACGAGGTCGcactgcc                                                                                                                                                                                                                        >  2:468907/1‑245 (MQ=255)
                           ataataCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCACGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGCGGTCGCCATTCCGa                                                                                                                                                                                                                      >  2:330703/1‑250 (MQ=255)
                             aataccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATa                                                                                                                                                                                                                    >  1:617306/5‑250 (MQ=255)
                                accTGTCACCAGCGGGGGAAAATCAAAAGGAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAGTTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGATGAGCGAGGGAACATTCGGCAATCGAGGCCGCCATGCCGATAAAcg                                                                                                                                                                                                                >  2:46497/2‑251 (MQ=255)
                                   ggTCACCAGCGGGGGAAAAGCAAAAGAAAAACGCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCgg                                                                                                                                                                                                             <  1:86498/250‑1 (MQ=255)
                                     tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGc                                                                                                                                                                                                            >  1:140285/1‑250 (MQ=255)
                                     tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGGCGCCCTGCCGATAAACGCGGc                                                                                                                                                                                                            >  2:508498/1‑250 (MQ=255)
                                             gggggAACATCAAACGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCCCGCAACTGCCCCTGGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGACAAACGAGGACGCCCTGCCGATAAACGCGGCCACCcacac                                                                                                                                                                                                   >  2:154775/1‑250 (MQ=255)
                                                  aaactCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCCCGCAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGCAGCCCTGCAGATAAACGCGTAAACCCACATcc                                                                                                                                                                                                 >  2:195973/5‑248 (MQ=255)
                                                       cAAAAGAAAAACTCAGGGCACGCGAAACGGCGTTCGCTTGAAGTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAAAAGAACCCCCATCCAGCGCACCCCCAGCCCGCGCGCCTTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGAc                                                                                                                                                                                          <  2:58691/250‑1 (MQ=255)
                                                        aaaaGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGCGCAAAAGGACGGCGAACAGAACGCCCCTCCAGCGCATCACCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGACACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGTGGTCGCCATGCCGCTAAACGCAGCAACCCACATCCAGAAGACgg                                                                                                                                                                                        >  1:320716/1‑251 (MQ=255)
                                                                             cccaccggcgtccgcctaaatcccGCTAAAAATTAGGCCGTGGGCGATGCGCAAAAAGCCGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGGGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTaa                                                                                                                                                                   <  2:672594/229‑1 (MQ=255)
                                                                               cAACGGCGTTCGCCTGAATTCCGATGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAATCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGGGGACGCCATGCAGTTAAACGCGGAAACCCACATCCAGACGACTGCCACAGGTCCACccgcggcaaat                                                                                                                                                                 >  2:657629/1‑246 (MQ=255)
                                                                                aaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCACCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg                                                                                                                                                                <  1:334784/251‑3 (MQ=255)
                                                                                aaCGGCGTTCGCCTGAATCCCGCGGAAAATTATGCCGTAGGCGATGACCAAAAAGACGGCGAACAGAACCCCCATCCAGCGCTCCCCCAGCCCGCGCGCCATATCCCAGCCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg                                                                                                                                                                <  1:263093/251‑3 (MQ=255)
                                                                                      gatCGCCTGAATTCCGCTGAAAATTATGCCGAAGGCGATGCGCAAAAAGACGGCGAACAGCACGCCCAACCGGCGCATCCCCAGCCCGCGCGCCATATACCATGACGGGCCGCCACGAAACTTCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTAGGCAACAGAAGTAGCCCTGCCGAACAACGCGACAACACACATACAGAAGCAGGCTCCACATCCACcgccgggaaaatacagc                                                                                                                                                           >  2:376219/3‑239 (MQ=255)
                                                                                           cctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTTCCCGCGC‑AAcc                                                                                                                                                    >  1:382952/3‑250 (MQ=255)
                                                                                           cctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAAtttcccacgcaacc                                                                                                                                                     >  1:632394/3‑238 (MQ=255)
                                                                                              gaattCCGCCGAAAATTATTCCGTCGGCGATGATCAAAAAGCCGGCGAACAGAACCCCCATCCAGCTCATCCCCAGCCCGCGCGCCATATCCCAAGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCTAGCTC‑AACGCCg                                                                                                                                                 <  1:166098/247‑1 (MQ=255)
                                                                                                           aaTTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGGTGCCGc                                                                                                                                    >  1:67150/1‑251 (MQ=255)
                                                                                                            aTTATGCCGTAGGCAATGAGAAAAAAGCCGGCGAACAGACCGCCCAGCCAGCGCATCCCCAACCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCAGTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGACGCCCTGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCt                                                                                                                                   <  1:46497/251‑1 (MQ=255)
                                                                                                                         cGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGcgcc                                                                                                                      <  1:154775/251‑1 (MQ=255)
                                                                                                                          gATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGACC‑AACGCCGGCCAGGGTGCCCCCACCCACGCGcgccc                                                                                                                     >  1:284770/1‑250 (MQ=255)
                                                                                                                                                     ccccctCCAGTGCACCCCCAGCCCGCGCGCCATATCCCACTCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAACCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATg                                                                                         <  1:565219/246‑1 (MQ=255)
                                                                                                                                                          ccccgccaacccccgccccGCGCCCCATATCCCATCACGGTCCGCCAAGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGGCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAGTGCCAGCGC‑AACGCCGGCCAGGGTGCCGCGACCCACGCGCGTCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑tt                                                                                   <  2:544714/235‑1 (MQ=255)
                                                                                                                                                           ccAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCGCGC‑AACGCCCGCCCGGGTGCCGCTACCCACGCGCGCCGCCATGCTGGGACA‑CAATGaagatcggaagcgcc                                                                                   >  1:480280/1‑237 (MQ=255)
                                                                                                                                                               caccacccccccccGCGCCCCATATCCCATGCGGGTCCGCAACAAAACTGCCCATTGGCACCACGTG‑CTTTATAAAGCTGTGGGAGGGACCATCCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACAGCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAAcc                                                                              <  2:629151/240‑1 (MQ=255)
                                                                                                                                                                gcATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCACGC‑AACGCCGGCCCGGGTGCCGCTAACCACGCGCGCCGCCAGACTGGGACA‑CAAAGACTGAAATGAGG‑TTaaaacc                                                                             >  1:20508/1‑248 (MQ=255)
                                                                                                                                                                       cAGCCCGCGCGCCATATCCCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGTTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCt                                                                      <  1:571211/251‑1 (MQ=255)
                                                                                                                                                                         ccccgcgccccaaaacacATGCCGGTCCACCACAAAGCTGCCCGTTGGCATCA‑GTTGCTTTATAAGGCGGTGCGGGGGAACATCCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCACCCCACATCCAGAAGACGGCACCAGGTCAACCGGCGGTAATTGCCAGCGT‑AACGCCGGCCAGGTGGCCGCTCCCCACGCGCGCCGAAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCtgt                                                                    <  2:107214/234‑1 (MQ=255)
                                                                                                                                                                                  ccATATCCCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATg                                                           <  1:603899/251‑1 (MQ=255)
                                                                                                                                                                                   cATAGACCATGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATGCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGAc                                                                                                               <  1:86826/200‑1 (MQ=255)
                                                                                                                                                                                      ataCCATGCCGTCCCCCCACGAACCTGCCCATGGGCATCACGTT‑CTTAATAAACCGGTGCGAGGTACCATCCGCCAAACGAGGTCGCCTTGCCGATAAACGCGGCAACCCACACCCAGAATACGGCTCCAGGTCCACCGGCGGTAATGGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCCCGCGCGCCGTAAGACTTGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTAt                                                       <  1:340879/251‑1 (MQ=255)
                                                                                                                                                                                          cATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAAATTCCAGCGCCAAC‑CCGGACAGGGTGCCGCCACCCCCACGGCCCGCAAGGCTGGGCAA‑CACTGAATGAAATGAGG‑GTAAACCCCCTGGCTGTGGATGAATGCTCttttt                                                   >  2:415830/1‑251 (MQ=255)
                                                                                                                                                                                               cGGTCCGCCACGAACCTGCCCATTGGCATCCCTTC‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGAc                                              <  1:144357/251‑1 (MQ=255)
                                                                                                                                                                                               cGGTCCGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTa                                                                                  <  1:668087/215‑1 (MQ=255)
                                                                                                                                                                                                    cGCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAATTGCCAGCGC‑AAAGCCGGCCAGGATGCCGCTACCCACGCGCGCCGCCAGACTGGTACA‑CCATGACTGGAATGAGG‑GTAAACCGCCTGGCTTTGGGTGTATGCTATTTTTAAggcctttg                                         >  1:480929/1‑244 (MQ=255)
                                                                                                                                                                                                     gCCACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCTCGC‑AACGCCGGCCAGGGTGCCGCCACCCCCCCGCGCCGCATGACTGGGACA‑CAATGGCTGAACAGAGG‑TTAAACCGCACGGCTTGGGATGAATGATCTTTTTAAAACtttttc                                        >  2:526969/1‑249 (MQ=255)
                                                                                                                                                                                                       cACGAAACTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGC‑AACGCCGGCCCGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGGACA‑CCATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGGTGAATGCTATTTTTAAGACTTTTGccc                                      >  1:364508/1‑250 (MQ=255)
                                                                                                                                                                                                          caaaaccgcccattggcatcccttacTTTATAAAGCTGTGCGAGGGAACATGCGGCAAATGGGGTCGCCAGGCCGAAAAACGGGGCAACCCACACCCAGAAGACGGCTCCTGGTCCACGGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAGGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCaaa                                    <  2:675250/226‑1 (MQ=255)
                                                                                                                                                                                                             aCTGCCCATTGGCATCACGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAATTGCCAGCGC‑AAAGCCGGCCAGGGTGCCGCTACCCCCGCGCGCCGCAAGGCTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCGGGGTGTGGGAGAATGCTATTTTTAAACCTTTTGCCAAACTgg                                >  1:227743/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGTT‑CTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCGCCC‑AACGCCGGCCAGGTTGCCGCTACCCCCGCTCGCAGCAAGACTGG‑ACAACAATGCTTGAAATTGGG‑TTTAACCGCCTCGATGTGGGTGAATGCTCATTTTAACACTTTTTCCAAACTGGCggggctagaggaaaa                >  2:7721/1‑237 (MQ=255)
                                                                                                                                                                                                                                     ttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTTATAGCCAGCCC‑CACCCCGGCCAGGGTGCCGCCACCCACCCGCGCCGCCAGAACGGGACA‑CAATAGCTGCAATTAGGGTTAAACC‑CCCGGCTGGGGAGGAAGGCGATTTTTAAAAAATTTGACCAACTGGCgggggggaaggaaacacacac          >  2:120823/1‑231 (MQ=255)
                                                                                                                                                                                                                                            gctgtgCGAGGGAACATGCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCCCCGGCGGTAATTGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTg  <  1:40028/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
ATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTT‑CTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGC‑AACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACA‑CAATGACTGAAATGAGG‑TTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTG  >  NC_000913/7382‑7867

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑