Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,037 T→C V77A (GTC→GCC)  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,0370TC98.1% 177.4 / ‑5.9 52V77A (GTC→GCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/30);  minor base G (1/0);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAAC  >  NC_000913/4548799‑4549278
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
atgatTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCGGATGGGCTGGTCGGTAAATGCTTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACg                                                                                                                                                                                                                                        <  1:341228/251‑1 (MQ=255)
      gTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGAGGCAGCGCCCATTGTTATGTAAAACCTGAGCCCGTCGTGTATGTGGGGCCAAACCTGGGCGTGGGTCTTTCGACGCCAAACTTTTGCCATTACGAATAT‑CCGGGAAACAATAAAAACTATGCCACCCTGAAACGaagacc                                                                                                                                                                                                                                  >  2:182944/1‑246 (MQ=255)
       tAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTGAGGGGCGGGTCGGTAAATGCCTGTTCATGCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTc                                                                                                                                                                                                                                  <  1:152715/250‑1 (MQ=255)
           gaacctaGTTATTCCCTTGTTTGCTGTACTCGTGAGTGGCTGGTGGGTAATTCCTTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCt                                                                                                                                                                                                                              <  2:310122/244‑1 (MQ=255)
                 aGTTATTCCCCTTTTTGCTTTCCGGGTGAGGGGTTGGTAGGTAAATGCCTGGTCATCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAG‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATggcg                                                                                                                                                                                                                       <  2:95559/251‑1 (MQ=255)
                   ttattcCCCTGTTTCTTGTACGGCTGATGGGTTGGTCGGTAAATGCCTGGTAATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATggcgg                                                                                                                                                                                                                      <  1:276462/250‑1 (MQ=255)
                    tattaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGGCGTGGGTCTTTCGACGCCAATCTTTTGCCATAACGATTAT‑CCGGAAAACATTACAGACTATGCCACACTGCCACGAGGCTCGGCtttttggggcg                                                                                                                                                                                                                    >  1:15337/1‑242 (MQ=255)
                                tgcgttACTGCTGAGGGGCGGGGCGATAAATGCCTGGTCATCCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑t                                                                                                                                                                                                        <  2:118839/245‑1 (MQ=255)
                                        tgctgATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCACAATTTTGTAAACCTTGCGCCCGTCGTGAAATTTGGGCAAAACCTTGTCGTGGATCTTTTTACGCAAATTTTTTTCCATAACAATTTT‑CCCGGAACCATTTCACACTAATCCACACTGCCACACGGCTCGGCTTTTTGCGGCTGGTTATTTAACTTTTTCCgg                                                                                                                                                                                                >  2:305571/1‑251 (MQ=255)
                                        tactgAGGGGCTGGGCGGTAAATGCCTTGTCATCCGCTTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGATTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTGTCCgg                                                                                                                                                                                                <  2:138158/248‑1 (MQ=255)
                                             aTGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGGGCCCATGTTTTTGAAAACCTTGCGCCCCTCGGGGATGTGGGGCAAAAACAGGGCGTGGACCTTTTGACGCAAATCTTTTGCCATAACGATTAT‑TCGGTAACCAATACAGACTATGGCACACCGCAACGAGGCTCGGCTTTTTGCGGCGCGGTATCTAAATTTTTCCGGGCCCg                                                                                                                                                                                           >  1:199056/1‑251 (MQ=255)
                                             aTGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGAGGCAGCCCCCAAGTTTATGTAAACCTTGCGCCCGGCGTGAATGTGGGGCAAAACCTTGTCGGGGATATTTTGTCGCAAAACTTTTTCCATAAAGGTTAT‑CCGGAAAACATTACAGACTATGCCACAATGCAACGAGGACCGGCTGTTGGCGGGGGGTTATTTAA‑tttttttggggcccg                                                                                                                                                                                          >  1:298237/1‑241 (MQ=255)
                                                ggCTGGTCGGTAAGTGCCTGTTCATACGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTaaa                                                                                                                                                                                       <  1:8190/251‑1 (MQ=255)
                                                           aaaTGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAGGGCCACTGTTTATGTAAACATTGCGCCCGTCGGGGATGTGGGGGAAAACCAGGTCGTGGGTCTTTTTGCGCCAAACTTTTTTCCTAACGATTATTCCGGAAAC‑ATTTCAGCCCATGCCACACTGCAACGAGGCTCGGCTTATTGCGGCGCGGTTTCTAT‑TTTTTcctgggccgggaaatattaatgg                                                                                                                                                                             >  2:30766/1‑229 (MQ=255)
                                                                         tcgactgaaaaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTAGGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAtt                                                                                                                                                               <  2:182788/242‑1 (MQ=255)
                                                                                   acccGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTAccacca                                                                                                                                                    <  2:15337/249‑1 (MQ=255)
                                                                                                  cGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGTGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                                                                                                                     <  1:30766/251‑1 (MQ=255)
                                                                                                                      gCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACg                                                                                                                 <  1:319139/251‑1 (MQ=255)
                                                                                                                       cAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAACACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAACTATAGTGGCCGTAGCTATCCCTTTCCTCCCACCAGCGAACCCCCTCGCGTTGTTTATACTTCGAGAACgg                                                                                                                >  2:23247/1‑251 (MQ=255)
                                                                                                                        aGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCGTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGa                                                                                                               <  1:443146/251‑1 (MQ=255)
                                                                                                                           gcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTGTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAGTCTAGAACGGATaa                                                                                                            <  2:417284/251‑1 (MQ=255)
                                                                                                                             cAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACATGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGc                                                                                                          >  1:284602/1‑251 (MQ=255)
                                                                                                                                 gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGt                                                                                                       <  1:288112/250‑1 (MQ=255)
                                                                                                                                   ttATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTgg                                                                                                     <  1:316513/250‑1 (MQ=255)
                                                                                                                                     aTGTAAACCTGGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCTAGAACGGATAAGCCGTGGcct                                                                                                  <  2:433014/251‑2 (MQ=255)
                                                                                                                                                gcgcCCGGCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCttt                                                                                        <  1:23247/250‑1 (MQ=255)
                                                                                                                                                                 ggggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGc                                                                       <  1:305571/250‑1 (MQ=255)
                                                                                                                                                                     cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTgc                                                                  <  1:257376/251‑1 (MQ=255)
                                                                                                                                                                         aCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTgcgggc                                                              <  1:182944/251‑1 (MQ=255)
                                                                                                                                                                          ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGGAACGCCGCGCGTTGTTTATAATTCgag                                                                                                                     >  1:247722/1‑195 (MQ=255)
                                                                                                                                                                          ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                                                                                                                     <  2:247722/195‑1 (MQ=255)
                                                                                                                                                                                cGTGGATCTTTCGACGCAAATCGTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTg                                                        <  2:176305/250‑1 (MQ=255)
                                                                                                                                                                                   ggATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCAACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGGCCGAAAAATATAGTGGCAAGAGCTATCCATTTCCTAACACCAGAGCAAAGCCGCGCGTTGTTTATAATATGAGAACGGATAAGCCGTGGCCGGGGGCGCTTTTTTTGTAGCCTGTGAGCAGGGCGGGCggggggggg                                                     >  2:450341/1‑245 (MQ=255)
                                                                                                                                                                                         ttCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTaa                                                >  1:413916/1‑249 (MQ=255)
                                                                                                                                                                                                 aaagcTTTTGCCATAACGATTAG‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTTTTATCTAA‑TTTTGCCGGGACCGTAAAATATAGTGGCAGTAGCGCTCCAtt                                                                                                                                                               <  2:310863/126‑1 (MQ=255)
                                                                                                                                                                                                 aaaTCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGCCCGTAAAATATAGTGGCAGTAGATATccctt                                                                                                                                                               >  1:310863/1‑127 (MQ=255)
                                                                                                                                                                                                     cTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTa                                   >  1:92243/1‑250 (MQ=255)
                                                                                                                                                                                                       tttGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAAtt                                >  1:460289/1‑251 (MQ=255)
                                                                                                                                                                                                         tGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTg                               <  1:111722/250‑1 (MQ=255)
                                                                                                                                                                                                            cATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGt                           >  1:159230/1‑251 (MQ=255)
                                                                                                                                                                                                            cATAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGt                           >  1:376537/1‑251 (MQ=255)
                                                                                                                                                                                                             aTAACGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTg                          <  1:379989/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAc                                                                                                                                             >  2:325755/1‑132 (MQ=255)
                                                                                                                                                                                                                 cGATTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAc                                                                                                                                             <  1:325755/132‑1 (MQ=255)
                                                                                                                                                                                                                   aTTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCggg                                                                                                      >  2:169537/1‑167 (MQ=255)
                                                                                                                                                                                                                   aTTAT‑CCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTg                                                                                                      <  1:169537/169‑1 (MQ=255)
                                                                                                                                                                                                                     tAT‑CCGGCAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTTTCCATTTCCTACCACCCGCGTAACGCCGCGCGGTGGTTATAATTCGAGAACGGATAAGCCCTGGCCGGGGGCGCTTTATTTGCCGCCCGTGAGCAGGTCGGGCGGGGGGGCGGTTAAAGCTGGCGCCTTAATTTCCGTGATTAtttt                  >  2:49561/1‑251 (MQ=255)
                                                                                                                                                                                                                          cGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCgc              >  1:209036/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGGAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGGTTAAAGCTGGCTCATTAATTGCCGGGCTTATTTTGCGACAGAcc        >  2:136146/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cATTACAGACTATGCCACTCTGCAACGAGGCTCGTCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGTCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTGATTTTGCGACAGACca       <  2:43631/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cATTACAGACTATGCCACACTGCACCGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAACATATAGTGGCAGTAGGTATCCATTTCCTACCAACAGCGAAAAGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTTACGCCTGTGAGCACTGCGGGCGGGGGGGCGATTAAAGATGGACCCTTAATTGCCGTGCTTATTTTTCgaccgacc        >  2:165496/1‑250 (MQ=255)
                                                                                                                                                                                                                                      cAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCaacaac  <  2:284602/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
ATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTAT‑CCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAA‑TTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAAC  >  NC_000913/4548799‑4549278

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