Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 155,463 A→T *247K (TAA→AAA)  yadV ← putative fimbrial chaperone YadV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,4630AT98.1% 190.6 / ‑6.2 54*247K (TAA→AAA) yadVputative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (27/25);  minor base G (0/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGA  >  NC_000913/155223‑155708
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTAACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGTAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGcc                                                                                                                                                                                                                                             >  2:651357/1‑250 (MQ=255)
aaCACTGACGTCGTAAACAACCGGTATAGCGGGGTTACATTAACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGAGCGCAAAACGAGGCGAAACGGGTAAGATGATGATAATTTTTAGAAAAATCTATAGTCACGCTATGCTTCATGCGGAATTATATCCTTGCCTGATGTA‑AGCCTGGc                                                                                                                                                                                                                                              >  2:543308/1‑250 (MQ=255)
        cgtcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTGGCCTGAT‑TTCAGCCTGGCATTACCTTc                                                                                                                                                                                                                                     <  2:276435/251‑1 (MQ=255)
                               gggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                                                                                                                                                                                                               >  2:363028/1‑250 (MQ=255)
                                                    cTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg                                                                                                                                                                                         >  1:395381/1‑251 (MQ=255)
                                                                 tAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCtt                                                                                                                                                                             <  2:1048/250‑1 (MQ=255)
                                                                  aGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTg                                                                                                                                                                           >  1:307334/1‑251 (MQ=255)
                                                                      gCATTCTGCCCCCTCAGGACGGTCTGGCCATATTCAACGAGTTCAGCACTGAAAGAAGTATTGCAATACAGCAGCGCGCAACACGTGGCGACACGGGTCAGATGATGCTAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCaa                                                                                                                                                                                                                                   >  2:362776/1‑191 (MQ=255)
                                                                      gCATTCTGCCCCAGCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGAGGATAATTTTTAGTATATTCGATAGTCACGCTATGCGTCCTTCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCaa                                                                                                                                                                                                                                   <  1:362776/191‑1 (MQ=255)
                                                                       cATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCGCTTGTGCAGAATTCGCTTTGCCAtt                                                                                                                                                                      <  1:484405/251‑1 (MQ=255)
                                                                        aTTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAATTGCACTTTTGCAGAATTCGCTTTGCCAtt                                                                                                                                                                      >  1:24935/1‑250 (MQ=255)
                                                                          tCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTaa                                                                                                                                                                    <  2:499477/250‑1 (MQ=255)
                                                                             gCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGcc                                                                                                                                                                <  1:47248/251‑1 (MQ=255)
                                                                                ccATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCAtt                                                                                                                                                             <  1:471062/251‑1 (MQ=255)
                                                                                 cATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTg                                                                                                                                                            >  1:691066/1‑251 (MQ=255)
                                                                                                              gTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGAGCAACACGAGGCGATACGGGTCAGATGATGATAATTTTTAGTCTATTCTATAGTCACGCTCTGCTTCCTGCGGAATTATATCCTTGCCTGCT‑TTCAGCCTGGCATTACATTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGc                                                                                                                               >  1:690128/1‑251 (MQ=255)
                                                                                                                        gAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTTGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCAGCGCTAAATGGAGCAATCAtttt                                                                                                                      <  1:598931/250‑1 (MQ=255)
                                                                                                                        gAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTc                                                                                                                     >  1:577103/1‑251 (MQ=255)
                                                                                                                           aGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGTCCTCATCACTAAATGGTGCAATCATTTTcac                                                                                                                   >  1:592272/1‑250 (MQ=255)
                                                                                                                              aGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCa                                                                                                               <  1:377522/251‑1 (MQ=255)
                                                                                                                               gTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTAAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCa                                                                                                               >  1:621042/1‑250 (MQ=255)
                                                                                                                                tATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAAt                                                                                                             >  1:137900/1‑251 (MQ=255)
                                                                                                                                 aTTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATTATGATAATTTTTAGGATATGCTATAGGCTCGCTATGCTTCCTGCTGAATTATATCCTTGCCTGAT‑TTCAGCCGGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAAAATTCGCTTTGCCATTAAGTCCATCGACTTTCATGTCCTCATCACTAAATGTTGCAATCATTGTCACATCAATc                                                                                                            <  2:621042/251‑1 (MQ=255)
                                                                                                                                   tGCAATACAGCAGCGCGCAAAACGTGGCGATACGGTTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGTCCTCATCACTAAATGGTGCAAGCATTTTCACATCAATCg                                                                                                           <  1:543308/250‑1 (MQ=255)
                                                                                                                                   tGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCg                                                                                                           <  1:225399/250‑1 (MQ=255)
                                                                                                                                             cagcGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTa                                                                                                <  1:30050/251‑1 (MQ=255)
                                                                                                                                                     aaaaCGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGc                                                                                        <  1:174267/251‑1 (MQ=255)
                                                                                                                                                       aaCGTGGAGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCtt                                                                                      >  2:265705/1‑251 (MQ=255)
                                                                                                                                                        aCGTGGCGATACGGGTCAGATGATGATAAGTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTc                                                                                     <  2:691066/251‑1 (MQ=255)
                                                                                                                                                             gCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAAt                                                                                <  1:364274/251‑1 (MQ=255)
                                                                                                                                                                   cGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAAGcaccac                                                                          <  1:216456/251‑1 (MQ=255)
                                                                                                                                                                     ggTCAGATGATGATAATTTTTAGGATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTg                                                                        >  2:154601/1‑251 (MQ=255)
                                                                                                                                                                         aGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTa                                                                     >  1:218458/1‑250 (MQ=255)
                                                                                                                                                                              atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAATTGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAagag                                                               <  2:24935/251‑1 (MQ=255)
                                                                                                                                                                                 atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTGTATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAACGAGACg                                                            >  1:369647/1‑251 (MQ=255)
                                                                                                                                                                                         tAGTATATTCTATCGTAGAGGTATGCTTCCTGCGGAATTAGATACTTGCCTGCT‑TTAAACCTGGCGTTACCTTCAAGTGCGCTACCAAAGTCATTACTGGCGTAAAAATGCAATTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACATTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCATGATTCTAAATCACAACTGCTAAAAGAGACgtagtacgg                                                    >  2:656101/1‑248 (MQ=255)
                                                                                                                                                                                         tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAgg                                                    >  1:510984/1‑251 (MQ=255)
                                                                                                                                                                                           gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAgggg                                                  <  1:343107/251‑1 (MQ=255)
                                                                                                                                                                                               atTCTATAGTCACGCTCTGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACATTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTgg                                              >  1:247177/1‑251 (MQ=255)
                                                                                                                                                                                                    ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTgg                                         >  2:556037/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAACGTCATTAATGGCGTAAAACTGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGTCGTAGTAAGGGGTTGGATTGGTc                                       >  2:171095/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATGTTCACATCAATCGGATAGCGTGTACCGCTAGCTGCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGt                                        <  2:307334/250‑1 (MQ=255)
                                                                                                                                                                                                       gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTaa                                                                                                                                      >  1:172687/1‑155 (MQ=255)
                                                                                                                                                                                                       gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTaa                                                                                                                                      <  2:172687/155‑1 (MQ=255)
                                                                                                                                                                                                          ccGCTATGCTTCCTGCGGAATTATATCCTTGCCGGAG‑TGCAGCCTCGAAGGGGCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGGGGGGGATGACCTCATCACTAAGTGGTGCAATCATTGTCACATCAATCGGATAGCGGTTGCCGCTAGCTTCTGGATCAGAAAGGCTCGATGAGACGTAGTAAGGGGTTGGATTGGTCACTc                                   <  1:582350/250‑1 (MQ=255)
                                                                                                                                                                                                          aCGCTATGCTTCCTGCTGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCAGTACCTTCAATTGCGCCACCAAAGGCATTAATGGCGTAAAAATGCACTTTTGCAGAATGCGCTTGGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCGAGCGTCGAAATCACCACTGCTAAAGGAGACGTAGTAAGTCGTTGTATTGGTCACTc                                   <  2:690128/251‑1 (MQ=255)
                                                                                                                                                                                                          aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTc                                   >  2:303860/1‑251 (MQ=255)
                                                                                                                                                                                                             cTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAggggcgcgaaaaaggtaatat                                 >  2:297182/1‑233 (MQ=255)
                                                                                                                                                                                                             cTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGt                                 >  1:110782/1‑250 (MQ=255)
                                                                                                                                                                                                               aTGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTa                                <  2:391281/249‑1 (MQ=255)
                                                                                                                                                                                                                                  atCCTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAAc           >  1:356362/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccTTGCCTGAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTGTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCGGAAcc          <  2:285395/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cTTGCCTGAT‑TTCAGACTGGAATTACCTTCACTTGCGCCACCAAAATCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTACATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGCACCTg        >  2:414606/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gAT‑TTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGa  <  1:46344/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
AACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGAT‑TACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGA  >  NC_000913/155223‑155708

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