Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 884,712 G→A G347D (GGT→GAT)  mdfA → multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,7120GA97.6% 128.0 / NA 41G347D (GGT→GAT) mdfAmultidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/23);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑GTACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTA  >  NC_000913/884467‑884948
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
cGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACg                                                                                                                                                                                                                                           <  1:640811/250‑1 (MQ=255)
   ggCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGttt                                                                                                                                                                                                                                       <  2:189218/251‑1 (MQ=255)
     cGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCt                                                                                                                                                                                                                                     <  1:313981/251‑1 (MQ=255)
        tAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGT‑ATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGcc                                                                                                                                                                                                                                  <  1:44660/250‑1 (MQ=255)
                ggTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATggg                                                                                                                                                                                                                          >  2:366799/1‑251 (MQ=255)
                   agcTTGCTGTTAGTGCGTCTGTCCTCGGGCCGCACCGGACTGTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTTGTCTATTGGTCGATGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAATTAGTTATGCTTTCGGTATGGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGTCAGTGATATGAGTAAGG‑ATACGGTGTCTCCCGCGATggttat                                                                                                                                                                                                                       <  2:265266/249‑5 (MQ=255)
                     cTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGANtgc                                                                                                                                                                                                                     >  2:595457/1‑251 (MQ=255)
                             tAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAAt                                                                                                                                                                                                             <  2:160962/250‑1 (MQ=255)
                               gcgcGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCt                                                                                                                                                                                                          <  1:94553/251‑1 (MQ=255)
                                  cgTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGAt                                                                                                                                                                                                       <  1:631286/251‑1 (MQ=255)
                                   gTCTGACCTCGCGACGCACCGAACGATCGCTGATTATAATGAGCGGCCGACCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTCTTTATGGATGACTGCCGGGATAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGAAATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATc                                                                                                                                                                                                      >  2:654135/1‑251 (MQ=255)
                                       gACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCtt                                                                                                                                                                                                    <  1:214284/249‑1 (MQ=255)
                                                             tCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCa                                                                                                                                                                             <  1:59284/250‑1 (MQ=255)
                                                                            ggcggcTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGa                                                                                                                                                              >  1:175575/1‑250 (MQ=255)
                                                                                       gATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACg                                                                                                                                                  >  1:594188/1‑251 (MQ=255)
                                                                                                                   acGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTGTCGGTATTGGTCTGGCGAATTCGGGACTGGTGCTATTAACCCTGTTTGCCAGCGATATGAGTAAAGAGT‑CGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGGTCTTTACCGTTGGTAGTGAAAGCAGCAAACATGCCTTGTTGAACGTGGGCAACGGACTTTTTAATCTCTTCAACCTTGTCa                                                                                                                       <  2:387331/249‑1 (MQ=255)
                                                                                                                           ggTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCaaa                                                                                                                                                                           <  1:201412/190‑1 (MQ=255)
                                                                                                                           ggTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCaaa                                                                                                                                                                           >  2:201412/1‑190 (MQ=255)
                                                                                                                           ggTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAAttt                                                                                                              >  2:539215/1‑251 (MQ=255)
                                                                                                                                   cATCGCACGCGTATTTATGGCTGACTGCCGGGTTAAGTATTTATGATTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTATTCTCTTCAACCTTGTCAACGGAATTTTGTGGCTg                                                                                                      >  2:447161/1‑251 (MQ=255)
                                                                                                                                         aCGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTg                                                                                                >  1:296100/1‑251 (MQ=255)
                                                                                                                                         aCGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTg                                                                                                <  2:175575/251‑1 (MQ=255)
                                                                                                                                           gcgTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGAt                                                                                              <  1:595457/251‑1 (MQ=255)
                                                                                                                                             gTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATgg                                                                                            <  1:24193/251‑1 (MQ=255)
                                                                                                                                             gTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGGCTGGCGCATGCGGGACTGGTGCGATTAACCCTGTTTGACAGCGATATGAGAAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATACTGAGCTTGACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGGCTGTTAAATCTATTCAACCTGGTCAACGGACTTTTGTGGCTGTCGCTGATgg                                                                                            >  1:255247/1‑251 (MQ=255)
                                                                                                                                               aTTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGtt                                                                                          <  2:287952/251‑1 (MQ=255)
                                                                                                                                                     ggNTGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCtt                                                                                     <  2:608405/250‑1 (MQ=255)
                                                                                                                                                               cGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGAt                                                                           >  1:333096/1‑250 (MQ=255)
                                                                                                                                                                 ggTTAAGTATCTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGTTCCGATTAACCCTGTTTGCCAGCTATATGAGTAAAG‑ATACGGTTTCTGCCGCGATCATAATGC‑TACAAATGCTGATCTTTACCTTTGGTATTTAAATCAGCAAACATGCCTGGATGATCGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTCTTAAAAGATaaa                                                                        <  2:612018/251‑1 (MQ=255)
                                                                                                                                                                         aTTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATggg                                                                <  1:247333/251‑1 (MQ=255)
                                                                                                                                                                           ttATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGa                                                               >  1:126899/1‑250 (MQ=255)
                                                                                                                                                                            tATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACATATGGGaa                                                              >  1:240908/1‑250 (MQ=255)
                                                                                                                                                                                 tttCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATtctc                                                        >  1:307233/1‑251 (MQ=255)
                                                                                                                                                                                           ggTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCTATATGAGTAAAG‑ATACGGTTTCGGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGt                                               <  1:348210/250‑1 (MQ=255)
                                                                                                                                                                                              cTGGAGAATGCGGGACTGGTGAGATTAACCCTGTATGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAAAGGGGGAAACGGACTGTTTAATATCTTCAACCTTGTCAACGGGATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTaaa                                            >  2:658611/1‑250 (MQ=255)
                                                                                                                                                                                                   gAATGCGGGACTGGTGCGATTAACCCTGTTTTCCAGCGATATGAGTAAAG‑ATACGGTTTTTGCCGCGATGGGAATGC‑TGCAAATGCTTATCTTTTCCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCATCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATTGGAAATTCTCACGAAGGGTaaaaaaagg                                      <  2:213092/251‑3 (MQ=255)
                                                                                                                                                                                                     aTGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAGCGGACTGTTTAATCTCTTCATCCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGcc                                    <  1:657043/251‑1 (MQ=255)
                                                                                                                                                                                                         gggACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGAc                                >  2:570002/1‑251 (MQ=255)
                                                                                                                                                                                                               ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCt                            <  2:170676/249‑1 (MQ=255)
                                                                                                                                                                                                                                ttNCCAGCGATATGAGTAAAG‑ATGCGTTTTCTGCCGCGGTGGGAATGCGTG‑AAATGCTGATCTTTACCGTTTGTAGTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACTGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGTAAATTCTCTCGAAGGCGAAAAAAATGCCTGACTGCTTTGTGCTATCAGGCATtc          <  2:590234/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gATATGAGTAAAG‑ATACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTa  >  1:9080/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
CGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAG‑GTACGGTTTCTGCCGCGATGGGAATGC‑TGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTA  >  NC_000913/884467‑884948

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