Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 892,459 G→C 10.9% D165H (GAC→CAC)  rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,4590GC10.9% 127.1 / 5.5 42D165H (GAC→CAC) rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base C (2/3);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGCCATTACCTTTTATGGGACGGCGGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGCCGTTGGTGATGGAGGTGAATGCGT  >  NC_000913/892166‑892756
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
ccGCCATTACCTTTTATGGGACGGCGGCACTGCGCCAGTTCGCGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCCAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTTGCGGAGACGCGTCAGGAGGCGGAAAGCGTGATTGACGCttt                                                                                                                                                                                                                                                                                                    >  1:212912/1‑301 (MQ=255)
      ttACCTTTTATGGGACGGCGGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCCGCGGAAAGCGTCAATCACGCTTTCcgccg                                                                                                                                                                                                                                                                                              <  2:411643/301‑3 (MQ=255)
         tcTTTTATGGGACGGCGGTACTTCGCCAGTTCGGGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcgccgcc                                                                                                                                                                                                                                                                                            <  2:328336/300‑5 (MQ=255)
          cTTTTATGGGACGGCGGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACCTGGTCGGTGGTGCGCCGCTGGTGGGCAAGTTGGTTGAAGGCACGCCGGGCATTGGCGTCGTGATGGCGGAGACGCGTCAGGCGGCGGCACACGTCCATGACGCTTTCcgcaggccg                                                                                                                                                                                                                                                                                          >  1:599705/1‑295 (MQ=255)
                                cgcCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCACGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAACACGGGATTCACGGTTTCCCCCGTCTACACGCGCATATTACGGTgacgg                                                                                                                                                                                                                                                                    >  1:452186/1‑297 (MQ=255)
                                          ggttGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGTCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                                                                                                                                                                                                                                                          <  2:35867/298‑1 (MQ=255)
                                            aTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAag                                                                                                                                                                                                                                                        <  2:486307/301‑1 (MQ=255)
                                                 ggggAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCCCGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCACGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGCAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCCTATTCTGGTGCCGGACTATACCAAAGAGgcc                                                                                                                                                                                                                                                   >  1:310909/1‑300 (MQ=255)
                                                       ttATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAggg                                                                                                                                                                                                                                             <  2:452750/300‑1 (MQ=255)
                                                                 aaTGAGTCGGTTGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATcc                                                                                                                                                                                                                                   <  2:311313/301‑1 (MQ=255)
                                                                 aaTGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATcc                                                                                                                                                                                                                                   <  2:539203/301‑1 (MQ=255)
                                                                       tcggtcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTc                                                                                                                                                                                                                             <  2:600466/301‑1 (MQ=255)
                                                                              ccATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGTGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgt                                                                                                                                                                                                                      >  2:461685/1‑301 (MQ=255)
                                                                                    cccGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGa                                                                                                                                                                                                                <  2:162030/301‑1 (MQ=255)
                                                                                                 cAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATCTTCTGGTGCAGGAATCTATCCACGAGGCGCAAGGGGGCGATATCCGCTGTCTTGTAGTTGGCGATGAAGTGGTCGCt                                                                                                                                                                                                   >  1:630582/1‑301 (MQ=255)
                                                                                                        tgTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTg                                                                                                                                                                                            <  2:155275/300‑1 (MQ=255)
                                                                                                                 cAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGACGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGCAGCTTTCCGCGGTCTGACCGCGCATATTCTGGTGCAGGAATCTCTCAAAGAGGCGCACGGGTGCGATATCCGCCTTCAGGTTGTTTGCGATTAAGTGGTcccacgcattgacccgcggg                                                                                                                                                                                   >  1:591530/1‑282 (MQ=255)
                                                                                                                         ggCGCGTGTGGGTATCGATCTGCCTGTCACGGGCATTGTGCATTCGCCGGATGATACCAGCGATTTAATGGACATGGTCGGTGGTGGGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAagag                                                                                                                                                                           <  1:461121/301‑1 (MQ=255)
                                                                                                                                ggcgtgcttctgttgttggttaCGGGTATTGGGCATTTGCCGGATGATACCAGCGTGTTATGGGGCAGGGTCGGTGGTGCGCGTCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGGGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATGCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAtt                                                                                                                                                                    <  1:143964/280‑1 (MQ=255)
                                                                                                                                 aGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGCAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGCATATATCAAAGAGGCGCAAGGGTGCGCTATCACCTGTCTGGTTGTTGGCGATGAAGGAGTCGCTGCGATTGAACGGCGGGCGAAAGCAGGCgcatt                                                                                                                                                                   >  1:493326/1‑297 (MQ=255)
                                                                                                                                   ggCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTc                                                                                                                                                                 >  2:493255/1‑301 (MQ=255)
                                                                                                                                     cATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGt                                                                                                                                                               <  2:493972/301‑1 (MQ=255)
                                                                                                                                       tCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc                                                                                                                                                             <  2:631459/301‑1 (MQ=255)
                                                                                                                                             tgcctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCAAGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGATTTCCGCGGTCTGAACGCGCAACTTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGCGGCAGTGGCCACTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAtat                                                                                                                                                       >  2:619350/1‑299 (MQ=255)
                                                                                                                                                ctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAttagatcggaagagcaca                                                                                                                                                    >  2:144090/1‑285 (MQ=255)
                                                                                                                                                  ctggagttcagacgtgtgctcttccgatctACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGAGTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATc                                                                                                                                                  <  2:137865/272‑1 (MQ=255)
                                                                                                                                                                     cGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGt                                                                                                                               <  2:213155/301‑1 (MQ=255)
                                                                                                                                                                               tACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTAta                                                                                                                     >  2:24323/1‑300 (MQ=255)
                                                                                                                                                                                aCCAGTGATTTAATCGACATGGTCGGGGGTGCGCTGCTGGTGGTCAAGTGGGTTGAAGGCACGAAGGGAATTGGCGGCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGgcgacgcttgccgcgccgccacgatgc                                                                                                                    <  2:329239/301‑21 (MQ=255)
                                                                                                                                                                                         ttgtgtgggaTGGTGGTTGGGGCGCCGCTGGTGGTCAAGTTGGTTTAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGATACGCGTCTGGCGGCGGAAAGCGTGATTGACGCTTGTCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATAGATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCTCTGCGATTGAACGGCGGGCGAAAGAGTGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAgg                                                                                                           <  1:24272/292‑1 (MQ=255)
                                                                                                                                                                                                             tgCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATaaaa                                                                                       <  2:241190/301‑1 (MQ=255)
                                                                                                                                                                                                              gcgcCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAg                                                                                      <  2:592291/301‑1 (MQ=255)
                                                                                                                                                                                                                gccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGCATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCCACTTGCATCGTGGCCGCGCGGCCAGCGTCCCCCGTAACACCACACAAGAGCggcaacccgccaacaaaacc                                                                                    >  1:508605/1‑282 (MQ=255)
                                                                                                                                                                                                                      ggtggtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCCAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGCTTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGCCCGCCCGGCAACCCCCGCCAGTACCCCCCCACAGGAgaccgaaaaccacaaaaacaccccccgc                                                                              >  1:634894/1‑274 (MQ=255)
                                                                                                                                                                                                                                     tgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGACAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGAATCCAACCACAGGAGCGTGAAAACCAGAAAAAAACACCacgaacgaagccacaggac                                                               >  1:75140/1‑290 (MQ=255)
                                                                                                                                                                                                                                     cgccagcacgacgccaattccctgcgtGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGCTTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCCGTATCACCCCCCAGGCGCGTGAAATCGCGATAAAAACCGCCGCAACGATGGCCCTGGag                                                               >  1:187901/25‑300 (MQ=255)
                                                                                                                                                                                                                                             cGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTgg                                                       <  2:75230/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gaCGCGTCAGGCGTCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCggg                           <  2:635771/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cgcgttCAGACCGCGGAAAGCGTCAATCACGCTTTCCGCCGCCTGACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGcc                         >  2:347894/6‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           aaaGCGTCAATCAGGGTTTCCGCCGCTTGACGGGCATTTTTTGGTGCCGGAATATATCAAAGGGGCGCAATGGTGCGATATCCGCTGTCTGGTTGTGGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAGGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGCCGTTGGTGATGGAGGTg         <  1:347431/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                gTGATTGTCGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGCCGTTGGTGATGGAGGTGAATg     <  2:509277/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   aTTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGCCGTTGGTGATGGAGGTGAATGCGt  >  2:193114/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
CCGCCATTACCTTTTATGGGACGGCGGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTGCATCGTGGCGGCGCGGCAAGCGTCGCCAGTATCACACCACAGGAGCGTGAAATCGCGATAAAAGCCGCGCGAACGATGGCGCTGGACGTTGCTGGTGTGGATATTCTGCGTGCTAATCGCGGGCCGTTGGTGATGGAGGTGAATGCGT  >  NC_000913/892166‑892756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: