Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,588,479 T→A 7.2% N121Y (AAT→TAT)  ydeS ← putative fimbrial protein YdeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,588,4790TA 7.2% 558.1 / 22.7 179N121Y (AAT→TAT) ydeSputative fimbrial protein YdeS
Reads supporting (aligned to +/- strand):  ref base T (73/93);  new base A (4/9);  total (77/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.

GGTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAA  >  NC_000913/1588330‑1588628
                                                                                                                                                     |                                                                                                                                                     
 gTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATAttt                                                                                                                                                      <  2:2127476/150‑1 (MQ=255)
  ttGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATg                                                                                                                                                     <  1:237841/150‑1 (MQ=255)
    gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCg                                                                                                                                                   <  1:343943/150‑1 (MQ=255)
       gAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTcc                                                                                                                                                <  1:954215/150‑1 (MQ=255)
         aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                              <  2:289027/150‑3 (MQ=255)
         aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                              <  1:316798/150‑3 (MQ=255)
          gCCCTCACCAGTCCCGGATTGACGGACTTCTTAGTTGACTTCAGACGTGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATACACTGCATCCCGGCATGAAGATCATTCAGTTTCACAGGACTCTAATATGCGTccggt                                                                                                                                             <  2:1693168/150‑4 (MQ=255)
               cagcAGTCCCGGATTGAAGGACTTCTGTTTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGGACCAGTGGGATCCACTGGATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                        <  1:2328924/147‑1 (MQ=255)
               cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                        <  1:1623646/150‑1 (MQ=255)
               cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                        <  1:1711630/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)
                       ccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTGTGGTACCAGTGGGATCCACTGCATCCCGTCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCa                                                                                                                                <  2:284921/150‑1 (MQ=255)
                        cGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCaa                                                                                                                               >  1:708742/1‑150 (MQ=255)
                             tGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCgg                                                                                                                          >  2:2291012/1‑150 (MQ=255)
                              gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                         >  2:1611625/1‑150 (MQ=255)
                              gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                         >  2:508302/1‑150 (MQ=255)
                               aCGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAg                                                                                                                        >  1:1277685/1‑150 (MQ=255)
                                      tCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAt                                                                                                                 <  2:1414617/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)
                                              gACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                         <  2:2161252/150‑1 (MQ=255)
                                              gACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                         >  2:431355/1‑150 (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)
                                                         aGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGttt                                                                                              >  2:1017952/1‑150 (MQ=255)
                                                          gCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTc                                                                                             <  1:2235058/150‑1 (MQ=255)
                                                             gAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaac                                                                                          >  2:1092382/1‑150 (MQ=255)
                                                             gAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCAGGAAGATCATTGAGTTTCACCGGACTCATATTTGCGTCCAATATTTCTATACCCTAACCGGAGGCCGTATTGCGTGCATCATCCAGttggaac                                                                                          <  2:1439286/150‑6 (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)
                                                                              gttctgttctGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGc                                                                         >  2:1121120/1‑150 (MQ=255)
                                                                               ttctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACTGTACGGATTTTTGCGTCCTATATTTCTTTACTCTAATCGGAGTCCTTATTGCTTCTTTCATCCAGTTTGAAGATAGTATTATTTTGTGct                                                                        <  2:1496095/150‑2 (MQ=255)
                                                                               ttctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCa                                                                        >  2:896083/1‑150 (MQ=255)
                                                                                  tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                     <  2:2365106/150‑1 (MQ=255)
                                                                                  tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                     <  1:2303641/150‑1 (MQ=255)
                                                                                  tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                     <  1:12591/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:1831981/1‑150 (MQ=255)
                                                                                        ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                               >  1:1958631/1‑150 (MQ=255)
                                                                                        ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                               >  1:798665/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:2097646/1‑150 (MQ=255)
                                                                                               gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCACATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaa                                                        >  2:381050/1‑150 (MQ=255)
                                                                                                tGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                       <  1:1417409/150‑1 (MQ=255)
                                                                                                 gggATCCTCTGCATCCCGGCATGAAGATCATTCAGTTTAACCGGACGCATATTTGCGTCCATTATTTCTATACCCAAACAGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATGCTCAATACCGTTTAAt                                                      >  2:1500596/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                                                    <  2:1518450/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)
                                                                                                             ctCCCGGCATGAAGATGATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCATACCGGACGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcc                                          <  2:1281083/149‑1 (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)
                                                                                                                              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:2100936/7‑150 (MQ=255)
                                                                                                                                     aatattggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCa                  >  2:1516646/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)
                                                                                                                                                     ttGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGaa  <  2:1556049/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAA  >  NC_000913/1588330‑1588628

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