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

ATTCCATTCGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATA  >  NC_000913/2782267‑2782706
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   
tttCCATTCGGCCCCCCCCCCTCTACGCGTGAAGCAACGCCCTTAATTTGCGTTCCATTGCTGATCCCCCCTTTATCAACGTCCTGATGGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCAGTCAGGAATTGACTTTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:382716/228‑1 (MQ=255)
tttCCATTCCCCCCCCCTCCCCCTCCGCGTGAAGCCCCGCCCTTAATTGTTGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTTCCTCTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:487335/228‑1 (MQ=255)
         gCACTCACTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTCATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCACCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAACTTGTTATGACTTCCCATTTCTTTTTATCCACCAGCTTCAATAACAATTCCATTaa                                                                                                                                                                                                            >  2:502943/1‑228 (MQ=255)
                       tACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCAACCGTTAATGGTTGTATTTTTAGCCTGCCCCACAAAATTGTTAAGACTTCCTATTACTTTAAAACCACCAGGTTCACAAAGATTTCCATAAAAACCCCCGccctaa                                                                                                                                                                                              >  2:313582/1‑225 (MQ=255)
                         cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGCTGAATTATTAGCCTTCCCCACAAAATTGTTATTACTTCCAATATCTTGTTATCATCCAGATTCAATAAGACTTCCATTAGATACcccccctttaa                                                                                                                                                                                             >  1:341145/1‑223 (MQ=255)
                         cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGCCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATGCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTACCCTGCCCCACAAAATTGTTATGACTCCCTATATCTTGATTTCACCCAGATTCAAAAAAATTTCCATTAAAAAACCCGCCATGaa                                                                                                                                                                                             >  2:350525/1‑228 (MQ=255)
                          ccgTGAAGCCCCGCCCTTAATTGTCGTTGCATTGCTGATCCCCCCTTTATAAACGTCCTTATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCAAAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACa                                                                                                                                                                                           <  1:680512/228‑1 (MQ=255)
                               aaGCCCCGCCCTTAATTGTCGTTCCATTGCTGAGACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTGCCTATATCTTGATACCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  1:544526/198‑1 (MQ=255)
                               aaGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAACCGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCAGCCGTTAATGTTTGTATTATTAGCCTGCCCCACACAATTGTTATGACTTCCTTTATCTTGATATCCCCCAGTATCAATCAGATCTCCTTTAGAtctgtttcttatatcaccttg                                                                                                                                                                                      >  2:544526/1‑209 (MQ=255)
                                             tttGTCGTTCCATTGCGTATACCCCCTTTAAAAACGTCCTGTTTGCCCCTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCCCCTTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGCCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                          <  1:176438/226‑1 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCACCCAGATTCAATAAGACTTACATTAGCTACCCCGCCATGAACATTCTTCTGGCCATGGTTGataat                                                                                                                                                                    >  2:247549/1‑229 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTAGGACTTCCTATTTCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCAGGAAAATTCTCCTGGTCATGGTTgtaaaa                                                                                                                                                                    >  2:606978/1‑224 (MQ=255)
                                                 tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTCCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGatcat                                                                                                                                                                    >  1:380817/1‑229 (MQ=255)
                                                     tCCGTTTCTGATCCCCCCTTAATAACCGCCCTGATTCCCACTCTCGTTTGTCGTACCTGTGGAAATACCCCTGTAATGAATTGATTTTCTGCCACCTTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGCCTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                                                                                                                                                <  1:269969/229‑1 (MQ=255)
                                                          ttttttcccccctTTTAAAACCCTCCTTATTGCCCTTTTTCTTTGTTGTCCCTGTGTAAATCCCCCCGTCAGTAATTGAATGTTTTCCACCTTTAATGGTTGTTTTATTCGCCTCCCCACAAAAATTGTTATGTCTTCATATATCTTGATATCCACCAGTTTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                           <  2:302701/221‑1 (MQ=255)
                                                          tttgaatccccccTTACTAAGCTTCCTGATTGCACCCCCCGTTTGTCGTACTTGTGGAAATCCCCCTACCATATATTGACTTTTGGCCACCGTTATTGCTTGTATTATTATCTTGCCCAAAAAATTTGTTATGTCTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                           <  1:312215/221‑1 (MQ=255)
                                                               aTACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  1:737188/166‑1 (MQ=255)
                                                               aTACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     >  2:737187/1‑166 (MQ=255)
                                                                 aCCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGAAACCCCTCCATGAACATTCGTCTGGCCATGGTTGACAATATGGTGGTTAttttt                                                                                                                                                     >  2:688808/1‑226 (MQ=255)
                                                                               cttCCTGATGGCCCCTCTCTATTGTCGTCCCTGTGAAAATCCCCCTGTCCTGAATTTACTGTCTGCCACCGTTAATGGTTGTATTATTACCCTGCCCCACAAAATTGTTATGACTTCCTATATCTGGATATCCACCAGATTCAATAAGACTTCCATGAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATc                                                                                                                                      <  1:451821/227‑1 (MQ=255)
                                                                                 tCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCCATAAGACTTCCCTTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTTCTAATATGAGTGTTCTTTGTTGTACCTCGTTAATCCa                                                                                                                                    >  1:738910/1‑229 (MQ=255)
                                                                                 tCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGAATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATATTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTCACGGCCATGGTTGATAATATGATTGTTATTTGTTGTACCTCGTTCATcc                                                                                                                                     >  1:349085/1‑228 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCTTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGCTTCACTAAGACTTCCCTTAGATACCCCCCCATGAACATTCTGCTGGCCATGATTGCTAATATCAGTGTTTTTTTTTGTACCCCGttaaccccc                                                                                                                                   >  2:215677/1‑221 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     >  1:493867/1‑146 (MQ=255)
                                                                                   cTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAc                                                                                                                                                                                                                     <  2:493867/146‑1 (MQ=255)
                                                                                                    ttttcgacgtttgtaaaaTCCCCCCTGTATGAATTGATTTTCTGCCACCATTAATGGTTGTATTATTGGCCGGCCCAACAAAATTATTATGACTTCCTTTATCTTGATTACCACCAGATTCAATAAGACTTCATTTAGATTCCCCGCCATGATCATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTAc                                                                                                                 <  1:313582/214‑1 (MQ=255)
                                                                                                         gttctggtgGAATTACCGCCGTCATGATGTTCTTTTTTTCCACCTTTAAGGTTTGTATTATTACCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                            <  1:606979/223‑1 (MQ=255)
                                                                                                            ctttggtattcccccccGTCATATATTAATTGTCTCCCCCTTTAAATGGTTATTTTATTACCCCTCCCCAAAAAATTGTTATGACTTCCTATATCTTGATACCCAACAGTTTAAATAAGACTTCCATTAGATTCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCAt                                                                                                         <  2:124778/217‑1 (MQ=255)
                                                                                                                          cccGCCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGCTCATCCACTTTGTGGTTGCCATCTACAGTCTCAACGTTTACCACACCa                                                                                           >  2:408158/1‑229 (MQ=255)
                                                                                                                            tatcttaaattaatttcttGCCACGGTTATTGGTTGATTATTAACCCGCCCCCAAAAAATTGTTAGAACTCCCTGTATTCAGATATCCACCAAATTAAATAAGATTTCTATTAGACACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  2:335214/211‑1 (MQ=255)
                                                                                                                             tgcaggaatgaacttcgtGCCTCCGTTATGGGTTGTATTATGAGCCTGCCCAACAAAATTGTTATGACTTCCTATACCTGGATACCCACCAGTTTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:185623/211‑1 (MQ=255)
                                                                                                                             gtcaaaaaTTGACTGGTGTCCACCTTTAATGTTTGTTTTATGAGCCTGCCCCACAAAATTGTTATGCCTTCCTATGTCTTGATATCCACCAGATTAAATAAGATTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:756604/222‑1 (MQ=255)
                                                                                                                             ctaagtaattgagtagttGCCACTGTTAAGTGTTGTATTATCAGCCTGCCCCAAAAATTTGTTATGACTTCCTATATCTGGAAAACCCACAGATTCAATAAGATTTCAATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTTTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  1:502943/211‑1 (MQ=255)
                                                                                                                              gttctaattttattttcctcACCGTTTAATGTTTGTATTATTCGCCCGCCCCACAAAATTGTTACGCGTCCCTATATCTTGTGATCCACCAGTTTAAAAAAGACTCCAATTAGTTCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                         <  2:519150/209‑1 (MQ=255)
                                                                                                                                        aCTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                               >  1:473223/1‑131 (MQ=255)
                                                                                                                                        aCTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTagtagt                                                                             >  1:654017/1‑229 (MQ=255)
                                                                                                                                             cTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATc                                                                                                                                      >  1:494870/1‑167 (MQ=255)
                                                                                                                                              tGCCACCGTTAATGGTTGTATTATTAGACTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTCCCTCGTTCATCCACTTTTTGGTTGACATGTACCTTCTCATCGTTTACAACACCAATAACATTAGTAGTGacgcc                                                                       >  2:42329/1‑225 (MQ=255)
                                                                                                                                                       ttAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATcc                                                               >  1:42603/1‑229 (MQ=255)
                                                                                                                                                             tttGTATTATTACCCTCCCCAAAAAAATTGTTATGCTTTCCTTATCCTTAATACCCCCCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                         <  1:522330/228‑1 (MQ=255)
                                                                                                                                                             gttttttttttccccTGCCCCCCAAAATGTTTATCACTTCCTATTTCTTAATAACCCCCATATTCAAAAAGATTTCCATTAGAACCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                         <  2:654016/216‑1 (MQ=255)
                                                                                                                                                              ttttttttttttcccctccccACAAAATTTGTTTAGCTTTCCTATCTATTGAAATCCCCGAGTTTAAGAAATATTTCAATGATATCCGCCGGCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggc                                                          <  2:218500/212‑1 (MQ=255)
                                                                                                                                                               tGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGCCATCCATCCcctgcggg                                                       >  1:545180/1‑229 (MQ=255)
                                                                                                                                                               tGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGATGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCACTCTACAGTCTCATCGTTTACCACACCACTAACATTATTAGTCAATGCAGCCACCCcggggtgg                                                       >  1:756157/1‑225 (MQ=255)
                                                                                                                                                                  ttttatgtccTGCCCCACAAAATTGTTATTCTTCCCATATGTTGGTATTCCCGCGGATCAAAGAAGTATTCCATTGGTTTCCCCGTCATGAACATTCTGCTTGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                    <  1:352706/221‑1 (MQ=255)
                                                                                                                                                                  attattAGCCTGCCCCACAAATTTGTTATGCCTTCCTATATCTTGATACCCACCAGATTCAATAAGATTCCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                    <  1:430026/229‑1 (MQ=255)
                                                                                                                                                                       taGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGATGATAATATGAGTGTTCTTTGTTGTACCTCGTTCATCCACTTTTTGTTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCggggggggaaaaaat                                               >  2:458160/1‑222 (MQ=255)
                                                                                                                                                                        aGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCACTAACATTAGTAGTGAAGGCAGCCATCCCGGGAGGGGCatatata                                              >  2:197712/1‑228 (MQ=255)
                                                                                                                                                                             tcccccAAAAAATTTTATTGCCTTCCTATTTTTTGATCACCCAAAGTTCAAATAAGGTTTCCTTGAAGTCCCCCGTCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAgg                                          <  2:738909/227‑1 (MQ=255)
                                                                                                                                                                             gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTATTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAgga                                         >  1:622963/1‑228 (MQ=255)
                                                                                                                                                                               ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTACTAGTGCAGGCAGCCTTCCCGGCCGGGGCATATATCCAGGCa                                        >  1:293568/1‑229 (MQ=255)
                                                                                                                                                                                          ttttttttgccTTCCTATATCTTGATACCCCCCAGATTCAAAAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAAt                               <  1:458160/219‑1 (MQ=255)
                                                                                                                                                                                              ttATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                  >  2:268032/1‑94 (MQ=255)
                                                                                                                                                                                              ttATCACTCCCTATATCTTGATATCCACCAGATTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                  <  1:268032/94‑1 (MQ=255)
                                                                                                                                                                                               tATGCCTTCATATACCATGATATCCCCAAGAGCTAATAAGCTTTCCTTGTGAAACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGa                          <  1:494833/229‑1 (MQ=255)
                                                                                                                                                                                               tATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAAATATAAAAAAGGa                          >  2:102281/1‑229 (MQ=255)
                                                                                                                                                                                                      tCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGGGAAGGCAGCCACCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAgacca                   >  2:569044/1‑225 (MQ=255)
                                                                                                                                                                                                           tatCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGTCGGGGCATATATCACGGCGGATATCATTAAGGAAAGTACTGAgcg              >  1:240914/1‑229 (MQ=255)
                                                                                                                                                                                                            ttCTAGATATCCCTCGGATTAAAGAAGTCTCGATTGGCTACCCCGGCATGAACATTCTGCTGGCCATGGTTTATAATATGCGTGTTATTTGTTGTCCCTCGTTATTCCACTTTTTGGTTGCCATCCACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:202376/228‑1 (MQ=255)
                                                                                                                                                                                                             tCTTTATATCCACAAATTTAAAAAAGATGTCCATTAGTTCCCCCGCCATGAACATTCTGCTGGCCATGTTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:215677/229‑1 (MQ=255)
                                                                                                                                                                                                              ttttGATACCCCCCAGATTCAATAAACTTTCCTTTGGTTACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc            <  1:112901/227‑1 (MQ=255)
                                                                                                                                                                                                                   ggtaaccacgaGATTTAATAAGTCTTCCATTAGATACCCCTCCATGAACATTCTGCTGGCCATGGTTTATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACGTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACa        <  1:688809/224‑1 (MQ=255)
                                                                                                                                                                                                                   atatCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGCTAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGCAAGTACTGAGCGGCGACaa       >  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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                   
ATTCCATTCGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATA  >  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: