Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,588,476 T→A 7.7% M122L (ATG→TTG)  ydeS ← putative fimbrial protein YdeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,588,4760TA 7.7% 549.5 / 23.1 180M122L (ATG→TTG) ydeSputative fimbrial protein YdeS
Reads supporting (aligned to +/- strand):  ref base T (73/92);  new base A (4/10);  total (78/102)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

AAGGTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCT  >  NC_000913/1588328‑1588625
                                                                                                                                                    |                                                                                                                                                     
aaGGTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCata                                                                                                                                                      <  1:898105/150‑1 (MQ=255)
  ggTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATAtt                                                                                                                                                    >  2:1417098/1‑150 (MQ=255)
   gTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATAttt                                                                                                                                                   <  2:2127476/150‑1 (MQ=255)
    ttGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATg                                                                                                                                                  <  1:237841/150‑1 (MQ=255)
      gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCg                                                                                                                                                <  1:343943/150‑1 (MQ=255)
       ccGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGt                                                                                                                                               >  2:1555989/1‑150 (MQ=255)
           aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                           <  1:316798/150‑3 (MQ=255)
           aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                           <  2:289027/150‑3 (MQ=255)
            gCCCTCACCAGTCCCGGATTGACGGACTTCTTAGTTGACTTCAGACGTGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATACACTGCATCCCGGCATGAAGATCATTCAGTTTCACAGGACTCTAATATGCGTccggt                                                                                                                                          <  2:1693168/150‑4 (MQ=255)
              ccTCACCAGTCCCGGATTGACGGACTTCTGAGTGGACTTCAGACGAGCGGAGTAACGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCAGGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccggtga                                                                                                                                        <  2:1497928/150‑6 (MQ=255)
                 cagcAGTCCCGGATTGAAGGACTTCTGTTTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGGACCAGTGGGATCCACTGGATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     <  1:2328924/147‑1 (MQ=255)
                 cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     <  1:9720/150‑1 (MQ=255)
                 cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     <  1:1711630/150‑1 (MQ=255)
                 cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     <  1:1623646/150‑1 (MQ=255)
                   ccAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCta                                                                                                                                   <  2:392501/150‑1 (MQ=255)
                   ccAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCta                                                                                                                                   <  2:2452737/150‑1 (MQ=255)
                      gTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAc                                                                                                                                <  1:430321/150‑1 (MQ=255)
                      gTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAc                                                                                                                                <  1:1075983/150‑1 (MQ=255)
                       tCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAcc                                                                                                                               >  2:1652888/1‑150 (MQ=255)
                        cccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAccc                                                                                                                              >  1:398116/1‑150 (MQ=255)
                        cccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAccc                                                                                                                              >  1:2161252/1‑150 (MQ=255)
                          cGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCaa                                                                                                                            >  1:708742/1‑150 (MQ=255)
                               tGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCgg                                                                                                                       >  2:2291012/1‑150 (MQ=255)
                                gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                      >  2:508302/1‑150 (MQ=255)
                                gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                      >  2:1611625/1‑150 (MQ=255)
                                 aCGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAg                                                                                                                     >  1:1277685/1‑150 (MQ=255)
                                   ggACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGc                                                                                                                   <  2:1101830/150‑1 (MQ=255)
                                         cTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAACATTTCTATACCCAAACCGGAGGCCGTAtt                                                                                                             <  2:729578/150‑1 (MQ=255)
                                             gTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGctt                                                                                                         >  2:654146/1‑150 (MQ=255)
                                              ttGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttc                                                                                                        <  1:429172/150‑1 (MQ=255)
                                               tGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcc                                                                                                       <  2:1349103/150‑1 (MQ=255)
                                                gACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                      >  2:431355/1‑150 (MQ=255)
                                                gACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                      <  2:2161252/150‑1 (MQ=255)
                                                  cTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                    >  1:2434686/1‑150 (MQ=255)
                                                   ttCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCa                                                                                                   >  2:1876833/1‑150 (MQ=255)
                                                         cGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGt                                                                                             >  1:1029995/1‑150 (MQ=255)
                                                           aGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGttt                                                                                           <  2:453009/150‑1 (MQ=255)
                                                            gCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTc                                                                                          <  1:2235058/150‑1 (MQ=255)
                                                                    aGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGt                                                                                  <  2:398116/150‑1 (MQ=255)
                                                                      gCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtat                                                                                >  2:676202/1‑150 (MQ=255)
                                                                      gCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtat                                                                                >  1:1123907/1‑150 (MQ=255)
                                                                             atTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTc                                                                         <  2:1972738/150‑1 (MQ=255)
                                                                             atTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTc                                                                         <  1:1871384/150‑1 (MQ=255)
                                                                              ttGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCt                                                                        >  1:1038398/1‑150 (MQ=255)
                                                                               tgttctgttcTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTg                                                                       <  2:1998436/150‑1 (MQ=255)
                                                                                gttctgttctGGTACCAGTTGGATCCACTGCATCCCGGCATGAATATCATTCAGTTTCACCGGACGCAAATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGc                                                                      <  2:976246/150‑1 (MQ=255)
                                                                                gttctgttctGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGc                                                                      >  2:1938642/1‑150 (MQ=255)
                                                                                 ttctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCa                                                                     >  2:896083/1‑150 (MQ=255)
                                                                                    tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                  <  2:2365106/150‑1 (MQ=255)
                                                                                    tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                  <  1:12591/150‑1 (MQ=255)
                                                                                    tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                  <  1:2303641/150‑1 (MQ=255)
                                                                                      ttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTc                                                                <  1:460408/150‑1 (MQ=255)
                                                                                       tctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCa                                                               <  1:2217622/150‑1 (MQ=255)
                                                                                       tctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCa                                                               <  2:1779362/150‑1 (MQ=255)
                                                                                        ctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCaa                                                              <  2:2486098/150‑1 (MQ=255)
                                                                                          ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            >  1:798665/1‑150 (MQ=255)
                                                                                          ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            >  1:1958631/1‑150 (MQ=255)
                                                                                          ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            >  1:1831981/1‑150 (MQ=255)
                                                                                           gTACCAGTGGGAAACACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAATGCTACCTTCATCCAGTTTCAACAAAGTATTATTTTCTCCATACTCAAAAc                                                           >  2:478798/1‑150 (MQ=255)
                                                                                             aCCAGTGGGATCCACTGCATCCTGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCg                                                         <  1:1756332/150‑1 (MQ=255)
                                                                                             aCCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCg                                                         <  1:2478432/150‑1 (MQ=255)
                                                                                              ccAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGt                                                        <  1:1555989/150‑1 (MQ=255)
                                                                                                 gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaa                                                     >  1:1365470/1‑150 (MQ=255)
                                                                                                 gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaa                                                     <  2:1227562/150‑1 (MQ=255)
                                                                                                 gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaa                                                     >  2:2097646/1‑150 (MQ=255)
                                                                                                 gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCACATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaa                                                     >  2:381050/1‑150 (MQ=255)
                                                                                                  tGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                    <  1:1417409/150‑1 (MQ=255)
                                                                                                   gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAATTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   <  2:1078407/150‑1 (MQ=255)
                                                                                                   gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   >  2:291773/1‑150 (MQ=255)
                                                                                                   gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   <  1:499331/150‑1 (MQ=255)
                                                                                                   gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   <  1:508302/150‑1 (MQ=255)
                                                                                                    ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATg                                                  <  1:12047/150‑1 (MQ=255)
                                                                                                     gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                 <  1:2175092/150‑1 (MQ=255)
                                                                                                        ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAc                                              >  1:879793/1‑150 (MQ=255)
                                                                                                          aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAccc                                            >  2:2434405/1‑150 (MQ=255)
                                                                                                            tGCATCCCGGAATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACACAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAATTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGa                                          <  2:531302/150‑1 (MQ=255)
                                                                                                              cATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCAGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAc                                        >  1:35570/1‑150 (MQ=255)
                                                                                                               aTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                       <  2:798665/150‑1 (MQ=255)
                                                                                                               aTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                       <  1:1900897/150‑1 (MQ=255)
                                                                                                               aTCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                       <  1:2196873/150‑1 (MQ=255)
                                                                                                                      cATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACGAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                                <  1:2461489/150‑1 (MQ=255)
                                                                                                                      cATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                                <  2:2177750/150‑1 (MQ=255)
                                                                                                                        tGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTc                              <  2:1831981/150‑1 (MQ=255)
                                                                                                                        tGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTc                              <  1:291773/150‑1 (MQ=255)
                                                                                                                         gAAGATCATTCAGTTTCCCCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCACCCCGAACCCCCGTTGTcc                             >  1:1685597/1‑150 (MQ=255)
                                                                                                                         gAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTcc                             >  1:2088879/1‑150 (MQ=255)
                                                                                                                          aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGcccc                            <  1:531267/150‑1 (MQ=255)
                                                                                                                          aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTccc                            <  2:1444795/150‑1 (MQ=255)
                                                                                                                           aGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCt                           >  2:1508981/1‑150 (MQ=255)
                                                                                                                                aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgct                      <  2:594134/150‑1 (MQ=255)
                                                                                                                                 ttCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgctg                     >  1:1574677/1‑150 (MQ=255)
                                                                                                                                  tCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgctgc                    <  2:1277685/150‑1 (MQ=255)
                                                                                                                                     gTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCACCCCGAACCCCCGTTGTCCCTTTGCTGCAtt                 >  2:2467810/1‑150 (MQ=255)
                                                                                                                                     gTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCAtt                 >  1:2431390/1‑150 (MQ=255)
                                                                                                                                       aatattggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCa               >  2:1516646/7‑150 (MQ=255)
                                                                                                                                       aatattggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCa               >  2:2100936/7‑150 (MQ=255)
                                                                                                                                         cACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACt             <  1:587782/150‑1 (MQ=255)
                                                                                                                                          attggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACtt            >  1:1438836/4‑150 (MQ=255)
                                                                                                                                              gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACg        <  2:239283/150‑1 (MQ=255)
                                                                                                                                                    tatTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCt  <  2:691156/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGGTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCT  >  NC_000913/1588328‑1588625

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