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

CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  NC_000913/558945‑559414
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ccGTCAAGAATTCACCAGGCTTTTTAACTATTATGACTAAAATAAAGTTAGTGATTCTCATTTTATTTTATTTAATCATTTCGGCCAGCGCCCAAGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTaaaagaa                                                                                                                                                                                                                                         <  2:878277/239‑1 (MQ=255)
         aTACCCAAGGCGTTTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  2:642786/238‑1 (MQ=255)
              caAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATTCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTcatcat                                                                                                                                                                                                                           <  1:935376/239‑1 (MQ=255)
                      aTTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcatt                                                                                                                                                                                                                                 <  1:234665/225‑1 (MQ=255)
                                atgaCTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAACTAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGGGTAAAAGAAAAGTCATTCATCATTACACCGCCACtgtttgt                                                                                                                                                                                                         >  1:510862/1‑239 (MQ=255)
                                     tAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCATGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTAAATTCGTGGATTGAAAACAGCAGCGGTGTCAAAGAAAAGcaattca                                                                                                                                                                                                                               >  1:337849/1‑212 (MQ=255)
                                        aaTAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCCTGCGGTATTGATTAGACATAGCCATACCAATGAGCGCTTTCTGGCCAATTCGTGGATTGAAAACAGCAGCGCTGTAAAAGAAAAGTCATTCCTCATTACCCCGCCACTGTTTGTTAGAGaa                                                                                                                                                                                                  >  2:218412/1‑238 (MQ=255)
                                             aGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTAATTAATCATTACACCGCCACTGTTTGTTAGGAAACCCaaa                                                                                                                                                                                            >  2:580641/1‑239 (MQ=255)
                                                                    tatttaATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtattat                                                                                                                                                                     >  1:512416/1‑239 (MQ=255)
                                                                         aaTCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAATCGAAAAAACTTTGCGAATTATTTac                                                                                                                                                                 >  1:1012316/1‑238 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTAGACCGCCACTGTTTGTTAGTGAACACAAAAACGAAAAAAATTTACTTATAATTTACAcc                                                                                                                                                              >  2:462499/1‑239 (MQ=255)
                                                                           tCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCCCCACCGTTGTTTAGTGAACCCAAAAACAAAAAAACATTCGGAATTATTTAAAcc                                                                                                                                                              >  2:422718/1‑239 (MQ=255)
                                                                                tCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACAGCACCGTTTTGTTAGTGAACCCAAAAGCAAAAATACTTTGCGCATTATTTAAACAGGTcc                                                                                                                                                         >  2:396844/1‑239 (MQ=255)
                                                                                     cccccccctggctgcGGGAGGTGTCGTATGGGGTGCCGGGTTTATTATTTCCCCCCATGATGAAAACCATGCGGCGGTATGGATTAAAAATAGCCATCCAAAGGAGCCTTTTCTGGTAAATTCGTGGATTGAAAACAGCAGCGGTAAAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGc                                                                                                                                                    <  2:541189/229‑1 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTGAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCCAGTTCTGGTCAATTCGAGGATTGAAAACAGCAGCGGTGAAAAAGAAAAGTCATTCATCATTACACCTCAACTGTTTGTTAGTGAACCCAAAAGCGAAAAACCTTTGGGAATACTTTAAACAGGCCCACCCCGGgaagca                                                                                                                                            >  2:549988/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTACAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCCCTGgcagca                                                                                                                                            >  1:685708/1‑239 (MQ=255)
                                                                                             aTGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATCATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGGGAAAATACTTTGCGaattat                                                                                                                                                                     >  1:723445/1‑214 (MQ=255)
                                                                                                    cGGAGGCACCGCATGAGGGACCACGCTTATTATTTACCCCGCTGATGCAAAACAGACTGCGGTAGGGATTAGAAATAGCAATACAAATGAGCGCTTTCTGGACAATTCGGGGATTGAAAACAGCAGCGGTTAAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTg                                                                                                                                     <  1:218412/239‑1 (MQ=255)
                                                                                                           aTCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCACCGGGGACAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCAAAAAGCGAAAATACTTTTGGTATTTTTTAAACCGGTCCCCCGCTGgcagcaactccaggatc                                                                                                                                 >  2:946343/1‑225 (MQ=255)
                                                                                                                    ggTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATCCTTGGCGTATTATTTAAACCGGTCCACCCCTGGCAGCAGATCGGGAGTTTTTGTTTTGGAt                                                                                                                     >  2:581827/1‑239 (MQ=255)
                                                                                                                     gTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTTGGTATTATTTACACCGGTCCACCGCGGGCGGCCAACCGTGAGTCTCTGTTTTGGGTg                                                                                                                    >  2:23258/1‑239 (MQ=255)
                                                                                                                           ggggttttttttacCCCGTTGTTGCTAAACAGATTGCGGTATGGATTAAAAAAAGCCATACAAATGGGCGTTTTCTGTTCAATTCGGGGATTGAAAACAGCAGCGGTGAAAAAGAAAAGTCATCAATATTTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGtt                                                                                                              <  1:422718/230‑1 (MQ=255)
                                                                                                                                attattTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAgtc                                                                                                         >  1:539917/1‑237 (MQ=255)
                                                                                                                                      tatCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTTAACCCAAAAGCGAAAATACTTTGCGACTTATTTACACCGGTCACCCGCTGGCGGAAGATCGTGGGTTCCTGTTCTGGATGAATGTTAAAAGGATc                                                                                                     >  2:663237/1‑237 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  1:721949/1‑239 (MQ=255)
                                                                                                                                       atCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCt                                                                                                  >  2:310216/1‑239 (MQ=255)
                                                                                                                                         cccGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTc                                                                                                >  1:176396/1‑239 (MQ=255)
                                                                                                                                             cTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTa                                                                                            >  1:1059627/1‑239 (MQ=255)
                                                                                                                                              tGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAg                                                                                           >  1:288935/1‑239 (MQ=255)
                                                                                                                                                    aaaaCAGACGGGGGTATGGATTAAAAATAGCCATCCCAATGAGCGTTTTCTGGTCAATGCGTGGATTAAAAACAGCAGCGTTGTAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATaaaa                                                                                     <  1:23258/238‑1 (MQ=255)
                                                                                                                                                          gACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGGATTATTTACACCGGTCCACCGCGTGCAGAAGATCGGGAGTCTATGTTCTGGATGAATGTTAAAACGATCCCTGCGGGAGATAAAAAAGCAt                                                                               >  2:474820/1‑239 (MQ=255)
                                                                                                                                                             tGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTAAACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTTTGGCTGAATGTTAAGACGACCCCTCCGGGAGATAAAAATGCAttaa                                                                            >  2:563696/1‑237 (MQ=255)
                                                                                                                                                                  tttgtataaaaaaTAGCCAACCAAATGGGTGTTTTTGTGTAAATTCGGGAATTAAAACCAGAAGGGGTAAAAAAAAAAGTTCATTCACCATTACCCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  2:288935/228‑1 (MQ=255)
                                                                                                                                                                  tttGGATTAAAAAGAGCCAAACAATGGAGCGCTTTCTGTCAATTTCGTGGATTGAAAACGGCAGGGGTGAAAAAAAAAAGTAATACATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACgg                                                                        <  1:751638/236‑1 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGaa                                                                  >  1:997208/1‑239 (MQ=255)
                                                                                                                                                                       aTTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTATTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAAAGGCCGGaa                                                                  >  1:389613/1‑239 (MQ=255)
                                                                                                                                                                                aGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCa                                                         >  1:1015822/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTGGATAAAAATGCATTGAACGACAGGAATGTTTTGCaa                                                        >  1:104855/1‑239 (MQ=255)
                                                                                                                                                                                 gCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCAACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGAACCCTTCGGTCGATAAAAATGCATTGAACGGCAGGGAATTTTTGCaa                                                        >  2:154567/1‑239 (MQ=255)
                                                                                                                                                                                           tGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGTTAGGTCAAAATGCATTGAACGGCAGGAATGTTTTGCAAAGGGCGAttt                                              >  2:657496/1‑239 (MQ=255)
                                                                                                                                                                                               cgcTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAAGGCATTTAACGGCAGGAATGTTTTGCAAATGGGGATTTTATc                                          >  2:395473/1‑239 (MQ=255)
                                                                                                                                                                                                  tttCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGAATTGAACGGCAGGAAATTTTTGCAACTGGGGTTTTTATcccc                                       >  2:575809/1‑236 (MQ=255)
                                                                                                                                                                                                    tgtggaaaaTTTGTGCATGAAAAACATCAGGGGTAAAAAAAAAAAGTCATTCATCTTTCCTCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCa                                     <  1:563696/232‑1 (MQ=255)
                                                                                                                                                                                                        gttttttttgtgtATTAAAAACAGCGGGGGTGAAAAAAAAAATGTAGTCACCATCACCCCGCCCCTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGa                                  <  2:510862/230‑1 (MQ=255)
                                                                                                                                                                                                        gTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGaa                                 >  2:876431/1‑239 (MQ=255)
                                                                                                                                                                                                            aTTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGCAAGGAATGTTTTGCAACTGGCGATTTTACCGCGCATGAAatta                             >  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)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
CCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAA  >  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: