Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,5940TG32.3% 91.9 / 14.4 68T247P (ACC→CCC) ecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base T (16/28);  new base G (20/1);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 3.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCGTAAATCCCGTACGGTAGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTA  >  NC_000913/308345‑308835
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
aaTCGTAAACCCCGTACGGTAGCCCCCGGGTACCCCCTCCATGATTGCCCATAGTAAAGTTCGGAACGCTTAGTAGCCGCCCATCCCGGGGGGGATGTCCTCCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGgtg                                                                                                                                                                                                                                                    <  2:232083/250‑1 (MQ=255)
 aTCGTAAACCCCGCACGGTTGGCCCGGGTTATCCACTTCATGAGTCCCCATGGTGAATTTCTTAACGCTTAGTGTCCGCCCATGACGGGTGAGATTTTCTTCGCCCGCCGCCGGTAAAAAGGCGATCACTTGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGGTTGGTTTCCCCATGAATGGCCGTAAATCTTCCGTGTTGAAATTTCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGgtgt                                                                                                                                                                                                                                                   <  1:595021/251‑1 (MQ=255)
    gTAAATCCCGTACGGTAGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATAACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGGCATCCGCCCCAAGGGCCGCAAAGTCCAGGGGTCgg                                                                                                                                                                                                                                                >  1:297701/1‑250 (MQ=255)
                 ggTAGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGACCGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCACCGGCCCCTAAGGCTGCACAGTCCCGGGGTCGGGCATTTCACCGGc                                                                                                                                                                                                                                  >  2:49549/1‑251 (MQ=255)
                  gTAGACCCGGGGTATCCACTTCATTATTGCCCATAGTGAAGTTCGGAACGCTTAGTAGCCGCCCATCACGGGTGAGATTTCCTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCa                                                                                                                                                                                                                                 <  2:25074/251‑1 (MQ=255)
                   tAGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGACAGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGCTCACTGGCGTGGATGACTGGCTGCTGTCGAAGATGGTGGAGATGGCCTGGTTTCCCCAGGAAAGGCCGTAAATATTCCCTGCTGAAATGGCGGTCCTCGGCCCCCAAGACTGCAAGGTTCACGGGTCGGGCATTccccccgcc                                                                                                                                                                                                                                  >  1:389980/1‑243 (MQ=255)
                    aGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGGGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGGCATCGGCCCTTAGGACTGCAAGTTCCCGGGGTCGAGCAtttcccccgcca                                                                                                                                                                                                                                 >  1:495425/1‑240 (MQ=255)
                       ccccGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCCAGTTCCAGGGGGCGAGCACTCCACCGGCAAATCg                                                                                                                                                                                                                            >  1:127267/1‑251 (MQ=255)
                          cgggttcccccttcataatgccccAGAGTGAAGTCCTGAACGCTTGGTATCCGCCAATCCCGTGTGAGAGGTACTCCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTg                                                                                                                                                                                                                         <  2:324267/230‑1 (MQ=255)
                                 ccACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCACCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTTACCGGc                                                                                                                                                                                                                  <  1:478610/250‑1 (MQ=255)
                                  cACTTCATGATTGCCCATCGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCCAATCGGGGACCGGcc                                                                                                                                                                                                                 >  1:186467/1‑249 (MQ=255)
                                       cATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGGCATCGGCCCTAAGGACTGCAAGTTCCCGGGGTCGAGCATTCCACCGGCAAAACGGTGACCGGCCAAATcg                                                                                                                                                                                                           >  1:109753/1‑251 (MQ=255)
                                           aTTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCAACGGCCCCAAGGGCTGCAAGTTCCAGGGGTCGAGCATTCCACCCGCAAATCGGTGAACGGGAAAATcgcctt                                                                                                                                                                                                       >  2:482185/1‑248 (MQ=255)
                                                  aTAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGATGAAATGGCGGTCATCGGCCCCAAGGACTGCAAGGTCCCGGGGGTGAGCCTTTCACCCGCAAAACGGGGAACGGCAAAATCGCGGTCATaaac                                                                                                                                                                                                >  1:363127/1‑248 (MQ=255)
                                                      ggAAGTTCTGAACCCTTAGTAGCCGCCAATCGCGGGGGGGAGGTCCTGCCCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAAGCGGTGACCGGCAAAATCGCGTTCATACATCGCt                                                                                                                                                                                            <  2:304056/250‑1 (MQ=255)
                                                       gAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAACAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGCGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGGCAACGGCCCCCAAGACTGCAAGTTCCACGGGTGGCGGCTTCCACCCGCCAATCGGGGGCCCGCCAAATCGCGTTCCTAAATCGCtt                                                                                                                                                                                           >  2:224658/1‑251 (MQ=255)
                                                               gAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGACGCCGGTCAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGACGTAAATCTTCCCTGCTGAAATGGCGGGCATAGGCCCTAAGGACTGCCAGAACCCGGTGTAGAGGATACCCCCGGCAAAACGGTGGCCGGCAAAAACGGGGTTATACATCGCTTTATATAAt                                                                                                                                                                                   >  1:582797/1‑251 (MQ=255)
                                                                aaCGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAAGACTGCAAGGTCCAGGGGTCGAGACTTCCCCCGGCCAATCGGGGACCCGGAAAAACGCGGTCCTACATCGCTTTATATAAtt                                                                                                                                                                                  >  1:319800/1‑251 (MQ=255)
                                                                 aCGCTTAGTAGCCACCCCTCCAGGGTGAGATGCACTTCGCCAGCCGACGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGCCCCTCAGGACTGCCAATTCCCGGGGGCGGGCCTTCACCCCGCAAATCGGCGTCCGGCCCAATCGAGTGCCTACATCGCTTTATATAAttt                                                                                                                                                                                 >  2:279456/1‑250 (MQ=255)
                                                                   ccTTAGGAGCCGCCCAACACGGTTGAGAGGTTCGTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAAGCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGGAAAATCGCGTTCAGACATCGCTTTATATAATTCAc                                                                                                                                                                               <  2:51006/250‑1 (MQ=255)
                                                                              gccCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAACGGCGATCAATGGCGTGGCTGACTCGCTGATGTCGCAGATGGTGGAGCTCGCCCGGCTTCCACAGGAAAGGCCGTAAATCTTCCCTGCTGACATGGCGGTCTTCGGCCCCAATGAGCGCTAGATCCCGGGGTCGTGTCTTTCACCGGCAAATCGGAGACCGGGCAAAACGCGATCACACAACGCTATATATAAtgttcttttctgttg                                                                                                                                                                     >  1:399345/1‑236 (MQ=255)
                                                                                ccATCAGGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGAcc                                                                                                                                                                  <  1:497750/251‑1 (MQ=255)
                                                                                         gggagaggtcttccgcccgccGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGAt                                                                                                                                                         <  2:299744/238‑1 (MQ=255)
                                                                                         gTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGGGTCGAGCATTCCACCCGCAAATCGGGGACCGGCAAAATCGCGTTCCTACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGAt                                                                                                                                                         >  1:472444/1‑251 (MQ=255)
                                                                                                gatccgcgccagccGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGGCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTAc                                                                                                                                                  <  2:511777/245‑1 (MQ=255)
                                                                                                     tCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGc                                                                                                                                               <  1:367092/249‑1 (MQ=255)
                                                                                                       gcccgccGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCCTCCGCCCAAAGGGCTGCAAAGTTCAGGTGTCGGGGATTCCACCCGGCAAACGGGGAACGGGCAAATCCAGGTCATACATCACTTTAAATAAATCACCGTCCTGTTTACCCCTaaccagcccgtaaagagcgg                                                                                                                                           >  1:532767/1‑232 (MQ=255)
                                                                                                        ccctccgccgGTAAAAAGTCGATCTCTGGGGTGGCTGACTGGCTGCTGTCGAAGATGTTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCGGCTGAAATGGCGGTCATCGGCCCTAAGGACTTCAAGTTCCAGGTGTCGAGCATTTCTCCGGCAAATCGTTGACCGGCAAAATCGCGTTCATACATCTCGTTATGTAAGTCACTGGCCTTGTGACCGCGACCTATCCCGTACACCGAGcc                                                                                                                                          <  2:566738/247‑1 (MQ=255)
                                                                                                         ccgccgccgGTAAAAAGGCGATCACTGGCGTGGCTGATTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAAGGCTGCAAGTTCCAGGGGTCGGGCAATCCACCGGCCAATCGGGTACCGGCAAAAACGCGTTCATACATCGCTTTCTATAATTCACTGTCCTTTTGTCCGCTTACGCTCCCGGACAGCggggcg                                                                                                                                         >  1:522073/1‑246 (MQ=255)
                                                                                                                 ggTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCCAAAGGCTGGCAAGTCCCGGGGTCGGGGCTTCCCCCGGCCAAACGGTGACCCGCCAAATCGCGTTCCTACAACCCTTTTTATAATTTACTTTCCCGTTTACCGCTTCCGCTcccccacagggcggccttcggcac                                                                                                                                 >  2:462090/1‑230 (MQ=255)
                                                                                                                    aaaaaGGCGATAACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCCGCCCTAAGGACTGCAAGATTCCGGGGTCGAGCCTTTCACCGGCAAAACGGGGACCGGCAAAATTGCGTTCATACATAGCTTTATATAATTCACTTTCCCGTTGACCGCTCCCGATCCCGGACAGCGCGCCGTCGAgacccac                                                                                                                              >  1:98185/1‑245 (MQ=255)
                                                                                                                     aaaaGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCCTCGGCCCCAAGGGCTGCCAAGTTCCGGTGGCGGGCATTCCCCCGGGCAATCGGGGTCCGGCAAAATCGCGTTCCTACCTCGCTTTTTATAATTTCCTGTCCTGTTTACCCCTACCGATCCCGTaagtacgggggaggggcgggtgg                                                                                                                             >  2:65675/1‑229 (MQ=255)
                                                                                                                     aaaaGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGTCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTa                                                                                                                                                   <  1:65675/229‑1 (MQ=255)
                                                                                                                             aTCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCCAGTTTCAGGTGTCGAGCATTCCACCGGCAAATCGGGGACCGGCAAAATCGCGGTCATACATCGCTTTATATAATTCACTGTCCTGTTGAACGCTACCGATCCCGTAAAGCGAGCCGTCGAGGACCACCTGATGttt                                                                                                                     >  1:64608/1‑250 (MQ=255)
                                                                                                                                   ggCGTGGCTGACTGGCTGCTGGCGAAGATGGTGGATCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCgca                                                                                                                <  2:173174/250‑1 (MQ=255)
                                                                                                                                    gcgtggcttcatGGCTGCTGTCTAAGATGGTGGAGCTGGCCTGGTTTCGCCAAGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGTCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCATAATCGCGTACAGACATCGCTTTAGATAATTCACTGTCCTGTTGACCGGTACCGATCCCGTACAGCGAGCCGTGGAGCACCACATGATGTTCGCgcatt                                                                                                              <  2:278355/240‑3 (MQ=255)
                                                                                                                                            gACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCGGGTTTCCCCAGGAAATGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAAc                                                                                                      <  2:193401/251‑1 (MQ=255)
                                                                                                                                             aCTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACg                                                                                                     <  1:224658/251‑1 (MQ=255)
                                                                                                                                             aCTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCACAGGAACAGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGGCCTAAAGACTCCAAGTTCCCGGCGTCAAGGATTCCACCGGCCAATCGGTGACCGGCCAAATCGCGCTCATACATCCCATTATATAAATCCCTGGCCTGTTGACCGCTACCGATCCCGTAAAGCGAGCCGTCGAGAACCACATGATTGTCGGGCAGTGCAGTaa                                                                                                       >  1:32868/1‑249 (MQ=255)
                                                                                                                                               tGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGGGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGtt                                                                                                   >  1:141303/1‑251 (MQ=255)
                                                                                                                                               tGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGGCATCGGCCCCAAGGACTGCAAGTTCCCGGTGTCGGGCCTTCCACCGGCAAATCGGCGACCGGCAAAATCGCGTTCATACATTGCTTTTTATAATTCACTGTCCTGTTGACCGCTACCGATTCCGTACAGCGGGCCATCGAGCACCCCATGATGTTCGCGCCGTGCAGTAACGtt                                                                                                   >  1:528925/1‑251 (MQ=255)
                                                                                                                                                ggCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGtt                                                                                                   <  1:276441/250‑1 (MQ=255)
                                                                                                                                                ggCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTa                                                                                                  <  2:450241/251‑1 (MQ=255)
                                                                                                                                                ggCTGCTGTCGAAGTGTGGGTAGCGGGCCTCGTTTCCCCAGCAAAGGTCGTAAATCTTCCCGGCTGAAATGGCGGTCATCGGCCCGAAGGACTGCAAGTTCTAGGTGTCTAGGATTCCACCGGCAAATCGGTGACCGGCAAAATCGCTTTCATACATCGCGTTATATAAGTCACTGTCCTGTTGGCCGCTACCGAGCCCGTACTGCGGACCGTCGAGGATGACATCATGGTCGCGCAGTTCAGTAACGt                                                                                                    <  2:371800/249‑1 (MQ=255)
                                                                                                                                                         cGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGAc                                                                                         <  1:420980/251‑1 (MQ=255)
                                                                                                                                                         cGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCAATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGAc                                                                                         <  1:459849/251‑1 (MQ=255)
                                                                                                                                                          gAAGATGGTGGCGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCCTCGGCCCCAAAGGACGCCAAGTCCCGGGGTAGAGCCTTCCACCGGCAAATAGGTGACCGGCAAAATAGCGTTCATCCATCGCTTTTTATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCCTCGAGCACCACCTGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGAc                                                                                         >  2:536173/1‑250 (MQ=255)
                                                                                                                                                          gAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGTGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGATCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACa                                                                                        <  1:556112/251‑1 (MQ=255)
                                                                                                                                                           aaGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGGAACGTTATTCAGCGACAg                                                                                       <  1:155918/251‑1 (MQ=255)
                                                                                                                                                                     gAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTACCTGCTGAAATGGCGGTCATCGGCCCCAAGGACTGAAAGCTTCACGAGTCGGGCATTCCACCCGGAAATCGGTGACCGGCCAAATCGCGTTCATACATCGCTTTTTATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACCGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGGAACGTTATTCAGCGACCGAAAGCTgg                                                                               >  1:241451/1‑250 (MQ=255)
                                                                                                                                                                      aGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGGCATAAAGAATGCAAGTTTCAGGTGTCGCGACTTCCACCGGCCAATCGGGGACCGGCACAATAGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGGGCACCACATGAAGTTCGCGCAGTGCAGTAACGTTATTAAGCGACCGATAGCTggct                                                                             >  1:406083/1‑249 (MQ=255)
                                                                                                                                                                          ggCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGAtgtg                                                                          <  1:527679/250‑1 (MQ=255)
                                                                                                                                                                          ggCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGAtgt                                                                           <  1:239695/249‑1 (MQ=255)
                                                                                                                                                                               gagttCCCCAGGAAAGGCCGGAAATCTTCCCGGCTGATAGGTCGGTGATCGGCCCTAAGGACTGCAAGTTCCAGGGGTCGAGCATTGCACCGGCAAAGCGGGGACCGGCAAAATCGCGTTCATACATCGCTATATATAATGCACTGTCCTGTTGACCGCTACCGAGCCCGTACAGCGAGCCGACGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGtttct                                                                     <  2:26977/247‑4 (MQ=255)
                                                                                                                                                                                  ttCCCCAGGAAAGGCCGTAAATCGTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTGCCAGGTGTCGAGCATTCCACCGGCAAAGCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGGTCGTGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTccc                                                                  <  1:537693/250‑1 (MQ=255)
                                                                                                                                                                                    ccccAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCAGACATCGCTTTATATAATTCACTGTCCTGTGGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTTGATGTGTTGCTCccgcc                                                               <  2:297701/251‑1 (MQ=255)
                                                                                                                                                                                     cccAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCCTCGGCCCCAAGGCCTGCAAGTTCCAGGGGTCGGGCATTCCCCCGGCAAATCGGGGACCGGCAAAATCGCGTTCATACATCCCTTTTTATAATTCCCTGTCCTGTTGACCGCCACCGATCCCGGCCAGCGAGCCGTCGAGCACTCCCTGATTTTCGCGCAGTGCAGGAACGTTATTCCGCGACAGATAGCTGGATGTGTTGCTCccgccg                                                              >  2:338092/1‑251 (MQ=255)
                                                                                                                                                                                                  gTAAATCTTCCCTGCTGAAATGGCGGTCATCGGGCCTAAAGAATACCAGTTCCCGGTGTCGAGCATTCCACCGGGCAATCGGGGAACGGCAAAATCGCGTTCCTACATCCCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTTCCGCGAGCCGTCGAGCACCACCTTATGTTCGCGCAGGGCAGTAACGTTACTTAGCGACAGATAGGTGGTTGAGATGCTCCCGCCGTTACCCAACTggg                                                 >  2:104630/1‑250 (MQ=255)
                                                                                                                                                                                                         ttCCCTGCTGAAATGGCGGTCATCGGCCCCAAGGACTGCAAGGTCCAGGGGTCGAGCATTCCACCGGCAAATCGGGGACCGGCAAAATCGCGCTTATACATCGCTTTATATAATTCAATGTCCTGTTGACCGCTACCGATCCCCTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGCAACGTTATTCAGCGACAGATAGCTGGAGGTGTTGCTCCCGCCGTTACGCAACTGGGTGTtatt                                          >  1:484678/1‑250 (MQ=255)
                                                                                                                                                                                                           cccTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCCCTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTtat                                           >  1:422007/1‑248 (MQ=255)
                                                                                                                                                                                                             ctgctgAAATGGCGGTCATCGGCCCTAAGGACTGCAAGGTCACGGAGTAGAGCATTCCACCGGCAAATCGGTGACCCGCAAAATCGCGTTCATACATCGCTTTTTTTAATTCAATGTCCTGTTGACCGCCACCCATCCCGCACAACGACCCGTCGAGCCCCCCATGATGTTCGCGCAGTGCCGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTTCTCCCGCCGGTCCGCAACTGGTTGTTATAGACg                                      >  1:353049/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaaaTGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGTTGTCGAGCATTCCACCGGCAAATCGGTGACCGACAAAATCGCTTTCATCCATCGCTTTATATAATTCCCTGTCCTGTTGACCGCTACCGATCCCTTACAGCGAGCCGTCGAGCACCACATAATGTTCGCGCATTGCAGTAACGTTATTCAGCGACAGATAGCTGGATTGGTTGCTCCCGCCGTACGCAACTGGTTGTTATAGACGCCCAAg                                <  1:77831/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGAACGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGc                         >  1:132646/1‑251 (MQ=255)
                                                                                                                                                                                                                             cATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTc                       >  1:313698/1‑249 (MQ=255)
                                                                                                                                                                                                                                cGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGCTGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTc                       >  2:45021/1‑247 (MQ=255)
                                                                                                                                                                                                                                cGGCCCTAAGGACAGCGATTTCCAGGTGTCGATCATTCCACCGGCAAATCAGAGACCGGCAACATCGCGTTCAGACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGTTCCCGTACAGCGCGCCGTCGAGCCCCACATGCTGTTCGCGCAGTGCAGTACCGTTATGCAGCGACAGATAGCTGGAGGTGTTGCTCCCGCCGTTACGCACCTGGTTGTGTTAGACGCCCAAGTTATGGCTCCGAt                   >  2:510652/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cccTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTAc                >  2:233798/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   aaGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTccc                                                                  <  2:339452/188‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTccc                                                                  >  1:339452/1‑188 (MQ=255)
                                                                                                                                                                                                                                                       gTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTa  >  1:344292/1‑250 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
CATCGTAAATCCCGTACGGTAGACCCCGGGTATCCACTTCATGATTGCCCATAGTGAAGTTCTGAACGCTTAGTAGCCGCCCATCACGGGTGAGATGTACTTCGCCCGCCGCCGGTAAAAAGGCGATCACTGGCGTGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTA  >  NC_000913/308345‑308835

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