Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 559,176 T→A 16.8% V68E (GTA→GAA)  sfmC → putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913559,1760TA16.8% 115.5 / 5.0 49V68E (GTA→GAA) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base T (34/6);  new base A (2/7);  total (36/13)
Fisher's exact test for biased strand distribution p-value = 5.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  NC_000913/558945‑559414
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ccGTCAAGAATTCACCAGGCTTTTTAACTATTATGACTAAAATAAAGTTAGTGATTCTCATTTTATTTTATTTAATCATTTCGGCCAGCGCCCAAGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTaaaagaa                                                                                                                                                                                                                                         <  2:878277/239‑1 (MQ=255)
         aTACCCAAGGCGTTTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  2:642786/238‑1 (MQ=255)
              caAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATTCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTcatcat                                                                                                                                                                                                                           <  1:935376/239‑1 (MQ=255)
                      aTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  1:234665/225‑1 (MQ=255)
                                atgaCTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAACTAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGGGTAAAAGAAAAGTCATTCATCATTACACCGCCACtgtttgt                                                                                                                                                                                                         >  1:510862/1‑239 (MQ=255)
                                     tAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCATGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTAAA‑TTCGTGGATTGAAAACAGCAGCGGTGTCAAAGAAAAGcaattca                                                                                                                                                                                                                               >  1:337849/1‑212 (MQ=255)
                                        aaTAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCCTGCGGTATTGATTAGACATAGCCATACCAATGAGCGCTTTCTGGCCAA‑TTCGTGGATTGAAAACAGCAGCGCTGTAAAAGAAAAGTCATTCCTCATTACCCCGCCACTGTTTGTTAGAGaa                                                                                                                                                                                                  >  2:218412/1‑238 (MQ=255)
                                             aGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTAATTAATCATTACACCGCCACTGTTTGTTAGGAAACCCaaa                                                                                                                                                                                            >  2:580641/1‑239 (MQ=255)
                                                                    tatttaATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtattat                                                                                                                                                                     >  1:512416/1‑239 (MQ=255)
                                                                         aaTCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAATCGAAAAAACTTTGCGAATTATTTac                                                                                                                                                                 >  1:1012316/1‑238 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTAGACCGCCACTGTTTGTTAGTGAACACAAAAACGAAAAAAATTTACTTATAATTTACAcc                                                                                                                                                              >  2:462499/1‑239 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCCCCACCGTTGTTTAGTGAACCCAAAAACAAAAAAACATTCGGAATTATTTAAAcc                                                                                                                                                              >  2:422718/1‑239 (MQ=255)
                                                                                tCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACAGCACCGTTTTGTTAGTGAACCCAAAAGCAAAAATACTTTGCGCATTATTTAAACAGGTcc                                                                                                                                                         >  2:396844/1‑239 (MQ=255)
                                                                                     cccccccctggctgcGGGAGGTGTCGTATGGGGTGCCGGGTTTATTATTTCCCCCCATGATGAAAACCATGCGGCGGTATGGATTAAAAATAGCCATCCAAAGGAGCCTTTTCTGGTAAA‑TTCGTGGATTGAAAACAGCAGCGGTAAAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGc                                                                                                                                                    <  2:541189/229‑1 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTGAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCCAGTTCTGGTCAA‑TTCGAGGATTGAAAACAGCAGCGGTGAAAAAGAAAAGTCATTCATCATTACACCTCAACTGTTTGTTAGTGAACCCAAAAGCGAAAAACCTTTGGGAATACTTTAAACAGGCCCACCCCGGgaagca                                                                                                                                            >  2:549988/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTACAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCCCTGgcagca                                                                                                                                            >  1:685708/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATCATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGGGAAAATACTTTGCGaattat                                                                                                                                                                     >  1:723445/1‑214 (MQ=255)
                                                                                                    cGGAGGCACCGCATGAGGGACCACGCTTATTATTTACCCCGCTGATGCAAAACAGACTGCGGTAGGGATTAGAAATAGCAATACAAATGAGCGCTTTCTGGACAA‑TTCGGGGATTGAAAACAGCAGCGGTTAAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTg                                                                                                                                     <  1:218412/239‑1 (MQ=255)
                                                                                                           aTCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCACCGGGGACAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCAAAAAGCGAAAATACTTTTGGTATTTTTTAAACCGGTCCCCCGCTGgcagcaactccaggatc                                                                                                                                 >  2:946343/1‑225 (MQ=255)
                                                                                                                    ggTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATCCTTGGCGTATTATTTAAACCGGTCCACCCCTGGCAGCAGATCGGGAGTTTTTGTTTTGGAt                                                                                                                     >  2:581827/1‑239 (MQ=255)
                                                                                                                     gTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTTGGTATTATTTACACCGGTCCACCGCGGGCGGCCAACCGTGAGTCTCTGTTTTGGGTg                                                                                                                    >  2:23258/1‑239 (MQ=255)
                                                                                                                           ggggttttttttacCCCGTTGTTGCTAAACAGATTGCGGTATGGATTAAAAAAAGCCATACAAATGGGCGTTTTCTGTTCAA‑TTCGGGGATTGAAAACAGCAGCGGTGAAAAAGAAAAGTCATCAATATTTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGtt                                                                                                              <  1:422718/230‑1 (MQ=255)
                                                                                                                                attattTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAgtc                                                                                                         >  1:539917/1‑237 (MQ=255)
                                                                                                                                      tatCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTTAACCCAAAAGCGAAAATACTTTGCGACTTATTTACACCGGTCACCCGCTGGCGGAAGATCGTGGGTTCCTGTTCTGGATGAATGTTAAAAGGATc                                                                                                     >  2:663237/1‑237 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  1:721949/1‑239 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  2:310216/1‑239 (MQ=255)
                                                                                                                                         cccGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTc                                                                                                >  1:176396/1‑239 (MQ=255)
                                                                                                                                             cTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTa                                                                                            >  1:1059627/1‑239 (MQ=255)
                                                                                                                                              tGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAg                                                                                           >  1:288935/1‑239 (MQ=255)
                                                                                                                                                    aaaaCAGACGGGGGTATGGATTAAAAATAGCCATCCCAATGAGCGTTTTCTGGTCAA‑TGCGTGGATTAAAAACAGCAGCGTTGTAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATaaaa                                                                                     <  1:23258/238‑1 (MQ=255)
                                                                                                                                                          gACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGGATTATTTACACCGGTCCACCGCGTGCAGAAGATCGGGAGTCTATGTTCTGGATGAATGTTAAAACGATCCCTGCGGGAGATAAAAAAGCAt                                                                               >  2:474820/1‑239 (MQ=255)
                                                                                                                                                             tGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTAAACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTTTGGCTGAATGTTAAGACGACCCCTCCGGGAGATAAAAATGCAttaa                                                                            >  2:563696/1‑237 (MQ=255)
                                                                                                                                                                  tttgtataaaaaaTAGCCAACCAAATGGGTGTTTTTGTGTAAA‑TTCGGGAATTAAAACCAGAAGGGGTAAAAAAAAAAGTTCATTCACCATTACCCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  2:288935/228‑1 (MQ=255)
                                                                                                                                                                  tttGGATTAAAAAGAGCCAAACAATGGAGCGCTTTCT‑GTCAATTTCGTGGATTGAAAACGGCAGGGGTGAAAAAAAAAAGTAATACATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  1:751638/236‑1 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGaa                                                                  >  1:997208/1‑239 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTATTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAAAGGCCGGaa                                                                  >  1:389613/1‑239 (MQ=255)
                                                                                                                                                                                aGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCa                                                         >  1:1015822/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTGGATAAAAATGCATTGAACGACAGGAATGTTTTGCaa                                                        >  1:104855/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGAACCCTTCGGTCGATAAAAATGCATTGAACGGCAGGGAATTTTTGCaa                                                        >  2:154567/1‑239 (MQ=255)
                                                                                                                                                                                           tGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGTTAGGTCAAAATGCATTGAACGGCAGGAATGTTTTGCAAAGGGCGAttt                                              >  2:657496/1‑239 (MQ=255)
                                                                                                                                                                                               cgcTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAAGGCATTTAACGGCAGGAATGTTTTGCAAATGGGGATTTTATc                                          >  2:395473/1‑239 (MQ=255)
                                                                                                                                                                                                  tttCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAACGGCAGGAAATTTTTGCAACTGGGGTTTTTATcccc                                       >  2:575809/1‑236 (MQ=255)
                                                                                                                                                                                                    tgtggaaaa‑TTTGTGCATGAAAAACATCAGGGGTAAAAAAAAAAAGTCATTCATCTTTCCTCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCa                                     <  1:563696/232‑1 (MQ=255)
                                                                                                                                                                                                        gttttttttgtgtATTAAAAACAGCGGGGGTGAAAAAAAAAATGTAGTCACCATCACCCCGCCCCTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGa                                  <  2:510862/230‑1 (MQ=255)
                                                                                                                                                                                                        gTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGaa                                 >  2:876431/1‑239 (MQ=255)
                                                                                                                                                                                                            a‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGCAAGGAATGTTTTGCAACTGGCGATTTTACCGCGCATGAAatta                             >  1:121291/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGCTCCCTTCGGTAGATAAAACTCCATTGAACGGCAGGAATGTTTTGCAACTGGGCATTTTAACGCGCATGAAATTAtttat                        >  1:439891/1‑236 (MQ=255)
                                                                                                                                                                                                                                   cGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAAc                                                                          >  2:353194/1‑172 (MQ=255)
                                                                                                                                                                                                                                   cGGTGTAAAAAAAAATTCATTCACCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAAc                                                                          <  1:353194/172‑1 (MQ=255)
                                                                                                                                                                                                                                        tAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGaa  >  1:665045/1‑239 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAA‑TTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  NC_000913/558945‑559414

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