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

AACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGA  >  NC_000913/155223‑155708
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTAACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGTAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGcc                                                                                                                                                                                                                                             >  2:651357/1‑250 (MQ=255)
aaCACTGACGTCGTAAACAACCGGTATAGCGGGGTTACATTAACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGAGCGCAAAACGAGGCGAAACGGGTAAGATGATGATAATTTTTAGAAAAATCTATAGTCACGCTATGCTTCATGCGGAATTATATCCTTGCCTGATGTAAGCCTGGc                                                                                                                                                                                                                                              >  2:543308/1‑250 (MQ=255)
        cgtcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTGGCCTGATTTCAGCCTGGCATTACCTTc                                                                                                                                                                                                                                     <  2:276435/251‑1 (MQ=255)
                               gggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTa                                                                                                                                                                                                               >  2:363028/1‑250 (MQ=255)
                                                    cTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTg                                                                                                                                                                                         >  1:395381/1‑251 (MQ=255)
                                                                 tAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCtt                                                                                                                                                                             <  2:1048/250‑1 (MQ=255)
                                                                  aGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTg                                                                                                                                                                           >  1:307334/1‑251 (MQ=255)
                                                                      gCATTCTGCCCCCTCAGGACGGTCTGGCCATATTCAACGAGTTCAGCACTGAAAGAAGTATTGCAATACAGCAGCGCGCAACACGTGGCGACACGGGTCAGATGATGCTAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                                                                                                                                                                                                   >  2:362776/1‑191 (MQ=255)
                                                                      gCATTCTGCCCCAGCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGAGGATAATTTTTAGTATATTCGATAGTCACGCTATGCGTCCTTCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa                                                                                                                                                                                                                                   <  1:362776/191‑1 (MQ=255)
                                                                       cATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCGCTTGTGCAGAATTCGCTTTGCCAtt                                                                                                                                                                      <  1:484405/251‑1 (MQ=255)
                                                                        aTTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAATTGCACTTTTGCAGAATTCGCTTTGCCAtt                                                                                                                                                                      >  1:24935/1‑250 (MQ=255)
                                                                          tCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTaa                                                                                                                                                                    <  2:499477/250‑1 (MQ=255)
                                                                             gCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGcc                                                                                                                                                                <  1:47248/251‑1 (MQ=255)
                                                                                ccATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCAtt                                                                                                                                                             <  1:471062/251‑1 (MQ=255)
                                                                                 cATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTg                                                                                                                                                            >  1:691066/1‑251 (MQ=255)
                                                                                                              gTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGAGCAACACGAGGCGATACGGGTCAGATGATGATAATTTTTAGTCTATTCTATAGTCACGCTCTGCTTCCTGCGGAATTATATCCTTGCCTGCTTTCAGCCTGGCATTACATTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGc                                                                                                                               >  1:690128/1‑251 (MQ=255)
                                                                                                                        gAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTTGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCAGCGCTAAATGGAGCAATCAtttt                                                                                                                      <  1:598931/250‑1 (MQ=255)
                                                                                                                        gAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTc                                                                                                                     >  1:577103/1‑251 (MQ=255)
                                                                                                                           aGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGTCCTCATCACTAAATGGTGCAATCATTTTcac                                                                                                                   >  1:592272/1‑250 (MQ=255)
                                                                                                                              aGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCa                                                                                                               <  1:377522/251‑1 (MQ=255)
                                                                                                                               gTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTAAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCa                                                                                                               >  1:621042/1‑250 (MQ=255)
                                                                                                                                tATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAAt                                                                                                             >  1:137900/1‑251 (MQ=255)
                                                                                                                                 aTTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATTATGATAATTTTTAGGATATGCTATAGGCTCGCTATGCTTCCTGCTGAATTATATCCTTGCCTGATTTCAGCCGGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAAAATTCGCTTTGCCATTAAGTCCATCGACTTTCATGTCCTCATCACTAAATGTTGCAATCATTGTCACATCAATc                                                                                                            <  2:621042/251‑1 (MQ=255)
                                                                                                                                   tGCAATACAGCAGCGCGCAAAACGTGGCGATACGGTTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGTCCTCATCACTAAATGGTGCAAGCATTTTCACATCAATCg                                                                                                           <  1:543308/250‑1 (MQ=255)
                                                                                                                                   tGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCg                                                                                                           <  1:225399/250‑1 (MQ=255)
                                                                                                                                             cagcGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTa                                                                                                <  1:30050/251‑1 (MQ=255)
                                                                                                                                                     aaaaCGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGc                                                                                        <  1:174267/251‑1 (MQ=255)
                                                                                                                                                       aaCGTGGAGATACGGGTCAGATGCTGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCtt                                                                                      >  2:265705/1‑251 (MQ=255)
                                                                                                                                                        aCGTGGCGATACGGGTCAGATGATGATAAGTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTc                                                                                     <  2:691066/251‑1 (MQ=255)
                                                                                                                                                             gCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAAt                                                                                <  1:364274/251‑1 (MQ=255)
                                                                                                                                                                   cGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAAGcaccac                                                                          <  1:216456/251‑1 (MQ=255)
                                                                                                                                                                     ggTCAGATGATGATAATTTTTAGGATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTg                                                                        >  2:154601/1‑251 (MQ=255)
                                                                                                                                                                         aGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTa                                                                     >  1:218458/1‑250 (MQ=255)
                                                                                                                                                                              atgatAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAATTGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAagag                                                               <  2:24935/251‑1 (MQ=255)
                                                                                                                                                                                 atAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTGTATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAACGAGACg                                                            >  1:369647/1‑251 (MQ=255)
                                                                                                                                                                                         tAGTATATTCTATCGTAGAGGTATGCTTCCTGCGGAATTAGATACTTGCCTGCTTTAAACCTGGCGTTACCTTCAAGTGCGCTACCAAAGTCATTACTGGCGTAAAAATGCAATTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACATTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCATGATTCTAAATCACAACTGCTAAAAGAGACgtagtacgg                                                    >  2:656101/1‑248 (MQ=255)
                                                                                                                                                                                         tAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAgg                                                    >  1:510984/1‑251 (MQ=255)
                                                                                                                                                                                           gTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAgggg                                                  <  1:343107/251‑1 (MQ=255)
                                                                                                                                                                                               atTCTATAGTCACGCTCTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACATTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTgg                                              >  1:247177/1‑251 (MQ=255)
                                                                                                                                                                                                    ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTgg                                         >  2:556037/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAACGTCATTAATGGCGTAAAACTGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGTCGTAGTAAGGGGTTGGATTGGTc                                       >  2:171095/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATGTTCACATCAATCGGATAGCGTGTACCGCTAGCTGCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGt                                        <  2:307334/250‑1 (MQ=255)
                                                                                                                                                                                                       gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTaa                                                                                                                                      >  1:172687/1‑155 (MQ=255)
                                                                                                                                                                                                       gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTaa                                                                                                                                      <  2:172687/155‑1 (MQ=255)
                                                                                                                                                                                                          ccGCTATGCTTCCTGCGGAATTATATCCTTGCCGGAGTGCAGCCTCGAAGGGGCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGGGGGGGATGACCTCATCACTAAGTGGTGCAATCATTGTCACATCAATCGGATAGCGGTTGCCGCTAGCTTCTGGATCAGAAAGGCTCGATGAGACGTAGTAAGGGGTTGGATTGGTCACTc                                   <  1:582350/250‑1 (MQ=255)
                                                                                                                                                                                                          aCGCTATGCTTCCTGCTGAATTATATCCTTGCCTGATTTCAGCCTGGCAGTACCTTCAATTGCGCCACCAAAGGCATTAATGGCGTAAAAATGCACTTTTGCAGAATGCGCTTGGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCGAGCGTCGAAATCACCACTGCTAAAGGAGACGTAGTAAGTCGTTGTATTGGTCACTc                                   <  2:690128/251‑1 (MQ=255)
                                                                                                                                                                                                          aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTc                                   >  2:303860/1‑251 (MQ=255)
                                                                                                                                                                                                             cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAggggcgcgaaaaaggtaatat                                 >  2:297182/1‑233 (MQ=255)
                                                                                                                                                                                                             cTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGt                                 >  1:110782/1‑250 (MQ=255)
                                                                                                                                                                                                               aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTa                                <  2:391281/249‑1 (MQ=255)
                                                                                                                                                                                                                                  atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAAc           >  1:356362/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTGTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCGGAAcc          <  2:285395/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cTTGCCTGATTTCAGACTGGAATTACCTTCACTTGCGCCACCAAAATCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTACATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGCACCTg        >  2:414606/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGa  <  1:46344/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
AACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAGA  >  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: