Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,221,311 A→G pseudogene (927/1436 nt) mdtQ ← multidrug resistance outer membrane protein MdtQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,221,3110AG97.9% 165.2 / ‑6.0 47pseudogene (927/1436 nt)mdtQmultidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/28);  minor base C (0/1);  total (18/29)
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.99e-01

CCTGCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGACTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTT  >  NZ_CP009273/2221175‑2221444
                                                                                                                                        |                                                                                                                                     
ccTGCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCt                                                                                                                                     <  2:256216/139‑1 (MQ=255)
ccTGCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCt                                                                                                                                     <  1:230799/139‑1 (MQ=255)
ccTGCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCt                                                                                                                                     <  1:542342/139‑1 (MQ=255)
   gCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCaa                                                                                                                                  >  2:250686/1‑139 (MQ=255)
                  tgctgcGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGcgcc                                                                                                                    <  2:132978/138‑1 (MQ=255)
                       cgcgGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTg                                                                                                               >  1:650242/1‑138 (MQ=255)
                       cgcgGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTg                                                                                                               >  1:628739/1‑138 (MQ=255)
                        gcgGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCa                                                                                                             <  1:173928/139‑1 (MQ=255)
                                   gAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCccgc                                                                                                  <  1:250686/139‑1 (MQ=255)
                                       aGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              <  2:761500/139‑1 (MQ=255)
                                       aGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              <  2:360801/139‑1 (MQ=255)
                                        gATCGATTAAGGGCAACGCATCCGGTTGCAGGAAGGCCATCAGGTGGATGTCAGGATAAAATGCCTCTTTTGCCGCATCAATGGTGCTTAGCGTTGCCTGAACGTACCAGTTCGCCGCCTGCAAATCTGCCCGCCGGGc                                                                                             <  2:455748/139‑1 (MQ=255)
                                            gCTTAAGTGCAACGCATCCGGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGt                                                                                         <  2:37873/139‑1 (MQ=255)
                                             cTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTa                                                                                        <  2:313686/139‑1 (MQ=255)
                                             cTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTa                                                                                        <  1:314932/139‑1 (MQ=255)
                                                aaGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAg                                                                                      >  1:259729/1‑138 (MQ=255)
                                                 aGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAgg                                                                                     <  2:303426/138‑1 (MQ=255)
                                                  gTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAg                                                                                   <  1:374624/139‑1 (MQ=255)
                                                   tGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAg                                                                                   >  1:584737/1‑138 (MQ=255)
                                                      aaCGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTAcc                                                                               >  2:157823/1‑139 (MQ=255)
                                                        cGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTAccc                                                                              >  1:210427/1‑138 (MQ=255)
                                                                        aaGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCt                                                             >  2:545777/1‑139 (MQ=255)
                                                                              aTCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTg                                                       <  2:29535/139‑1 (MQ=255)
                                                                              aTCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTg                                                       <  2:259729/139‑1 (MQ=255)
                                                                                     tGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACttt                                                >  1:208637/1‑139 (MQ=255)
                                                                                             ggATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAAcg                                        >  2:757168/1‑139 (MQ=255)
                                                                                             ggATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAAcg                                        >  1:339810/1‑139 (MQ=255)
                                                                                                tAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGa                                     <  2:210427/139‑1 (MQ=255)
                                                                                                  aaaTGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGAc                                    <  2:584737/138‑1 (MQ=255)
                                                                                                  aaaTGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGAc                                    <  2:254920/138‑1 (MQ=255)
                                                                                                  aaaTGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGAc                                    <  1:757168/138‑1 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              <  2:116925/74‑1 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              >  1:116925/1‑74 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGt                                                                                         <  2:299048/79‑1 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGt                                                                                         >  1:299048/1‑79 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTg                                                       <  1:723003/113‑1 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTg                                                       >  2:723003/1‑113 (MQ=255)
                                                                                                        cGCTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGttt                              >  1:184523/1‑138 (MQ=255)
                                                                                                          cTTTTGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAg                           <  2:208637/139‑1 (MQ=255)
                                                                                                             ttGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGt                                                                                         <  1:101348/74‑1 (MQ=255)
                                                                                                             ttGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGt                                                                                         >  2:101348/1‑74 (MQ=255)
                                                                                                             ttGCCGCATCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGc                          <  1:545777/137‑1 (MQ=255)
                                                                                                                     tCAATGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGCTTCAATGACt                <  2:184523/139‑1 (MQ=255)
                                                                                                                         tGGTGCTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGCTTCAATGACTTTGt            <  2:650242/139‑1 (MQ=255)
                                                                                                                              cTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              <  1:145904/52‑1 (MQ=255)
                                                                                                                              cTTAGCGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCggg                                                                                              >  2:145904/1‑52 (MQ=255)
                                                                                                                                   cGATGGCTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTAttt  >  2:22560/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CCTGCCGTAACGCCCATTTGCTGCGCGGAATGACGGAACAGATCGCTTAAGTGCAACGCATCCTGTTGCAGGAAGGCCATCAGGTTGATGTCAGGATAAAATGCCGCTTTTGCCGCATCAATGGTGCTTAGCGATGACTCAACGTACCAGTGCGCCGCCTGCAAATCTGCCCGCCGGGCCAGTAAGGAGTACCCCAGTTCATCAGGAAGCTGGCTTGCCACTTTCGGCAACGCGACCGGTTTAAGCTTCAATGACTTTGTCTGGTTATTT  >  NZ_CP009273/2221175‑2221444

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