Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,223 T→G V139G (GTG→GGG)  fimH → type 1 fimbriae D‑mannose specific adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,2230TG93.5% 112.6 / ‑1.1 46V139G (GTG→GGG) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (16/27);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCAC  >  NC_000913/4548976‑4549469
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
gtgtaTCTTTCTACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATTGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGc                                                                                                                                                                                                                                                        <  1:97692/247‑1 (MQ=255)
gTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAGCCGGAAACCATTACAGACTATGTCACACAGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAGCGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGc                                                                                                                                                                                                                                                        <  1:502946/251‑1 (MQ=255)
  ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTACTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAACACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCgg                                                                                                                                                                                                                                                      >  1:115314/1‑250 (MQ=255)
         gcGACGCAAATCTGTGGCCATATCGATTAGCCGGAAACCATTACAGACTATGTAACACTGCAACGAGGATCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTTACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAg                                                                                                                                                                                                                                                <  2:419954/249‑1 (MQ=255)
              gaaaaTCGTTTTCCCTAACGATTATCCTGAAACCATTACAGACTAGGTCACACTGCAACGAGGCTCGGCGTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCt                                                                                                                                                                                                                                          <  2:322959/249‑1 (MQ=255)
               cAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGCGCTTTTTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTc                                                                                                                                                                                                                                         >  1:265457/1‑251 (MQ=255)
                      tttGCCATAACTATTATCCGGAAACCATTACAGGCTATGTCACACTGCAACGAGGCTCGGCTTATTGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTTCCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTa                                                                                                                                                                                                                                     <  1:85651/248‑1 (MQ=255)
                      tttGCCATAACGATTATCCGGAAACCCTTCCAGCCTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAAtt                                                                                                                                                                                                                                  >  1:391785/1‑251 (MQ=255)
                          ccATATCGATTATCCGGAAACCATGTCATACTATGTCACACTGCAGCGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCAATCCATTTCCTGCCACCAGCGAAACGCCTCGCGTTGTTTATAATACGAGAACGGATAAGCCGTTTCCGGTGGCGCTTTAGTTGACGCCTTTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCg                                                                                                                                                                                                                              <  1:473468/251‑1 (MQ=255)
                                        cGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCg                                                                                                                                                                                                                 >  1:210117/1‑250 (MQ=255)
                                            aaCCATTACAGACTATGTCACACTGCAACGAGTCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGa                                                                                                                                                                                                            <  2:491116/251‑1 (MQ=255)
                                             aCCATTCCAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGGAGGTATCCATTTCATAACACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAACCGGTGGCGGTGGCGCCTTATTTTCCGCCTGAGAGCAGGAGGGGGGGGGGGGCGTTTAAAGCTGGATAATTTATTTCAGTGCTTATTTTgttacagg                                                                                                                                                                                                            >  2:210822/1‑243 (MQ=255)
                                                                                      tGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCCCCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGGCCAACAACTATAACAGCGATGCTTTACAGtttagttcggaa                                                                                                                                                                  >  2:105519/1‑241 (MQ=255)
                                                                                      tGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGttt                                                                                                                                                                           <  1:105519/241‑1 (MQ=255)
                                                                                          ggcgTGTTATCTAATTTTTCCGGGAACGTAAAATATAGTTGCAGTAGCTATCCATTTCCTACCACCAGCGAAACTCCGCGCGTTGTTTATAATTCGAGAACTGATAAGCCGTGGGCGGTTGCGCTTTAGTTGACGCCTGTGAGCCGTGCGGGCGGGGGGTCGATTAAAGCTGGCTCATTAAGTGCCGTGCTTATTTTGCGACAGATCAACAACTAGAACAGCGATGATTTCCAGTTTGTGTGGCATTTTTa                                                                                                                                                             <  2:381145/251‑1 (MQ=255)
                                                                                            cgTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACg                                                                                                                                                           <  2:398459/251‑1 (MQ=255)
                                                                                                       aTTTTGCCGTGACTGTAAAATATAGTTGCAGTAGCTATACATTTCCTACCACCAGCGAAACGCCGCGTGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTAGTTGTCGCCTGTGAGCAGTGCGGGTGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCATTAACTATAACAGCGACGATTTCCAGTTTGTGTGGAATATTTACGCCAATAatgat                                                                                                                                                <  1:57793/251‑1 (MQ=255)
                                                                                                             ccGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATgtggtg                                                                                                                                          <  1:406109/251‑1 (MQ=255)
                                                                                                                  aCCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGcc                                                                                                                                     <  1:292785/251‑1 (MQ=255)
                                                                                                                    cgtaaaagtaaGTGGCAGTTGCTTTCCATTTCCTACCACCAGCGAAACGCCGCGCGGTGTTTAGTATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTACGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGGTTGTGTGGAATATTTCCGCCAATAATGATGTGGTCGTGCCTa                                                                                                                                   <  2:468149/241‑1 (MQ=255)
                                                                                                                    cGTAACATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTCAGCCAATAATGATGTGGTGGTGCCt                                                                                                                                    >  2:36644/1‑250 (MQ=255)
                                                                                                                            ataGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAAt                                                                                                                                                                                                                                   <  2:16927/148‑1 (MQ=255)
                                                                                                                            ataGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAAt                                                                                                                                                                                                                                   >  1:16927/1‑148 (MQ=255)
                                                                                                                                  gcaggagctcaccagttCCTACCACCAGCGCAACGCCGCGCTTGGTTGAGAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTGTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGCCTGCGAt                                                                                                                     <  2:29028/236‑1 (MQ=255)
                                                                                                                                    aGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGCTGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGAGGCcgcggttt                                                                                                                   >  2:86827/1‑247 (MQ=255)
                                                                                                                                       aGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTg                                                                                                                                                                                                                   <  2:277117/153‑1 (MQ=255)
                                                                                                                                       aGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTg                                                                                                                                                                                                                   >  1:277117/1‑153 (MQ=255)
                                                                                                                                                                            ttgttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTTACCGTTAATCTGCCGGACTACCCt                                                                          >  1:348850/1‑251 (MQ=255)
                                                                                                                                                                                  ataATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTc                                                                     <  1:86827/250‑1 (MQ=255)
                                                                                                                                                                                    aaTTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGt                                                                  <  1:207040/251‑1 (MQ=255)
                                                                                                                                                                                    aaTTCGAGAACGGATAAGCCGGGGCCGGTGGCGCTTGATTTGGCGCCGGTGAGCAGTGCTGGCGGGGGGGCGATTAAAGCTGGTTCATTAATTGCCGTGCTTATTGTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGCGCCTACTGGCGGCTGCGAGGTTTCTGCTCGGGATGTCACCGTTACTCTGCCGGACGGCCCTGGTTCAg                                                                   <  2:168853/250‑1 (MQ=255)
                                                                                                                                                                                     aTTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTTCCGGACTACCCTGGTTCAGTg                                                                 >  1:341738/1‑251 (MQ=255)
                                                                                                                                                                                               ggATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGAGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGGGCTTATTTTGCGACAGACAACAAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCCATAATGATGTGGTGGTTCCTCCTGGCGGCTGGGTTGTTTCTGCTCGTGATTGTCACGTTACTTTACCCGAGTCCCCTGGTTCAggtgacatttccc                                                       >  2:375239/1‑239 (MQ=255)
                                                                                                                                                                                                   aaGCCGTGGCCGGTTGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTa                                                    <  2:249376/250‑1 (MQ=255)
                                                                                                                                                                                                   aaGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTa                                                    <  1:356925/250‑1 (MQ=255)
                                                                                                                                                                                                     gCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCGCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCg                                                <  2:482222/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGTGGCCGGTGGCTCTTTATTTGACGACTGTGCGCAGTGCGTGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGGTTGTGTGGAATATTTACGCCAATAATGATGTGGTTGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGt                                               <  2:426371/251‑1 (MQ=255)
                                                                                                                                                                                                        gTGGCCGGTGGCGCTGTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGTCGATTAAAGCTTGCGCATTAATTGCCGTGCTTATTTTGCGACAGCCCAACAACTATAACAGCGAGGATGTCCAGTTTTTGTGGAATATTTACGTCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTTCCCTGGTTCAGTGCCAATTTCTCTTACCGttt                                             <  1:41067/251‑1 (MQ=255)
                                                                                                                                                                                                               ggTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATtgtg                                       >  1:232413/1‑251 (MQ=255)
                                                                                                                                                                                                                          agttGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCaa                            <  1:36644/249‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGTGAGCAGAGCGGGCGGGGGGGCGATTAAAGCTGGCTAATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGGTGGCATGTTTCTTCTCGTGATGTCAACGTTCATATGCCGGACTACCCTGGTTCATTGCCCAATGCTATTACCCGTTTTTGTGAGAACAGCCAAAAccatggg                   >  2:460542/1‑246 (MQ=255)
                                                                                                                                                                                                                                      tgAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGCTGATTTCCAGTTTGTGTGGACTATTTACGCCAATAATGATGTGGTGGTGCCTCCTGGCGGCTGCGATGTTTCTGCTCGTGTTGTCCCCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTTCCGTTTATTGTGCGAAAAGCCCAAACCCgggggct                >  1:235486/1‑248 (MQ=255)
                                                                                                                                                                                                                                      tgAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCCTTACTCTTCCGGACCACCCTGGTTCATTGCCAATTCCCCTTTCCGTTTTTTTTGCGAAAAGCCCAAACCTgggggat                >  2:486983/1‑248 (MQ=255)
                                                                                                                                                                                                                                       gAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGTCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGTTGTCACCGTTACTCTGCCGGGCTCCCCTGGTTCAGTGCCAATTCCTCTTTCCGTTTTTTGGGCGAAAAGCCAAAACCTGGGGGAtt               >  2:150027/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTAcc            <  1:394695/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gcgggcggGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCtctc        <  1:248208/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggggggGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGGCACCGTTAATCTGGCGGAATACCCTGGGTCAGTGCCAATTCCTCTTAACGTTTATTGTGCGAAAAGCCAAAACCTGGGGGATTAACTCTCCGGcac  >  1:363811/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
GTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCAC  >  NC_000913/4548976‑4549469

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