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

CGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGTCCTGTTCATAAATTTCATCCTCTGAAAAGT  >  NC_000913/2782275‑2782745
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
cGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTCTTATTAGCCTGCCCCACAAAATTGTTTTGACTTCCTATATCTTGATATCCACCAAATTCAAAAAAAC‑TTCCATTTAATCcccccc                                                                                                                                                                                                                                          >  2:660935/1‑237 (MQ=255)
             tccccgcgtcacccccctcccTTATTTGTCTTTCCTTTGCTGTTCCCTCTTTTAAAACCGTCT‑TGATTGCCTCTCTGGATTTGCGTACCTGTGGAAATCCCCCCGTTAGAAAGT‑GACTGTCGTCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCAATACCTTGATATCCACCAGATTAAATAAGCC‑TTCCATTAGATACCCCGCCATGAACATTctg                                                                                                                                                                                                                             <  2:519129/221‑1 (MQ=255)
                 cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTctcctg                                                                                                                                                                                                                          >  1:940466/1‑238 (MQ=255)
                 cgcgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAA‑TTCCATTAGATACCCCGCCATGAACATTctgctg                                                                                                                                                                                                                          >  1:983663/1‑238 (MQ=255)
                   cgtaaaacgccGCCTTTTGTTGTCTTTCCATTCTTTAT‑CCCCCTTTTTAAACCTCCGTTTTTCCCACTCTCGATGGTCGTACCTGTGGAAACACCCCCGTCATGAATT‑GCATGTCTGCCACCGTTAAT‑GGTTTTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGcc                                                                                                                                                                                                                       <  2:386778/230‑1 (MQ=255)
                   cgTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACGGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGcc                                                                                                                                                                                                                       >  1:327880/1‑239 (MQ=255)
                    ttaaccccccccccTTATTTTTCCTTCCTTTTTTTATCCCCCTTTTATAACCTTCC‑TGTTGCCCCCTCTGGATTGTTGTCCCTGTGGAAATACCCCCGTAATGAATT‑GCCTGTCTGCCCCCGTTAAT‑GGTTGTATTATTAGCCTCCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGCC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      <  2:1054004/228‑1 (MQ=255)
                       aaGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGTTATCCACAAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATgg                                                                                                                                                                                                                   >  1:590873/1‑239 (MQ=255)
                        atcccccactttAATTGCCTTTCCCATGGTCATCCCCCCTTTATAACCGGCC‑TGTTGCCCACCCTCGATTGTCGTGCCTGGGAAAATCCCCCCGTCATGAATT‑GACTGTCTCCCACCGTTAAT‑GGTTGTATGATTACCCTGCCCCACAAAATGTTTATAACTTCCTATATCTTGATACCCACCAGTTTCAATAAGAT‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATgg                                                                                                                                                                                                                   <  2:651298/228‑1 (MQ=255)
                         ccccccccttttaTTGTCTTTCCTTTGCTTATCCCCCCTTTAAAACCGTCC‑TTATTGCCCCTTCGGATTTGTGTCCTTGGGGAAATCCCCCCGTCATGAATT‑GACTGTCTGCCACCTTTAAG‑GGTTTTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATACCCACCAAATTCAATAAGAC‑TCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                                                                 <  2:350242/227‑1 (MQ=255)
                         accccccctttttttttctttCCGTTGCGGATACCCCCTTTATAACCGTCC‑GGATTGCCCCTCGCGATTGTCGTACCGGGGGAAATCCCCCCTCAATGAATT‑GGTTGTGTGCCCCCTTTAAT‑GTTTGTATTATTGCCCTGCCCTAAAAATTTTTTTTGACTTCCTATATCTTGATATCCCCCAGATTAAATAAGAC‑TTCCATTAGAAACCCCGCCATGAACATTCTGCTGGCCATGGtt                                                                                                                                                                                                                 <  2:627731/220‑1 (MQ=255)
                            ccGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCCCAAAATTGTTATGACTTCCTATTTCTTGCTATCCACCAGATTCAATAAAAC‑TTCCATTAGATACCCCGCCATGAACATCTTGCTGGCCATGGTTGaa                                                                                                                                                                                                              >  2:139996/1‑238 (MQ=255)
                                 cttttttGTCTTTCCATTGCTGATCCCCCCTTTAAAACCGCCC‑GTATGGCATCTCTCGATTGTCGTACTTGTGGAAATACCCTTGTCATGAATG‑TACTGTCTCCCCCCTTTAAT‑GGTTGTATTATTGGCCTGCCCCACAAAATTGTTATTCCTTCTTATATCTTGATTTCCACCAGATTCAATAACAC‑TTCCATTAGCTACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                                                         <  2:709624/234‑1 (MQ=255)
                                 ctgttttttgtttCCTTTGTGGATCCCCCTTTTATAATCGTCC‑TCATGCCCCCTCCGGATGGTCGTCCCTGTGGAATCTCCCCCTTCAAGACTT‑GAGTTGCTGCCGCGTTTAAT‑GGTTTTATTATTAGCCTGCCCCACAAATTTGTTATGACTTCTTATATTTTGATATCCACCAGTTTCAAAAAGAT‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                                                                                                                                         <  2:102710/228‑1 (MQ=255)
                                   tAATTGTCGTTCAATTGCTGATCCCCCCTTTAAAACCGTCC‑TGATTGCCCCTCTCGTTGGTCGTCCCTGTGAAAATACCCCCTTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCAAAAAATTGTTATGACTTCCTATATCTTAATACCCACCAGATTAAATAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                                                                                                                                                                                       <  1:468394/239‑1 (MQ=255)
                                        ttttttccattctctatcccctcTTTAAAATCGTCC‑TATTCCCCCCTCTTGTTTGCTGTACCTGTGGAAATCCCCCCGTCATGAATT‑GCTTGTCTGCCACCGTTAAT‑GTTTTTATTATTAGCCGCCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCATATTCAATAAGCT‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  2:636972/221‑1 (MQ=255)
                                        gTCGTTCCTTTGCTGATCCCCCCTTAAAAACCGTCC‑TGTTTGCCACTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTAAATAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  1:29356/239‑1 (MQ=255)
                                        cttgTTCCATTTGTGATCCCCCTTTTAAAAACGTCC‑TGATTGCCCCTCTCGATTGTCTTACCTGTGGAAATACCCCCGTCATGAATT‑GACTTTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATATCCCCCATATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                                                                                                                                  <  1:727920/236‑1 (MQ=255)
                                         tCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTCCTGGCCATGGTTGATAAAATGAggtttt                                                                                                                                                                                                 >  2:468394/1‑234 (MQ=255)
                                          cGTTCCAGTGCTGATCCCCCCTTTAAAACCGTCC‑TGATTGCCACTCGCGATTGTCGTCCCTGTGGAAATCCCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAt                                                                                                                                                                                                <  1:340335/239‑1 (MQ=255)
                                           gTTCCATTGCTGATACCCCTTTTAAAAACGTCC‑TGATTGCCACTCTGTATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGATTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAt                                                                                                                                                                                                <  1:334488/238‑1 (MQ=255)
                                                 ttGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCCTGGTTGATAATATGAGGGTTATTtgttgt                                                                                                                                                                                         >  2:340335/1‑239 (MQ=255)
                                                 ttGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAACC‑TTCCATTAGCTACCCCCCCATGAACATTCTCCTGGCCATTGGTGATAAAATAAGGGTTTTttgttt                                                                                                                                                                                          >  2:970212/1‑237 (MQ=255)
                                                     tGATTCCCCTTTTATAAACTTTC‑TTCTTACCCTTCTGATTTGTTGTACGTGGGGAATCTCCCCTCTAATGATTG‑GACGTTCTGCCACTGTTAAT‑GGTTTTATTAATAGCCCGCCCCACAAAATTGTTATGACTGCCTATATCTTGATAACCTCAAGAGTTAAAAAGAC‑CTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCt                                                                                                                                                                                     <  2:439493/239‑1 (MQ=255)
                                                           ccccTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCa                                                                                                                                                                               >  1:43705/1‑239 (MQ=255)
                                                              cTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGTTAATATGAGTGTTATTTGTTTTACCTCGTTCTTcc                                                                                                                                                                            >  1:269776/1‑239 (MQ=255)
                                                                     aaCGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAAAAAAAC‑TTCCATTAGATACCCCGCCCTGAACATTCTGCTGGCCAGGGTTGATAATATGAGGGTTCTTTTTTTTACCTCGTTACTCACCttttt                                                                                                                                                                     >  2:925343/1‑239 (MQ=255)
                                                                         tCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGTTAATATGAGGGTTATTTTTTGTACCTCGTTCATCCCCTTTTTGGtt                                                                                                                                                                 >  2:151207/1‑239 (MQ=255)
                                                                         tCC‑TGATTGCAACTCTAGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCTGACAGGATCATTCTTCTGGCCATGGTTGTAAATATGGAAGCGATTTGATGTACCTAGTTCATGCACTTTTTgga                                                                                                                                                                  >  1:157119/1‑237 (MQ=255)
                                                                             tGATTGCCCCTCTCGATTTCCGTCCCTTGGGAAATCCCCCCGTAATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                                                                                                                              <  1:935273/239‑1 (MQ=255)
                                                                              taTTGCCACTCGCGATTGTCGTATCTGTGGAAATCCCCCCGTCATGAATT‑GACTGTCTGCCACCTTTAAG‑GTTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCa                                                                                                                                                             <  1:139996/238‑1 (MQ=255)
                                                                                ttGCCCCTCTCGATTGTCGTACGTGTGGAAATCCCCCCTTCATGAATT‑AACTGTCTGCCCCCTTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATc                                                                                                                                                           <  1:797925/239‑1 (MQ=255)
                                                                                  gCCACTCTCTATTGCTGTCCTTGTGGAAGTACCCCCTTCATGAATT‑GACTGTTGGCCCCCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                                                                                                                                         <  1:595285/239‑1 (MQ=255)
                                                                                  gCCACTCTCGATTGTCGTACCTGTGGAAATCCCCCCGTCATGAATT‑GCCTTTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTTTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                                                                                                                                         <  1:627281/239‑1 (MQ=255)
                                                                                        tctttttttcttttctttgGAAATCCCCCCTTCTTTATTT‑ACTTGTCTGCCACCTTTATT‑GGTTGTTTTATCAGCCGCCCCCCAAAAATTGTTATGCTTTCCTATTCCTTGATATCCACCAGATTAAAAAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                                                                   <  2:421004/222‑1 (MQ=255)
                                                                                         tcGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCTTCCACTTTTGGGTTTCCCTCTACAGtctc                                                                                                                                                  >  1:927062/1‑239 (MQ=255)
                                                                                         tagaggggctgcccTGTGAAAATACCCCCGTTATGAATT‑GACTTTTTTCCA‑CGTTAATGGGTTGTATTATTAGCCTGCCCCACAAAATGTTTATGACTTCCTATATTTTGATACCCACCAGATTAAATAAGAC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                                                                                                                                   <  2:671327/226‑1 (MQ=255)
                                                                                                tCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACt                                                                                                                                                                         >  1:884047/1‑209 (MQ=255)
                                                                                                 cGTACCTGTGGAAATCCCCCCGTCATGAATT‑GACTGTCTTCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCAATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTa                                                                                                                                          <  1:992623/239‑1 (MQ=255)
                                                                                                    aCCTGTGGAAAACCCCCCGTAATGAATT‑GACTGTCTGCCACCTTTAAT‑GGTTGTTTTATTAGCCTGCCCCAAAAAATGGTTATGACTGCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCAATAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACca                                                                                                                                       <  1:970212/239‑1 (MQ=255)
                                                                                                         tgGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTAc                                                                                                                                         >  1:671327/1‑232 (MQ=255)
                                                                                                          ggAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTTCCATCTAAAGTCTCAACGTTTACCACACCaa                                                                                                                                 >  2:104530/1‑239 (MQ=255)
                                                                                                            aattACCCCCTGTATGAATG‑GACTGTCTTCACCCTTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATTCAATAAGCC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:117268/235‑1 (MQ=255)
                                                                                                            aaaTCCCCCCGTCATGAATT‑GCCTGTTCGCCCCCGTTAAG‑GTTTGTTTTATTAGCCTGCCCCAAAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:151207/238‑1 (MQ=255)
                                                                                                             aaTACCCCCCTTATGTATTGGACGGTTGTCCCCCGTTAAT‑TGTTTTATTTTTCGCCTCCCCCACAAAATTTTTAGGACTTCCTATTTCTTGATATCCCCCATATTCAATAAGTC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  2:327880/238‑1 (MQ=255)
                                                                                                                aCCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGGC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      <  1:897810/148‑1 (MQ=255)
                                                                                                                aCCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                                                                                                                                                                                      >  2:897810/1‑148 (MQ=255)
                                                                                                                      gTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATTAGTGTTATTTGTTGTACCTCGTTCATCCCCTTTTTGGTTGACATCTACAGTCTCATCGTTTccccacccaaaaaaattggaag                                                                                                                     >  2:138805/1‑220 (MQ=255)
                                                                                                                                  aCTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGcc                                                                                                          >  1:246365/1‑239 (MQ=255)
                                                                                                                                      tCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCTTc                                                                                                       >  1:913562/1‑238 (MQ=255)
                                                                                                                                            ccggTTATT‑TGTTGTATTATTAGCCCCCCCCAAAAAATTGTTTGGCTTTCCTATATCGTGATATCCACTAGATTCAAAAAGCC‑TGCCATAAATTTCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACGTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                                                                <  2:983663/236‑1 (MQ=255)
                                                                                                                                             cctttttttgtttttttttttaCCCTCCCCCACAAATTTTTTTTATTTCCCTATATCTGTATATCACCCTGTTTGAATAAGAC‑TTCCTTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCcgggcg                                                                                                <  2:378470/220‑1 (MQ=255)
                                                                                                                                                ttAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGAC‑TTCAATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAAt                                                                                                                                <  1:604342/203‑1 (MQ=255)
                                                                                                                                                ttAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTCCCACACCAAt                                                                                                                                >  2:604342/1‑203 (MQ=255)
                                                                                                                                                 tAAG‑GGTTTTATTATTGCCAGGCCCCACAAAATGTTTATGACTTCATATATTTGTATATCCCCCAAATTCAATAAGCT‑TTCCATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCa                                                                                           <  1:524266/239‑1 (MQ=255)
                                                                                                                                                           tattatTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAATATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATAAAg                                                                                  >  2:643324/1‑239 (MQ=255)
                                                                                                                                                                ttaGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAgc                                                                             >  1:29808/1‑238 (MQ=255)
                                                                                                                                                                  aGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGata                                                                           >  1:748455/1‑239 (MQ=255)
                                                                                                                                                                        ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATa                                                                     >  1:228990/1‑239 (MQ=255)
                                                                                                                                                                            aaaaaaTTTTTTTTATTTCTTATATCTTAATACCCCCCA‑ATTCAAAAAGACTTTCCTTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGa                                                                 <  2:246365/237‑1 (MQ=255)
                                                                                                                                                                                aaTTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTATATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtc                                                                                                                                                    <  2:417741/152‑1 (MQ=255)
                                                                                                                                                                                aaTTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtc                                                                                                                                                    >  1:417741/1‑152 (MQ=255)
                                                                                                                                                                                          ccTTCCTAATTCTTGAATCCCCCCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:660935/238‑1 (MQ=255)
                                                                                                                                                                                          aCTTCCTATGTCTTGATACCCACCAAATTCAATAAGAC‑TTCCATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:925343/239‑1 (MQ=255)
                                                                                                                                                                                          aCTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGa                                                                                                                                                                                                               >  1:102668/1‑83 (MQ=255)
                                                                                                                                                                                          aCTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGa                                                                                                                                                                                                               <  2:102668/83‑1 (MQ=255)
                                                                                                                                                                                           ctcccTATATCTTTATACCCACCAGGTTCAAAAAGCC‑TTCTATTAAATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAgcggc                                                   <  1:32030/235‑1 (MQ=255)
                                                                                                                                                                                             tCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCACTAACATTAGTAGTGACGGCAGCCATCCCGGCCGGGGCATATATAAAGTCAGATATCAATAAGGAAAGTAATGAGGGGGGAc                                                >  2:648594/1‑239 (MQ=255)
                                                                                                                                                                                                     gatggtaataagagaaagcAATAAGCC‑TTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAataata                                         <  1:842043/220‑1 (MQ=255)
                                                                                                                                                                                                      ttGATATCCCCCAGATTCAATAAGAC‑TTCCATTAAATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   <  1:156084/115‑1 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTCTTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   >  2:156084/1‑115 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   >  2:689453/1‑115 (MQ=255)
                                                                                                                                                                                                      ttGATATCCACCAGATTAAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                                                                                                                                   <  1:689453/115‑1 (MQ=255)
                                                                                                                                                                                                         atatCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAAAAAGGAAAGAACTGAGCGGCGCCAAAAAAAGGGAc                                    >  2:965422/1‑239 (MQ=255)
                                                                                                                                                                                                              caccaGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATCAGTAGTGACGCCAGCCATCCCGGCCGGGGCCTATATAAAGGAAGATTTCAATAAGGACAGTACTTAGCGACGACAATCATAGTGAATgg                                 >  1:1024347/1‑237 (MQ=255)
                                                                                                                                                                                                                              aC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAAAAGGGACTGGTCCTGTTCATAAAtttaa               >  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)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
CGCACTCCCTCCCTCTACGCGTGAAGCACCGCCCTTAATTGTCGTTCCATTGCTGATACCCCCTTTATAAACGTCC‑TGATTGCCACTCTCGATTGTCGTACCTGTGGAAATACCCCCGTCATGAATT‑GACTGTCTGCCACCGTTAAT‑GGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGAC‑TTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAGCCATCCCGGGCGGGGCATATATCAAGGCAGATATCAATAAGGAAAGTACTGAGCGGCGACAATAATAGGGACTGGTCCTGTTCATAAATTTCATCCTCTGAAAAGT  >  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: ‑