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

TCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  NC_000913/4548981‑4549464
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
actttcgaagataaTCTTGTGTCATAAAGAGTATGCGTAAACAATAACAGACCATGTCACACTGCAACGAGGCTCGGCTTATGGGGGCGTGTTATCTAATGTTTCCGGGACCGTAAAATATAGTGGTAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTGATAATGCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTGTACGCCTGTGAGCAGTTCGGGCGGGGGGGCGATTa                                                                                                                                                                                                                                           <  2:216974/239‑1 (MQ=255)
    tcgaatcaaagtttttGCCATAACGATTATTCTGAATCCATTACATACTATGTCACACTGCAACGAGGCTCGTCTTATGGCGGCGTGTGATCTAATTTTTCCGGGTCCGTAAAATATAGTTGCAGTAGCTATCCATTTCCTACCTCCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGATCAGTGCGGTCGTGGGGGCGATTAAAg                                                                                                                                                                                                                                        <  2:631737/238‑1 (MQ=255)
         gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGAGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCt                                                                                                                                                                                                                                  >  1:367002/1‑251 (MQ=255)
                 tttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCAtt                                                                                                                                                                                                                              <  1:85511/247‑1 (MQ=255)
                     ccATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCg                                                                                                                                                                                                                      >  1:281009/1‑251 (MQ=255)
                            gtatttCCGAAAACCTTTACAGACAAGGTCACACTGCAAATATGTTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCTGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAt                                                                                                                                                                                                               <  2:372456/246‑1 (MQ=255)
                             aTTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtt                                                                                                                                                                                                              >  1:230457/1‑251 (MQ=255)
                                aTCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGAGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCCGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGCTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGACTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGTTGGCCCATTAGTTGCCGTGCTTCttttt                                                                                                                                                                                                           >  1:328512/1‑250 (MQ=255)
                                    ggAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCAACAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAAAGGATAAGCCGGGG‑CCGGTGGCGCTTTATTTGACGCCTGTGGGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCggc                                                                                                                                                                                                       >  1:439893/1‑249 (MQ=255)
                                     gAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACa                                                                                                                                                                                                      <  1:249385/251‑1 (MQ=255)
                                           aTTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaa                                                                                                                                                                                                <  1:548281/251‑1 (MQ=255)
                                              ccAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGAGGCGTGTTATCTAATTTTTCCGGGCCCGGAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGGTTAAAATTCGAGAACGGATAAGCCGTGG‑CCAGGGGCGCTTTTTTTGACGCCTGTGAGCAGTGAGGGCGGGGGGGGGGTTAAAGCTGGCTCATTAATTGACACGCTTATTTTggggacgaccaaaa                                                                                                                                                                                              >  2:433895/2‑237 (MQ=255)
                                                       gTCACACTGCAACGAGGCTCGGCGTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCGCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACa                                                                                                                                                                                     <  2:230457/250‑1 (MQ=255)
                                                            acTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCgatg                                                                                                                                                                               >  1:121416/1‑251 (MQ=255)
                                                                  aCGAGGCTCGGCTTATTGCGGTGTGTTTTCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCCCCCTCGAAACGCCGTGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACATCGATGATTTcc                                                                                                                                                                         <  2:597535/251‑1 (MQ=255)
                                                                   cGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGACCGGATAAGCCGTGG‑CCGGTGGCGCTTTTTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGGAAGAACAACAACTATAACAACGATGATTTccc                                                                                                                                                                        >  1:391829/1‑250 (MQ=255)
                                                                     aGGCTCGGCTTATGGCGGCTTGTTATCTAATTTTTCCGGGACAGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAACCGAGG‑CCGGTGGCGCTTTATTTGACGCCTGGGAGCAGTGCGGGCGGGGGGGCGAGTAAAGCTGGCTCCATAATTGACGTGCTTTTTTTGTGAGAGACCCACAACtcacaaagcgatgctttcc                                                                                                                                                                         >  2:138165/1‑230 (MQ=255)
                                                                             cTTATGGCGGCGAGTTATCTAATTTTTCCGGGACCCTAAAATATAGTGGCAGTAGCTATCCCTTTACTTCCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACCGCTAAGCCGTGGACCGG‑GGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCAATAATTGCCGTGCGTATTTTTCGACCGACCAACAACTATAAAAAGGATGATTTCCACTTTGTGTgg                                                                                                                                                              >  2:9538/1‑251 (MQ=255)
                                                                               tATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATTGTGGCAGTAGCTATCCATTTCCTACCACCAGCGCAACGCCGCGCGTTGGTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGGGCTTTATTTGGAGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTAAAGGTTGTTCATTAATTTCCGTGTTTATTTTTCGACAGACCAACAACTATAACAGCgattgttttcaggttttggtgt                                                                                                                                                             >  2:41266/1‑231 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGGTTAAAGCTGGATCATTAATTGCCGGGCTTATTTTGCGACAGACca                                                                                                                                                                                                 >  2:522169/1‑190 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGTTCATTAATTGCCGTGCTTATTTTGCGACAGACca                                                                                                                                                                                                 <  1:522169/190‑1 (MQ=255)
                                                                                                       tCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATgtgg                                                                                                                                     <  1:83300/250‑1 (MQ=255)
                                                                                                        ccGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGCTgtggtg                                                                                                                                   >  2:165190/1‑251 (MQ=255)
                                                                                                                  aaaaTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGG‑CCGGTGTCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTg                                                                                                                         <  1:288106/251‑1 (MQ=255)
                                                                                                                             gCAGTAGCTATCCATTTCCTACCACCATCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGAt                                                                                                              <  2:455015/251‑1 (MQ=255)
                                                                                                                                  aGCTATCCATTTCCTCCCACCAGCGAAACGCCGAGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACCGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTggcggcatcggtgtttt                                                                                                         >  2:350375/1‑240 (MQ=255)
                                                                                                                                               ctcccaacaGCGAAACGCCGCGCGTTGTTTATCATTCGAGAACGGCTAAGCCGTGG‑CCGGTAGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACCACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATTATGGTGTGGTGCATACTGGCGGCTGCGATGTTTCTgtgcgttatttt                                                                                            >  2:26163/8‑240 (MQ=255)
                                                                                                                                               cTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATggc                                                                                            >  1:558058/1‑249 (MQ=255)
                                                                                                                                                  ccaccaGCGAAACGCCGCGGGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCAcc                                                                                         <  1:529044/251‑1 (MQ=255)
                                                                                                                                                     ccaGCGAAACGCCGCGCGTTGTTTATAATTCTAGAACGTATAAGCCGTGG‑TCGGTGGCTCTTAATTTGACGCCTGTGATCAGTGCGGGAGTGTGGGAGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTTCGACAGACCAACAACTATAACAGCGATGATGTCCAGTTTGTGTGGAATATTTACTCCAATAATGATGTGGTGGTGCCTACTGGCGGCTTCGATGTTTCTGCTCGGGATGTTACCGtt                                                                                      <  2:73793/251‑1 (MQ=255)
                                                                                                                                                     caaGCGAAACGCCGCGCGTATTTTATAATACTGGAACGGATAAGCCGTGG‑TCGGTGGTGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGTGCTATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTGCCAGTTTGTGTGGAATATTTACGCCAATAATGTTGTGGTGGTGCCTTCTGGCGGCTGCGATGTTTCTGCGCGttagggaaccgtt                                                                                      <  2:335824/249‑13 (MQ=255)
                                                                                                                                                       aGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTACTTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGGCTACTGGCGGCTGCGATGGTTCTGCTCGGGATGGCACCGttcc                                                                                    >  1:505889/1‑249 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTCTTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAGAGCTGGCTCATTAATTTCCGTGCTTATTTTGGGACCGACCAACAACTATAACAGCgatgc                                                                                                                                                                              >  2:132285/1‑139 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCgatga                                                                                                                                                                              <  1:132285/140‑1 (MQ=255)
                                                                                                                                                                            tataATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGtt                                                                <  1:9850/250‑1 (MQ=255)
                                                                                                                                                                                  tcgtgaACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGGGCGGGCGGGGGGGCGATTAACGCTGGCTCATTAAGTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTTGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGGTACTCTGCCGGACTACCCTGGTTCAGTGc                                                          <  2:30996/246‑1 (MQ=255)
                                                                                                                                                                                          ggATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCctc                                                 <  1:131184/251‑1 (MQ=255)
                                                                                                                                                                                          ggATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCctc                                                 >  1:515340/1‑251 (MQ=255)
                                                                                                                                                                                                 ccGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGt                                          >  1:29925/1‑251 (MQ=255)
                                                                                                                                                                                                      g‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTAtt                                     <  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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
TCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGG‑CCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCC  >  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: ‑