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,0370TC95.2% 131.3 / ‑4.9 43V77A (GTC→GCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/19);  minor base G (2/0);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAAC  >  NC_000913/4548802‑4549278
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
atTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCACAACATGGTCGTGGGTCTTTTGACGCAAAACTTTTTCCATAACGATTATCCGGAAACCATTTAAAACAATGGCACACCGCAACGAgg                                                                                                                                                                                                                                     >  1:270496/1‑251 (MQ=255)
  tGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCCCCAATGTTTATGTAAACCTTGCGCCCGTCGGGGATGTGGGGCAAAACCTGGGCGTGGATCTTTCTTCCCCAACATTTTGTCATTACGATTATTCGGAAACCCTTTCAGACTATGCCACTACGCAACGAGGcc                                                                                                                                                                                                                                   >  2:13136/1‑250 (MQ=255)
   gTAATGAAACGTGTTATGCCCCTGTTTGCTGTCCTGTTGATGGTTTGGTCGATAAATGCCTGTTAATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTc                                                                                                                                                                                                                                  <  2:123546/251‑1 (MQ=255)
   gTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAAACTTGGGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGGTCTTTTGACGCAAATCTTTTTCCCTAAAGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTc                                                                                                                                                                                                                                  >  1:112262/1‑251 (MQ=255)
        gAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAACGCCAATGTTTATGTAAACCCTGCGCCCCTCGTGAATGTGGGGCAAAACCCGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGCAAACATTACAGACTATGCCACACTGCAACGAGGCTCGGCtt                                                                                                                                                                                                                             >  1:144856/1‑251 (MQ=255)
           aCGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAATGTTTATTTAAAACTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGTCGTTGATCTTTTGACGCAAATCTTTTGCCCTAACGATTTTCCGGAAAACCTTTCAGACTATTCCACCCTGCAACGAGGCCCGGCTTAtt                                                                                                                                                                                                                          >  1:377278/1‑250 (MQ=255)
            cGAGTTATTACCCTTTTTGCTGTACTGCTGAGGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATTCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATgg                                                                                                                                                                                                                         <  1:379794/251‑1 (MQ=255)
                 tattaCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAAAGTTTATGTAAACCTTTCGCCCGTCGTGAATGTTGGGCCAAACCTGGGCGTGGATCTTTCGACCCAAATCTTTTTCCCTAAAGATTATCCGGAAACCATTTCAGACTATGCCCCCATTCAACAAGGCTCGGCCTATGgcggcg                                                                                                                                                                                                                    >  1:100580/1‑251 (MQ=255)
                    taCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTTTGTAAACCCTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAATTTTTTCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCCACGAGGCTCGGCTTTTggtggggggg                                                                                                                                                                                                                 >  1:87598/1‑243 (MQ=255)
                             ggcggTACTGCTGAGGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAtt                                                                                                                                                                                                       <  1:42725/247‑1 (MQ=255)
                                   aCTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGTAAACCTTGCCCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGGTCTTTCTACGCCAATCTTTTTTCATAACGATTTTCCGGAAAACATTACAGACTATTGCACACTCCAACGAGGCCCGGCTTATGGCGGCGGGtttttttattttttcc                                                                                                                                                                                                 >  2:118358/1‑237 (MQ=255)
                                           tGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCCGCGCCAAAGTTTTTGTAAACCCTTCGCCCCTCGTGAATGTGGGGGAAAAACCGGTCGGGGGTGTTTTGACGCAAAAATTTTGCCATAACGCGTATCCGGAAAACATTACCAACCATGCCACACCGCAAACAGGCCCGGCTCTTGGCGGCGGGGTTTTTTAATTTTTCCgggccgcg                                                                                                                                                                                          >  2:170278/1‑246 (MQ=255)
                                                            gCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATACCTATTCGCGGCGGCAGCCCCCATATTTATGTAAAACTTGCGCCCGTCGGGCATGTGGGGGAAAAACCGGGCGTGGAACTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCCTTACAGACTATGCCCCACTGCAACCCGGCACGGGTTTTGGCGGCGGGTTATTTAATTTTTCCGGGACCGTAAAATATAGTGGCAgg                                                                                                                                                                          >  1:25718/1‑248 (MQ=255)
                                                                    aTTCGCCTATAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAGATAGTGGCAGTAGCTATCCa                                                                                                                                                                 <  2:457768/250‑1 (MQ=255)
                                                                      gcGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAttt                                                                                                                                                              <  1:76632/250‑1 (MQ=255)
                                                                                     ccAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTTTGTAAAACTTGCGCCCGTCGTGAATGTGGGGCAAAAACTGGTCGTGGGTCTTTCGCCGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCCACGAGGCTCGGCTTTTGGGGGCGGGGTATCTAATTTTTTCGGGGCCGTAAAATATTGTGGGAGGAGCTAttccttttctcccccccacga                                                                                                                                                >  1:421965/1‑230 (MQ=255)
                                                                                       aaTGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAc                                                                                                                                             <  2:140380/251‑1 (MQ=255)
                                                                                                           tGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAt                                                                                                                          <  1:280946/250‑1 (MQ=255)
                                                                                                            ggCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAtt                                                                                                                         <  2:375662/250‑1 (MQ=255)
                                                                                                               ggTGGCAGCGCCAATGTGTATGTAAACCTTGCGCACGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTATCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag                                                                                                                     <  2:244909/251‑1 (MQ=255)
                                                                                                                     aGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCATACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACgg                                                                                                                >  1:384337/1‑250 (MQ=255)
                                                                                                                                    gTAAACCGTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg                                                                                                 <  2:144856/250‑1 (MQ=255)
                                                                                                                                          cTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCt                                                                                          >  1:61081/1‑251 (MQ=255)
                                                                                                                                                            gtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCtgtg                                                                          >  1:150272/1‑249 (MQ=255)
                                                                                                                                                            gtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAg                                                                        <  2:232773/251‑1 (MQ=255)
                                                                                                                                                                           gTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGCGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCggggg                                                         >  2:257368/1‑250 (MQ=255)
                                                                                                                                                                             cGTGGATCTTTAGCCGCAAATCTTTTGCCATAACGATTATCCAGAAACCATTCCAGACTCTGCCACACTGCAACGAAGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTACTTCCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGCTAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGTGCATTGCGGGCggggggg                                                       >  2:223884/1‑248 (MQ=255)
                                                                                                                                                                                ggATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCgggggggggg                                                    >  2:4759/1‑245 (MQ=255)
                                                                                                                                                                                         gACGCAAATCTTTTGCCATAACGATTATCCGGAACCCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCCGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGGGGGGCTTTATTTGACGCCTGTGAGAAGTGCGGGCGGGGTGGCGGTTAAAgatg                                           >  2:373933/1‑248 (MQ=255)
                                                                                                                                                                                         gACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTg                                           <  1:270/251‑1 (MQ=255)
                                                                                                                                                                                           cGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGc                                         <  2:265462/251‑1 (MQ=255)
                                                                                                                                                                                             cAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAAGTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTc                                       <  2:421965/251‑1 (MQ=255)
                                                                                                                                                                                                            aaCGATTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGc                         <  1:165960/250‑1 (MQ=255)
                                                                                                                                                                                                             aCGATTATCCGGAAACCATTACAGACTACGCCACACTGCCACGAGGCTCGGCTTATGGCGGCGTGTTATCTCATTTTTCCGGGACCGCAAACTATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGGGGCCGGTGTCGCTTTTTTTGACCCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCtt                       >  2:33388/1‑251 (MQ=255)
                                                                                                                                                                                                                aTTATCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTAtt                    >  1:181727/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTCCGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGGGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTg                 >  1:80315/1‑251 (MQ=255)
                                                                                                                                                                                                                      cGGAAACCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGAGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCg               >  1:78752/1‑250 (MQ=255)
                                                                                                                                                                                                                         aaaCCATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAg           <  2:371060/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccATTACAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGAcc        <  1:114848/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccATTACAGACTATGACACACTGCAACGAGGCTCGGCTTCTGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGGGGCAGTAGCTATCCATTTCCTACCAACAGCGAAACGCCGCGCGTTGTTTAAAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTGTTTGTAGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAACTGGCTCATTAAATGCCGTGCTTATTTTGAgaccgacc        >  2:261344/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aCAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGGGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaacaa   <  1:91908/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aCAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACcaaca    >  1:313974/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cAGACTATGCCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCaacaac  <  1:174456/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
ATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAAC  >  NC_000913/4548802‑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: