Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,782,495 C→T 21.9% S77N (AGT→AAT)  ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,782,4950CT21.9% 115.7 / 4.8 60S77N (AGT→AAT) ypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base C (28/16);  new base T (3/11);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCATTCGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATA  >  NC_000913/2782267‑2782706
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   
tttCCATTCGGCCCCCCCCCCTCTACGCGTGAAGCAACGCCCTTAATTTGCGTTCCATTGCTGATCCCCCCTTTATCAACGTCCTGATGGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCAGTCAGGAATTGACTTTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:382716/228‑1 (MQ=255)
tttCCATTCCCCCCCCCTCCCCCTCCGCGTGAAGCCCCGCCCTTAATTGTTGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTTCCTCTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:487335/228‑1 (MQ=255)
         gCACTCACTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTCATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCACCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAC‑TTGTTATGACTTCCCATTT‑CTTTTTA‑TCCACCAGCTTCAATAACAATTCCATTaa                                                                                                                                                                                                            >  2:502943/1‑228 (MQ=255)
                       tACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCAACCG‑TTAATGGTTGTATTTTTAGCCTG‑CCCCAC‑AAAA‑TTGTTAAGACTTCCTATTA‑CTTTAAA‑ACCACCAGGTTCACAAAGATTTCCATAAAAACCCCCGccctaa                                                                                                                                                                                              >  2:313582/1‑225 (MQ=255)
                         cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGCTGAATTATTAGCCTT‑CCCCAC‑AAAA‑TTGTTATTACTTCCAATAT‑CTTGTTA‑TCATCCAGATTCAATAAGACTTCCATTAGATACcccccctttaa                                                                                                                                                                                             >  1:341145/1‑223 (MQ=255)
                         cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGCCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATGCCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTACCCTG‑CCCCAC‑AAAA‑TTGTTATGACTCCCTATAT‑CTTGATT‑TCACCCAGATTCAAAAAAATTTCCATTAAAAAACCCGCCATGaa                                                                                                                                                                                             >  2:350525/1‑228 (MQ=255)
                          ccgTGAAGCCCCGCCCTTAATTGTCGTTGCATTGCTGATCCCCCCTTTATAAACGTCCTTATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAA‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACa                                                                                                                                                                                           <  1:680512/228‑1 (MQ=255)
                               aaGCCCCGCCCTTAATTGTCGTTCCATTGCTGAGACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTGCCTATAT‑CTTGATA‑CCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  1:544526/198‑1 (MQ=255)
                               aaGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAACCGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCAGCCG‑TTAATGTTTGTATTATTAGCCTG‑CCCCAC‑ACAA‑TTGTTATGACTTCCTTTAT‑CTTGATA‑TCCCCCAGTATCAATCAGATCTCCTTTAGAtctgtttcttatatcaccttg                                                                                                                                                                                      >  2:544526/1‑209 (MQ=255)
                                             tttGTCGTTCCATTGCGTATACCCCCTTTAAAAACGTCCTGTTTGCCCCTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCCCCT‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGCCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                          <  1:176438/226‑1 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCACCCAGATTCAATAAGACTTACATTAGCTACCCCGCCATGAACATTCTTCTGGCCATGGTTGataat                                                                                                                                                                    >  2:247549/1‑229 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTAGGACTTCCTATTT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCAGGAAAATTCTCCTGGTCATGGTTgtaaaa                                                                                                                                                                    >  2:606978/1‑224 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTC‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGatcat                                                                                                                                                                    >  1:380817/1‑229 (MQ=255)
                                                     tCCGTTTCTGATCCCCCCTTAATAACCGCCCTGATTCCCACTCTCGTTTGTCGTACCTGTGGAAATACCCCTGTAATGAATTGATTTTCTGCCACCT‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGCCTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                                                                                                                                                <  1:269969/229‑1 (MQ=255)
                                                          ttttttcccccctTTTAAAACCCTCCTTATTGCCCTTTTTCTTTGTTGTCCCTGTGTAAATCCCCCCGTCAGTAATTGAATGTTTTCCACCT‑TTAATGGTTGTTTTATTCGCCT‑‑CCCCACAAAAA‑TTGTTATGTCTTCATATAT‑CTTGATA‑TCCACCAGTTTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                           <  2:302701/221‑1 (MQ=255)
                                                          tttgaatccccccTTACTAAGCTTCCTGATTGCACCCCCCGTTTGTCGTACTTGTGGAAATCCCCCTACCATATATTGACTTTTGGCCACCG‑TTATTGCTTGTATTATTATCTTG‑CCCAAA‑AAAT‑TTGTTATGTCTTCCTATAT‑CTTGATA‑TCCCCCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                           <  1:312215/221‑1 (MQ=255)
                                                               aTACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  1:737188/166‑1 (MQ=255)
                                                               aTACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     >  2:737187/1‑166 (MQ=255)
                                                                 aCCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGAAACCCCTCCATGAACATTCGTCTGGCCATGGTTGACAATATGGTGGTTAttttt                                                                                                                                                     >  2:688808/1‑226 (MQ=255)
                                                                               cttCCTGATGGCCCCTCTCTATTGTCGTCCCTGTGAAAATCCCCCTGTCCTGAATTTACTGTCTGCCACCG‑TTAATGGTTGTATTATTACCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTGGATA‑TCCACCAGATTCAATAAGACTTCCATGAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATc                                                                                                                                      <  1:451821/227‑1 (MQ=255)
                                                                                 tCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCCATAAGACTTCCCTTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTTCTAATATGAGTGTTCTTTGTTGTACCTCGTTAATCCa                                                                                                                                    >  1:738910/1‑229 (MQ=255)
                                                                                 tCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGAATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑ATTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTCACGGCCATGGTTGATAATATGATTGTTATTTGTTGTACCTCGTTCATcc                                                                                                                                     >  1:349085/1‑228 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCTTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGCTTCACTAAGACTTCCCTTAGATACCCCCCCATGAACATTCTGCTGGCCATGATTGCTAATATCAGTGTTTTTTTTTGTACCCCGttaaccccc                                                                                                                                   >  2:215677/1‑221 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     >  1:493867/1‑146 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:493867/146‑1 (MQ=255)
                                                                                                    ttttcgacgtttgtaaaaTCCCCCCTGTATGAATTGATTTTCTGCCACCA‑TTAATGGTTGTATTATTGGCCGG‑CCCAAC‑AAAA‑TTATTATGACTTCCTTTAT‑CTTGATT‑ACCACCAGATTCAATAAGACTTCATTTAGATTCCCCGCCATGATCATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTAc                                                                                                                 <  1:313582/214‑1 (MQ=255)
                                                                                                         gttctggtgGAATTACCGCCGTCATGATGTTCTTTTTTTCCACCT‑TTAAGGTTTGTATTATTACCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                            <  1:606979/223‑1 (MQ=255)
                                                                                                            ctttggtattcccccccGTCATATATTAATTGTCTCCCCCTT‑TAAATGGTTATTTTATTACCCCT‑CCCCAA‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑CCCAACAGTTTAAATAAGACTTCCATTAGATTCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCAt                                                                                                         <  2:124778/217‑1 (MQ=255)
                                                                                                                          cccGCCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGCTCATCCACTTTGTGGTTGCCATCTACAGTCTCAACGTTTACCACACCa                                                                                           >  2:408158/1‑229 (MQ=255)
                                                                                                                            tatcttaaattaatttcttGCCACGG‑TTATTGGTTG‑ATTATTAACCCGCCCCCAA‑AAAA‑TTGTTAGAACTCCCTGTAT‑TCAGATA‑TCCACCAAATTAAATAAGATTTCTATTAGACACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  2:335214/211‑1 (MQ=255)
                                                                                                                             tgcaggaatgaacttcgtGCCTCCG‑TTATGGGTTGTATTATGAGCCTG‑CCCAAC‑AAAA‑TTGTTATGACTTCCTATAC‑CTGGATA‑CCCACCAGTTTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:185623/211‑1 (MQ=255)
                                                                                                                             gtcaaaaaTTGACTGGTGTCCACCT‑TTAATGTTTGTTTTATGAGCCTG‑CCCCAC‑AAAA‑TTGTTATGCCTTCCTATGT‑CTTGATA‑TCCACCAGATTAAATAAGATTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:756604/222‑1 (MQ=255)
                                                                                                                             ctaagtaattgagtagttGCCACTG‑TTAAGTGTTGTATTATCAGCCTG‑CCCCAA‑AAAT‑TTGTTATGACTTCCTATAT‑CTGGAAA‑ACCCACAGATTCAATAAGATTTCAATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTTTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:502943/211‑1 (MQ=255)
                                                                                                                              gttctaattttattttcctcACCGTTTAATGTTTGTATTATTCGCCCG‑CCCCAC‑AAAA‑TTGTTACGCGTCCCTATAT‑CTTGTGA‑TCCACCAGTTTAAAAAAGACTCCAATTAGTTCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  2:519150/209‑1 (MQ=255)
                                                                                                                                        aCTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                               >  1:473223/1‑131 (MQ=255)
                                                                                                                                        aCTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTagtagt                                                                             >  1:654017/1‑229 (MQ=255)
                                                                                                                                             cTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATc                                                                                                                                      >  1:494870/1‑167 (MQ=255)
                                                                                                                                              tGCCACCG‑TTAATGGTTGTATTATTAGACTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTCCCTCGTTCATCCACTTTTTGGTTGACATGTACCTTCTCATCGTTTACAACACCAATAACATTAGTAGTGacgcc                                                                       >  2:42329/1‑225 (MQ=255)
                                                                                                                                                       ttAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATcc                                                               >  1:42603/1‑229 (MQ=255)
                                                                                                                                                             tttGTATTATTACCCTC‑CCCAAA‑AAAA‑TTGTTATGCTTTCCT‑TATCCTTAATA‑CCCCCCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                         <  1:522330/228‑1 (MQ=255)
                                                                                                                                                             gttttttttttccccTG‑CCCCCC‑AAAA‑TGTTTATCACTTCCTATTT‑CTTAATA‑ACCCCCATATTCAAAAAGATTTCCATTAGAACCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                         <  2:654016/216‑1 (MQ=255)
                                                                                                                                                              ttttttttttttcccctccccAC‑AAAATTTGTTTAGCTTTCCTATCT‑ATTGAAA‑TCCCCGAGTTTAAGAAATATTTCAATGATATCCGCCGGCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggc                                                          <  2:218500/212‑1 (MQ=255)
                                                                                                                                                               tGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGCCATCCATCCcctgcggg                                                       >  1:545180/1‑229 (MQ=255)
                                                                                                                                                               tGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGATGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCACTCTACAGTCTCATCGTTTACCACACCACTAACATTATTAGTCAATGCAGCCACCCcggggtgg                                                       >  1:756157/1‑225 (MQ=255)
                                                                                                                                                                  ttttatgtccTG‑CCCCAC‑AAAA‑TTGTTAT‑TCTTCCCATAT‑GTTGGTATTCCCGCGGATCAAAGAAGTATTCCATTGGTTTCCCCGTCATGAACATTCTGCTTGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                    <  1:352706/221‑1 (MQ=255)
                                                                                                                                                                  attattAGCCTG‑CCCCAC‑AAAT‑TTGTTATGCCTTCCTATAT‑CTTGATA‑CCCACCAGATTCAATAAGATTCCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                    <  1:430026/229‑1 (MQ=255)
                                                                                                                                                                       taGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGATGATAATATGAGTGTTCTTTGTTGTACCTCGTTCATCCACTTTTTGTTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCggggggggaaaaaat                                               >  2:458160/1‑222 (MQ=255)
                                                                                                                                                                        aGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCACTAACATTAGTAGTGAAGGCAGCCATCCCGGGAGGGGCatatata                                              >  2:197712/1‑228 (MQ=255)
                                                                                                                                                                             tcccccAA‑AAAA‑TTTTATTGCCTTCCTATTT‑TTTGATC‑ACCCAAAGTTCAAATAAGGTTTCCTTGAAGTCCCCCGTCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAgg                                          <  2:738909/227‑1 (MQ=255)
                                                                                                                                                                             g‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTATTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAgga                                         >  1:622963/1‑228 (MQ=255)
                                                                                                                                                                               ccccAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTACTAGTGCAGGCAGCCTTCCCGGCCGGGGCATATATCCAGGCa                                        >  1:293568/1‑229 (MQ=255)
                                                                                                                                                                                          ttttttttgccTTCCTATAT‑CTTGATA‑CCCCCCAGATTCAAAAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAAt                               <  1:458160/219‑1 (MQ=255)
                                                                                                                                                                                              ttATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                  >  2:268032/1‑94 (MQ=255)
                                                                                                                                                                                              ttATCACTCCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                  <  1:268032/94‑1 (MQ=255)
                                                                                                                                                                                               tATGCCTTCATATAC‑CATGATA‑TCCCCAAGAGCTAATAAGCTTTCCTTGTGAAACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGa                          <  1:494833/229‑1 (MQ=255)
                                                                                                                                                                                               tATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAAATATAAAAAAGGa                          >  2:102281/1‑229 (MQ=255)
                                                                                                                                                                                                      tCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGGGAAGGCAGCCACCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAgacca                   >  2:569044/1‑225 (MQ=255)
                                                                                                                                                                                                           tat‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGTCGGGGCATATATCACGGCGGATATCATTAAGGAAAGTACTGAgcg              >  1:240914/1‑229 (MQ=255)
                                                                                                                                                                                                            tt‑CTAGATA‑TCCCTCGGATTAAAGAAGTC‑TCGATTGGCTACCCCGGCATGAACATTCTGCTGGCCATGGTTTATAATATGCGTGTTATTTGTTGTCCCTCGTTATTCCACTTTTTGGTTGCCATCCACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:202376/228‑1 (MQ=255)
                                                                                                                                                                                                             t‑CTTTATA‑TCCACAAATTTAAAAAAGATGTCCATTAGTTCCCCCGCCATGAACATTCTGCTGGCCATGTTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:215677/229‑1 (MQ=255)
                                                                                                                                                                                                              ttttGATA‑CCCCCCAGATTCAATAAACTTTCCTTTGGTTACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:112901/227‑1 (MQ=255)
                                                                                                                                                                                                                   ggtaaccacgaGATTTAATAAGTCTTCCATTAGATACCCCTCCATGAACATTCTGCTGGCCATGGTTTATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACGTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACa        <  1:688809/224‑1 (MQ=255)
                                                                                                                                                                                                                   ata‑tCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGCTAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGCAAGTACTGAGCGGCGACaa       >  1:663080/1‑229 (MQ=255)
                                                                                                                                                                                                                         caccaGATTCAATAAGTTTTCTATCAGCTCCCCCGCCATGAACATTTTGCTTGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCAGCTACAGTCGCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAataata  <  2:622962/229‑1 (MQ=255)
                                                                                                                                                                                                                         caccaGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATtagta                                                                               >  2:144022/1‑152 (MQ=255)
                                                                                                                                                                                                                         caccaGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGCGAAGGCAGCCACCCCGGGCGGGGCATATATAAAGCTAGCTATCAATAAGGAAAGCACCGAGCGGCGACaataaaa  >  2:112901/1‑227 (MQ=255)
                                                                                                                                                                                                                          accaGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCCGGCGGGGCATATATCAAGGCAGATATCAATAAGGATAGTACTGATCGGCGGCagtaat   >  1:513354/1‑227 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   
ATTCCATTCGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCG‑TTAATGGTTGTATTATTAGCCTG‑CCCCAC‑AAAA‑TTGTTATGACTTCCTATAT‑CTTGATA‑TCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATA  >  NC_000913/2782267‑2782706

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