Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 562,573 G→A V80I (GTT→ATT)  sfmH → putative fimbrial adhesin protein SfmH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913562,5730GA98.0% 185.9 / NA 49V80I (GTT→ATT) sfmHputative fimbrial adhesin protein SfmH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/20);  minor base C (1/0);  total (29/20)
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

TATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGACAGGTGATGTTTAC  >  NC_000913/562331‑562817
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
tattattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGa                                                                                                                                                                                                                                               >  1:82036/1‑250 (MQ=255)
    attGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACt                                                                                                                                                                                                                                           >  2:697617/1‑250 (MQ=255)
    attGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACt                                                                                                                                                                                                                                           >  1:517561/1‑250 (MQ=255)
         cGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGAt                                                                                                                                                                                                                                     <  2:542316/251‑1 (MQ=255)
                               aaaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTAc                                                                                                                                                                                                               <  2:82036/251‑1 (MQ=255)
                                    tatCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGAt                                                                                                                                                                                                          >  1:223150/1‑251 (MQ=255)
                                                      gtgCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGAt                                                                                                                                                                                         <  1:655029/250‑1 (MQ=255)
                                                           aaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCt                                                                                                                                                                                   <  1:154095/251‑1 (MQ=255)
                                                            aTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGCGCGATGAGTCTg                                                                                                                                                                                  <  1:483148/251‑1 (MQ=255)
                                                                 aCTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTg                                                                                                                                                                              <  1:697617/250‑1 (MQ=255)
                                                                 aCTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTg                                                                                                                                                                              >  1:377185/1‑250 (MQ=255)
                                                                   tGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                                                                                                            >  1:227767/1‑250 (MQ=255)
                                                                   tGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                                                                                                            >  1:120479/1‑250 (MQ=255)
                                                                                 aTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATAtt                                                                                                                                                             <  1:62678/251‑1 (MQ=255)
                                                                                 aTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGCAAGATTGTGGTTTTACCGGAAAACTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCAACCGGTACGCTGGTGAATTATACATACCGTAGTTACGTCACCAACTTTATTCTTCAGGAACCTATCGATAATTTTAAATATATGCAATTACATAATTATCTATTAGGTGCGTTGAGTCTGGtggttagggcgtgggttatt                                                                                                                                                             >  1:560166/1‑232 (MQ=255)
                                                                                   gAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTc                                                                                                                                                            >  1:694210/1‑250 (MQ=255)
                                                                                                        atCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTatat                                                                                                                                      >  2:396652/1‑251 (MQ=255)
                                                                                                              aaTAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGAt                                                                                                                                <  2:517561/251‑1 (MQ=255)
                                                                                                                               gccagGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACg                                                                                                               >  2:364421/1‑251 (MQ=255)
                                                                                                                               gccagGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACg                                                                                                               >  2:24645/1‑251 (MQ=255)
                                                                                                                               gccaaGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGAAGGAGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGACTTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACg                                                                                                               >  2:249883/1‑251 (MQ=255)
                                                                                                                                  agGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGttt                                                                                                            >  2:518169/1‑251 (MQ=255)
                                                                                                                                      cAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGc                                                                                                         >  1:435844/1‑250 (MQ=255)
                                                                                                                                           tgttgtTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                                                                                                            <  1:306986/178‑1 (MQ=255)
                                                                                                                                           tgttgtTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAAc                                                                                                    <  1:434887/250‑1 (MQ=255)
                                                                                                                                           ngttgttTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                                                                                                            >  2:306986/2‑178 (MQ=255)
                                                                                                                                               gttTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTcc                                                                                               >  1:139878/1‑251 (MQ=255)
                                                                                                                                               gttTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTcc                                                                                               >  2:106657/1‑251 (MQ=255)
                                                                                                                                                 tttACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCAt                                                                                             >  1:13763/1‑251 (MQ=255)
                                                                                                                                                    aCCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTc                                                                                           <  2:13763/250‑1 (MQ=255)
                                                                                                                                                    aCCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCg                                                                                          >  2:666423/1‑251 (MQ=255)
                                                                                                                                                         aaaaaTCCGGCTGGGTAGGGGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTg                                                                                     >  2:294074/1‑251 (MQ=255)
                                                                                                                                                           aaaTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGAt                                                                                   >  1:171715/1‑251 (MQ=255)
                                                                                                                                                                cGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGaa                                                                                                                           <  2:528962/206‑1 (MQ=255)
                                                                                                                                                                cGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGaa                                                                                                                           >  1:528962/1‑206 (MQ=255)
                                                                                                                                                                       gTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGttt                                                                       <  1:364421/251‑1 (MQ=255)
                                                                                                                                                                           gtgtCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGAtt                                                                              <  2:188635/240‑1 (MQ=255)
                                                                                                                                                                           gtgtCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGAtt                                                                              >  1:188635/1‑240 (MQ=255)
                                                                                                                                                                                  aGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGtt                                                            >  2:632756/1‑251 (MQ=255)
                                                                                                                                                                                         tGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAgg                                                      <  2:227767/250‑1 (MQ=255)
                                                                                                                                                                                             cACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCa                                                                                              >  1:454751/1‑206 (MQ=255)
                                                                                                                                                                                             cACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTAtt                                                 >  1:47766/1‑251 (MQ=255)
                                                                                                                                                                                             cACCCGGTACGCTGGTGAATTATACATACCGTAGGTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCa                                                                                              <  2:454751/206‑1 (MQ=255)
                                                                                                                                                                                                  ggTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTcc                                            <  2:435844/251‑1 (MQ=255)
                                                                                                                                                                                                                        acCGTAGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACAtggtgg                           >  1:451170/1‑246 (MQ=255)
                                                                                                                                                                                                                             aGTTATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAg                 <  1:518169/251‑1 (MQ=255)
                                                                                                                                                                                                                                tATGTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGan               <  2:328682/250‑2 (MQ=255)
                                                                                                                                                                                                                                   gTCACCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGACAGGt           >  1:203671/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCCAACTTTATTATTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGGTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGACAGGTGAt        <  1:396652/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ttattaTTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGACAGGTGATGTTTAc  <  2:223150/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
TATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCGGGGTGATGGATTCTCGTTTAATATTTCGTTTAAAGGTTATTCGTCCCTTTATTAACATGGTGGAGATCCCCAGACAGGTGATGTTTAC  >  NC_000913/562331‑562817

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