Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,782,507 A→C 17.8% V73G (GTA→GGA)  ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,782,5070AC17.8% 176.2 / 2.8 81V73G (GTA→GGA) ypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (37/25);  new base C (2/16);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 3.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGTCCTGTTCATAAATTTCATCCTCTGAAAAGT  >  NC_000913/2782275‑2782745
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
cGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTCTTATTAGCCTGCCCCACAAAATTGTTTTGACTTCCTATATCTTGATATCCACCAAATTCAAAAAAACTTCCATTTAATCcccccc                                                                                                                                                                                                                                          >  2:660935/1‑237 (MQ=255)
             tccccgcgtcacccccctcccTTATTTGTCTTTCCTTTGCTGTTCCCTCTTTTAAAACCGTCTTGATTGCCTCTCTGGATTTGCGTACCTGTGGAAATCCCCCCGTTAGAAAGTGACTGTCGTCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCAATACCTTGATATCCACCAGATTAAATAAGCCTTCCATTAGATACCCCGCCATGAACATTctg                                                                                                                                                                                                                             <  2:519129/221‑1 (MQ=255)
                 cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctcctg                                                                                                                                                                                                                          >  1:940466/1‑238 (MQ=255)
                 cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAATTCCATTAGATACCCCGCCATGAACATTctgctg                                                                                                                                                                                                                          >  1:983663/1‑238 (MQ=255)
                   cgtaaaacgccGCCTTTTGTTGTCTTTCCATTCTTTATCCCCCTTTTTAAACCTCCGTTTTTCCCACTCTCGATGGTCGTACCTGTGGAAACACCCCCGTCATGAATTGCATGTCTGCCACCGTTAATGGTTTTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGcc                                                                                                                                                                                                                       <  2:386778/230‑1 (MQ=255)
                   cgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACGGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGcc                                                                                                                                                                                                                       >  1:327880/1‑239 (MQ=255)
                    ttaaccccccccccTTATTTTTCCTTCCTTTTTTTATCCCCCTTTTATAACCTTCCTGTTGCCCCCTCTGGATTGTTGTCCCTGTGGAAATACCCCCGTAATGAATTGCCTGTCTGCCCCCGTTAATGGTTGTATTATTAGCCTCCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGCCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      <  2:1054004/228‑1 (MQ=255)
                       aaGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGTTATCCACAAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATgg                                                                                                                                                                                                                   >  1:590873/1‑239 (MQ=255)
                        atcccccactttAATTGCCTTTCCCATGGTCATCCCCCCTTTATAACCGGCCTGTTGCCCACCCTCGATTGTCGTGCCTGGGAAAATCCCCCCGTCATGAATTGACTGTCTCCCACCGTTAATGGTTGTATGATTACCCTGCCCCACAAAATGTTTATAACTTCCTATATCTTGATACCCACCAGTTTCAATAAGATTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATgg                                                                                                                                                                                                                   <  2:651298/228‑1 (MQ=255)
                         ccccccccttttaTTGTCTTTCCTTTGCTTATCCCCCCTTTAAAACCGTCCTTATTGCCCCTTCGGATTTGTGTCCTTGGGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCTTTAAGGGTTTTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATACCCACCAAATTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                                                                 <  2:350242/227‑1 (MQ=255)
                         accccccctttttttttctttCCGTTGCGGATACCCCCTTTATAACCGTCCGGATTGCCCCTCGCGATTGTCGTACCGGGGGAAATCCCCCCTCAATGAATTGGTTGTGTGCCCCCTTTAATGTTTGTATTATTGCCCTGCCCTAAAAATTTTTTTTGACTTCCTATATCTTGATATCCCCCAGATTAAATAAGACTTCCATTAGAAACCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                                                                 <  2:627731/220‑1 (MQ=255)
                            ccGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCCCAAAATTGTTATGACTTCCTATTTCTTGCTATCCACCAGATTCAATAAAACTTCCATTAGATACCCCGCCATGAACATCTTGCTGGCCATGGTTGaa                                                                                                                                                                                                              >  2:139996/1‑238 (MQ=255)
                                 cttttttGTCTTTCCATTGCTGATCCCCCCTTTAAAACCGCCCGTATGGCATCTCTCGATTGTCGTACTTGTGGAAATACCCTTGTCATGAATGTACTGTCTCCCCCCTTTAATGGTTGTATTATTGGCCTGCCCCACAAAATTGTTATTCCTTCTTATATCTTGATTTCCACCAGATTCAATAACACTTCCATTAGCTACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                                                         <  2:709624/234‑1 (MQ=255)
                                 ctgttttttgtttCCTTTGTGGATCCCCCTTTTATAATCGTCCTCATGCCCCCTCCGGATGGTCGTCCCTGTGGAATCTCCCCCTTCAAGACTTGAGTTGCTGCCGCGTTTAATGGTTTTATTATTAGCCTGCCCCACAAATTTGTTATGACTTCTTATATTTTGATATCCACCAGTTTCAAAAAGATTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                                                         <  2:102710/228‑1 (MQ=255)
                                   tAATTGTCGTTCAATTGCTGATCCCCCCTTTAAAACCGTCCTGATTGCCCCTCTCGTTGGTCGTCCCTGTGAAAATACCCCCTTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCAAAAAATTGTTATGACTTCCTATATCTTAATACCCACCAGATTAAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                                                                                                                                                                                       <  1:468394/239‑1 (MQ=255)
                                        ttttttccattctctatcccctcTTTAAAATCGTCCTATTCCCCCCTCTTGTTTGCTGTACCTGTGGAAATCCCCCCGTCATGAATTGCTTGTCTGCCACCGTTAATGTTTTTATTATTAGCCGCCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCATATTCAATAAGCTTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  2:636972/221‑1 (MQ=255)
                                        gTCGTTCCTTTGCTGATCCCCCCTTAAAAACCGTCCTGTTTGCCACTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTAAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  1:29356/239‑1 (MQ=255)
                                        cttgTTCCATTTGTGATCCCCCTTTTAAAAACGTCCTGATTGCCCCTCTCGATTGTCTTACCTGTGGAAATACCCCCGTCATGAATTGACTTTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATATCCCCCATATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  1:727920/236‑1 (MQ=255)
                                         tCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTCCTGGCCATGGTTGATAAAATGAggtttt                                                                                                                                                                                                 >  2:468394/1‑234 (MQ=255)
                                          cGTTCCAGTGCTGATCCCCCCTTTAAAACCGTCCTGATTGCCACTCGCGATTGTCGTCCCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAt                                                                                                                                                                                                <  1:340335/239‑1 (MQ=255)
                                           gTTCCATTGCTGATACCCCTTTTAAAAACGTCCTGATTGCCACTCTGTATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGATTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAt                                                                                                                                                                                                <  1:334488/238‑1 (MQ=255)
                                                 ttGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCCTGGTTGATAATATGAGGGTTATTtgttgt                                                                                                                                                                                         >  2:340335/1‑239 (MQ=255)
                                                 ttGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAACCTTCCATTAGCTACCCCCCCATGAACATTCTCCTGGCCATTGGTGATAAAATAAGGGTTTTttgttt                                                                                                                                                                                          >  2:970212/1‑237 (MQ=255)
                                                     tGATTCCCCTTTTATAAACTTTCTTCTTACCCTTCTGATTTGTTGTACGTGGGGAATCTCCCCTCTAATGATTGGACGTTCTGCCACTGTTAATGGTTTTATTAATAGCCCGCCCCACAAAATTGTTATGACTGCCTATATCTTGATAACCTCAAGAGTTAAAAAGACCTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCt                                                                                                                                                                                     <  2:439493/239‑1 (MQ=255)
                                                           ccccTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCa                                                                                                                                                                               >  1:43705/1‑239 (MQ=255)
                                                              cTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGTTAATATGAGTGTTATTTGTTTTACCTCGTTCTTcc                                                                                                                                                                            >  1:269776/1‑239 (MQ=255)
                                                                     aaCGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAAAAAAACTTCCATTAGATACCCCGCCCTGAACATTCTGCTGGCCAGGGTTGATAATATGAGGGTTCTTTTTTTTACCTCGTTACTCACCttttt                                                                                                                                                                     >  2:925343/1‑239 (MQ=255)
                                                                         tCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGTTAATATGAGGGTTATTTTTTGTACCTCGTTCATCCCCTTTTTGGtt                                                                                                                                                                 >  2:151207/1‑239 (MQ=255)
                                                                         tCCTGATTGCAACTCTAGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCTGACAGGATCATTCTTCTGGCCATGGTTGTAAATATGGAAGCGATTTGATGTACCTAGTTCATGCACTTTTTgga                                                                                                                                                                  >  1:157119/1‑237 (MQ=255)
                                                                             tGATTGCCCCTCTCGATTTCCGTCCCTTGGGAAATCCCCCCGTAATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                                                                                                                              <  1:935273/239‑1 (MQ=255)
                                                                              taTTGCCACTCGCGATTGTCGTATCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCTTTAAGGTTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCa                                                                                                                                                             <  1:139996/238‑1 (MQ=255)
                                                                                ttGCCCCTCTCGATTGTCGTACGTGTGGAAATCCCCCCTTCATGAATTAACTGTCTGCCCCCTTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATc                                                                                                                                                           <  1:797925/239‑1 (MQ=255)
                                                                                  gCCACTCTCTATTGCTGTCCTTGTGGAAGTACCCCCTTCATGAATTGACTGTTGGCCCCCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                                                                                                                                         <  1:595285/239‑1 (MQ=255)
                                                                                  gCCACTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATTGCCTTTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                                                                                                                                         <  1:627281/239‑1 (MQ=255)
                                                                                        tctttttttcttttctttgGAAATCCCCCCTTCTTTATTTACTTGTCTGCCACCTTTATTGGTTGTTTTATCAGCCGCCCCCCAAAAATTGTTATGCTTTCCTATTCCTTGATATCCACCAGATTAAAAAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                                                                   <  2:421004/222‑1 (MQ=255)
                                                                                         tcGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCTTCCACTTTTGGGTTTCCCTCTACAGtctc                                                                                                                                                  >  1:927062/1‑239 (MQ=255)
                                                                                         tagaggggctgcccTGTGAAAATACCCCCGTTATGAATTGACTTTTTTCCACGTTAATGGGTTGTATTATTAGCCTGCCCCACAAAATGTTTATGACTTCCTATATTTTGATACCCACCAGATTAAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                                                                   <  2:671327/226‑1 (MQ=255)
                                                                                                tCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACt                                                                                                                                                                         >  1:884047/1‑209 (MQ=255)
                                                                                                 cGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTTCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCAATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTa                                                                                                                                          <  1:992623/239‑1 (MQ=255)
                                                                                                    aCCTGTGGAAAACCCCCCGTAATGAATTGACTGTCTGCCACCTTTAATGGTTGTTTTATTAGCCTGCCCCAAAAAATGGTTATGACTGCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCAATAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACca                                                                                                                                       <  1:970212/239‑1 (MQ=255)
                                                                                                         tgGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTAc                                                                                                                                         >  1:671327/1‑232 (MQ=255)
                                                                                                          ggAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTTCCATCTAAAGTCTCAACGTTTACCACACCaa                                                                                                                                 >  2:104530/1‑239 (MQ=255)
                                                                                                            aattACCCCCTGTATGAATGGACTGTCTTCACCCTTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGCCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:117268/235‑1 (MQ=255)
                                                                                                            aaaTCCCCCCGTCATGAATTGCCTGTTCGCCCCCGTTAAGGTTTGTTTTATTAGCCTGCCCCAAAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:151207/238‑1 (MQ=255)
                                                                                                             aaTACCCCCCTTATGTATTGGACGGTTGTCCCCCGTTAATTGTTTTATTTTTCGCCTCCCCCACAAAATTTTTAGGACTTCCTATTTCTTGATATCCCCCATATTCAATAAGTCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  2:327880/238‑1 (MQ=255)
                                                                                                                aCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGGCTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      <  1:897810/148‑1 (MQ=255)
                                                                                                                aCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      >  2:897810/1‑148 (MQ=255)
                                                                                                                      gTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATTAGTGTTATTTGTTGTACCTCGTTCATCCCCTTTTTGGTTGACATCTACAGTCTCATCGTTTccccacccaaaaaaattggaag                                                                                                                     >  2:138805/1‑220 (MQ=255)
                                                                                                                                  aCTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGcc                                                                                                          >  1:246365/1‑239 (MQ=255)
                                                                                                                                      tCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCTTc                                                                                                       >  1:913562/1‑238 (MQ=255)
                                                                                                                                            ccggTTATTTGTTGTATTATTAGCCCCCCCCAAAAAATTGTTTGGCTTTCCTATATCGTGATATCCACTAGATTCAAAAAGCCTGCCATAAATTTCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACGTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                                                                <  2:983663/236‑1 (MQ=255)
                                                                                                                                             cctttttttgtttttttttttaCCCTCCCCCACAAATTTTTTTTATTTCCCTATATCTGTATATCACCCTGTTTGAATAAGACTTCCTTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                                                                <  2:378470/220‑1 (MQ=255)
                                                                                                                                                ttAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCAATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:604342/203‑1 (MQ=255)
                                                                                                                                                ttAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTCCCACACCAAt                                                                                                                                >  2:604342/1‑203 (MQ=255)
                                                                                                                                                 tAAGGGTTTTATTATTGCCAGGCCCCACAAAATGTTTATGACTTCATATATTTGTATATCCCCCAAATTCAATAAGCTTTCCATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                                                           <  1:524266/239‑1 (MQ=255)
                                                                                                                                                           tattatTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAATATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATAAAg                                                                                  >  2:643324/1‑239 (MQ=255)
                                                                                                                                                                ttaGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAgc                                                                             >  1:29808/1‑238 (MQ=255)
                                                                                                                                                                  aGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGata                                                                           >  1:748455/1‑239 (MQ=255)
                                                                                                                                                                        ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATa                                                                     >  1:228990/1‑239 (MQ=255)
                                                                                                                                                                            aaaaaaTTTTTTTTATTTCTTATATCTTAATACCCCCCAATTCAAAAAGACTTTCCTTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGa                                                                 <  2:246365/237‑1 (MQ=255)
                                                                                                                                                                                aaTTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTATATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtc                                                                                                                                                    <  2:417741/152‑1 (MQ=255)
                                                                                                                                                                                aaTTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtc                                                                                                                                                    >  1:417741/1‑152 (MQ=255)
                                                                                                                                                                                          ccTTCCTAATTCTTGAATCCCCCCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:660935/238‑1 (MQ=255)
                                                                                                                                                                                          aCTTCCTATGTCTTGATACCCACCAAATTCAATAAGACTTCCATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:925343/239‑1 (MQ=255)
                                                                                                                                                                                          aCTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGa                                                                                                                                                                                                               >  1:102668/1‑83 (MQ=255)
                                                                                                                                                                                          aCTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGa                                                                                                                                                                                                               <  2:102668/83‑1 (MQ=255)
                                                                                                                                                                                           ctcccTATATCTTTATACCCACCAGGTTCAAAAAGCCTTCTATTAAATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:32030/235‑1 (MQ=255)
                                                                                                                                                                                             tCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCACTAACATTAGTAGTGACGGCAGCCATCCCGGCCGGGGCATATATAAAGTCAGATATCAATAAGGAAAGTAATGAGGGGGGAc                                                >  2:648594/1‑239 (MQ=255)
                                                                                                                                                                                                     gatggtaataagagaaagcAATAAGCCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAataata                                         <  1:842043/220‑1 (MQ=255)
                                                                                                                                                                                                      ttGATATCCCCCAGATTCAATAAGACTTCCATTAAATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   <  1:156084/115‑1 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTCTTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   >  2:156084/1‑115 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   >  2:689453/1‑115 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTAAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   <  1:689453/115‑1 (MQ=255)
                                                                                                                                                                                                         atatCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAAAAAGGAAAGAACTGAGCGGCGCCAAAAAAAGGGAc                                    >  2:965422/1‑239 (MQ=255)
                                                                                                                                                                                                              caccaGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATCAGTAGTGACGCCAGCCATCCCGGCCGGGGCCTATATAAAGGAAGATTTCAATAAGGACAGTACTTAGCGACGACAATCATAGTGAATgg                                 >  1:1024347/1‑237 (MQ=255)
                                                                                                                                                                                                                              aCTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAAAAGGGACTGGTCCTGTTCATAAAtttaa               >  2:258753/1‑237 (MQ=255)
                                                                                                                                                                                                                                   ccATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTCTTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTATCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGAATGAGCCTGTTcctgtcccttataca           >  1:906905/1‑225 (MQ=255)
                                                                                                                                                                                                                                   ccATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGACCTGTTCAAAAATTTCaa              >  2:60354/1‑235 (MQ=255)
                                                                                                                                                                                                                                   ccATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGGACTGAGCGGCGCAAATAAAAGGGACGGGTccctgtctcttaaacaaa              >  2:842043/1‑220 (MQ=255)
                                                                                                                                                                                                                                    cATTAGATCCCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      <  2:946412/35‑1 (MQ=255)
                                                                                                                                                                                                                                         gATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAAAAGGGACTGGCCCTGTTCATAAATTTCATCCCGTGaaa     >  2:752450/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGACAGAACTGGGCGGGGACAataata                                         >  2:568846/1‑200 (MQ=255)
                                                                                                                                                                                                                                            aCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGTCCTGTTCATAAATTTCATCCTCTGAAAAGt  <  1:104530/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
CGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCCTGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGTCCTGTTCATAAATTTCATCCTCTGAAAAGT  >  NC_000913/2782275‑2782745

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