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,2230TG98.0% 134.6 / NA 50V139G (GTG→GGG) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (23/26);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  NC_000913/4548981‑4549464
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
actttcgaagataaTCTTGTGTCATAAAGAGTATGCGTAAACAATAACAGACCATGTCACACTGCAACGAGGCTCGGCTTATGGGGGCGTGTTATCTAATGTTTCCGGGACCGTAAAATATAGTGGTAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTGATAATGCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTGTACGCCTGTGAGCAGTTCGGGCGGGGGGGCGATTa                                                                                                                                                                                                                                           <  2:216974/239‑1 (MQ=255)
    tcgaatcaaagtttttGCCATAACGATTATTCTGAATCCATTACATACTATGTCACACTGCAACGAGGCTCGTCTTATGGCGGCGTGTGATCTAATTTTTCCGGGTCCGTAAAATATAGTTGCAGTAGCTATCCATTTCCTACCTCCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGATCAGTGCGGTCGTGGGGGCGATTAAAg                                                                                                                                                                                                                                        <  2:631737/238‑1 (MQ=255)
         gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGAGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCt                                                                                                                                                                                                                                  >  1:367002/1‑251 (MQ=255)
                 tttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCAtt                                                                                                                                                                                                                              <  1:85511/247‑1 (MQ=255)
                     ccATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCg                                                                                                                                                                                                                      >  1:281009/1‑251 (MQ=255)
                            gtatttCCGAAAACCTTTACAGACAAGGTCACACTGCAAATATGTTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCTGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAt                                                                                                                                                                                                               <  2:372456/246‑1 (MQ=255)
                             aTTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtt                                                                                                                                                                                                              >  1:230457/1‑251 (MQ=255)
                                aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGAGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCCGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGCTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGACTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGTTGGCCCATTAGTTGCCGTGCTTCttttt                                                                                                                                                                                                           >  1:328512/1‑250 (MQ=255)
                                    ggAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCAACAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAAAGGATAAGCCGGGGCCGGTGGCGCTTTATTTGACGCCTGTGGGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCggc                                                                                                                                                                                                       >  1:439893/1‑249 (MQ=255)
                                     gAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACa                                                                                                                                                                                                      <  1:249385/251‑1 (MQ=255)
                                           aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaa                                                                                                                                                                                                <  1:548281/251‑1 (MQ=255)
                                              ccAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGAGGCGTGTTATCTAATTTTTCCGGGCCCGGAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGGTTAAAATTCGAGAACGGATAAGCCGTGGCCAGGGGCGCTTTTTTTGACGCCTGTGAGCAGTGAGGGCGGGGGGGGGGTTAAAGCTGGCTCATTAATTGACACGCTTATTTTggggacgaccaaaa                                                                                                                                                                                              >  2:433895/2‑237 (MQ=255)
                                                       gTCACACTGCAACGAGGCTCGGCGTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCGCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACa                                                                                                                                                                                     <  2:230457/250‑1 (MQ=255)
                                                            acTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCgatg                                                                                                                                                                               >  1:121416/1‑251 (MQ=255)
                                                                  aCGAGGCTCGGCTTATTGCGGTGTGTTTTCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCCCCCTCGAAACGCCGTGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACATCGATGATTTcc                                                                                                                                                                         <  2:597535/251‑1 (MQ=255)
                                                                   cGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGACCGGATAAGCCGTGGCCGGTGGCGCTTTTTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGGAAGAACAACAACTATAACAACGATGATTTccc                                                                                                                                                                        >  1:391829/1‑250 (MQ=255)
                                                                     aGGCTCGGCTTATGGCGGCTTGTTATCTAATTTTTCCGGGACAGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAACCGAGGCCGGTGGCGCTTTATTTGACGCCTGGGAGCAGTGCGGGCGGGGGGGCGAGTAAAGCTGGCTCCATAATTGACGTGCTTTTTTTGTGAGAGACCCACAACtcacaaagcgatgctttcc                                                                                                                                                                         >  2:138165/1‑230 (MQ=255)
                                                                             cTTATGGCGGCGAGTTATCTAATTTTTCCGGGACCCTAAAATATAGTGGCAGTAGCTATCCCTTTACTTCCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACCGCTAAGCCGTGGACCGGGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCAATAATTGCCGTGCGTATTTTTCGACCGACCAACAACTATAAAAAGGATGATTTCCACTTTGTGTgg                                                                                                                                                              >  2:9538/1‑251 (MQ=255)
                                                                               tATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATTGTGGCAGTAGCTATCCATTTCCTACCACCAGCGCAACGCCGCGCGTTGGTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGGGCTTTATTTGGAGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTAAAGGTTGTTCATTAATTTCCGTGTTTATTTTTCGACAGACCAACAACTATAACAGCgattgttttcaggttttggtgt                                                                                                                                                             >  2:41266/1‑231 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTAAAGCTGGATCATTAATTGCCGGGCTTATTTTGCGACAGACca                                                                                                                                                                                                 >  2:522169/1‑190 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGTTCATTAATTGCCGTGCTTATTTTGCGACAGACca                                                                                                                                                                                                 <  1:522169/190‑1 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATgtgg                                                                                                                                     <  1:83300/250‑1 (MQ=255)
                                                                                                        ccGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGCTgtggtg                                                                                                                                   >  2:165190/1‑251 (MQ=255)
                                                                                                                  aaaaTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGGCCGGTGTCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTg                                                                                                                         <  1:288106/251‑1 (MQ=255)
                                                                                                                             gCAGTAGCTATCCATTTCCTACCACCATCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGAt                                                                                                              <  2:455015/251‑1 (MQ=255)
                                                                                                                                  aGCTATCCATTTCCTCCCACCAGCGAAACGCCGAGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACCGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTggcggcatcggtgtttt                                                                                                         >  2:350375/1‑240 (MQ=255)
                                                                                                                                               ctcccaacaGCGAAACGCCGCGCGTTGTTTATCATTCGAGAACGGCTAAGCCGTGGCCGGTAGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACCACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATTATGGTGTGGTGCATACTGGCGGCTGCGATGTTTCTgtgcgttatttt                                                                                            >  2:26163/8‑240 (MQ=255)
                                                                                                                                               cTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATggc                                                                                            >  1:558058/1‑249 (MQ=255)
                                                                                                                                                  ccaccaGCGAAACGCCGCGGGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCAcc                                                                                         <  1:529044/251‑1 (MQ=255)
                                                                                                                                                     ccaGCGAAACGCCGCGCGTTGTTTATAATTCTAGAACGTATAAGCCGTGGTCGGTGGCTCTTAATTTGACGCCTGTGATCAGTGCGGGAGTGTGGGAGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTTCGACAGACCAACAACTATAACAGCGATGATGTCCAGTTTGTGTGGAATATTTACTCCAATAATGATGTGGTGGTGCCTACTGGCGGCTTCGATGTTTCTGCTCGGGATGTTACCGtt                                                                                      <  2:73793/251‑1 (MQ=255)
                                                                                                                                                     caaGCGAAACGCCGCGCGTATTTTATAATACTGGAACGGATAAGCCGTGGTCGGTGGTGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGTGCTATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTGCCAGTTTGTGTGGAATATTTACGCCAATAATGTTGTGGTGGTGCCTTCTGGCGGCTGCGATGTTTCTGCGCGttagggaaccgtt                                                                                      <  2:335824/249‑13 (MQ=255)
                                                                                                                                                       aGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTACTTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGGCTACTGGCGGCTGCGATGGTTCTGCTCGGGATGGCACCGttcc                                                                                    >  1:505889/1‑249 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAGAGCTGGCTCATTAATTTCCGTGCTTATTTTGGGACCGACCAACAACTATAACAGCgatgc                                                                                                                                                                              >  2:132285/1‑139 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCgatga                                                                                                                                                                              <  1:132285/140‑1 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGtt                                                                <  1:9850/250‑1 (MQ=255)
                                                                                                                                                                                  tcgtgaACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGGGCGGGCGGGGGGGCGATTAACGCTGGCTCATTAAGTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTTGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGGTACTCTGCCGGACTACCCTGGTTCAGTGc                                                          <  2:30996/246‑1 (MQ=255)
                                                                                                                                                                                          ggATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCctc                                                 <  1:131184/251‑1 (MQ=255)
                                                                                                                                                                                          ggATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCctc                                                 >  1:515340/1‑251 (MQ=255)
                                                                                                                                                                                                 ccGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGt                                          >  1:29925/1‑251 (MQ=255)
                                                                                                                                                                                                      gCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAtt                                     <  1:239605/251‑1 (MQ=255)
                                                                                                                                                                                                        ccGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGGTGTTTCTGCTCGGGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAttt                                    >  1:121983/1‑250 (MQ=255)
                                                                                                                                                                                                         cGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGGTGCGATGATTCTTATCGGGATTTTACCGGTTCTCTTCCGGACTAACCTGGTTCAGTGCCCATTACTCTTACCCtttttttt                                   >  2:289585/1‑246 (MQ=255)
                                                                                                                                                                                                           gttgttcttttgttGACGCCTGGTAGCCGTGCGGGGGGGGGGGCGATGAAAGCTGGCTCATTAAGTTCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTGTGTGTGGAATATTTACGCCAATAATTATGTTGTGGTGCCTACTGGCGACTGCTATGTTTCTGCTCGTGATGTCACCGTTACTTTGCCGGTCTTCCCTGGTTCAGTGCCAATTCCTCTTTCCGTTTATTGTGc                                 <  2:444295/239‑1 (MQ=255)
                                                                                                                                                                                                                        tGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTTGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCGCTTACCGTTTATTGTGCGAAAAGCCaaaa                     <  1:291615/250‑1 (MQ=255)
                                                                                                                                                                                                                         gACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAAcc                   <  1:289585/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGGTGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTg                 >  1:550706/1‑251 (MQ=255)
                                                                                                                                                                                                                             cctttgAGCGGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACCGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAATCCTggg               <  2:29925/247‑1 (MQ=255)
                                                                                                                                                                                                                               cttgAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCGTATTTTGTGAAAGACCAACAACTATAACAGCCAGGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATTATGTGGGTGTGCCTACGGGCGGTTGCGATGTTTCTGCTCCTGATGTTCCCGTTATTCTGCCGGTCTACTCTGCTTCAGTGCCAATTTCTCTCACCGTTTATTGTGCGAAATGCCAAATCCTggg               <  2:680998/247‑1 (MQ=255)
                                                                                                                                                                                                                                  gAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTAtt          <  1:165190/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggcgggGGGGAGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTcc  <  2:515340/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
TCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  NC_000913/4548981‑4549464

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