Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,471,425 A→G 10.8% pseudogene (907/3495 nt) ydbA → putative outer membrane protein, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,471,4250AG10.8% 93.3 / 3.9 41pseudogene (907/3495 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base A (16/20);  new base G (2/3);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

CACAGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACTTAGTAAA  >  NC_000913/1471286‑1471564
                                                                                                                                           |                                                                                                                                           
cacaGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTATGCAAATAGGGGAATGAGGACAATAATACAGATGATTTAAGATCAATGGCAGATACAGGTGGCGATAAGATTACACGAAGAGGAGAAAATGTaaagt                                                                                                                                          >  2:330762/1‑140 (MQ=255)
 acaGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCt                                                                                                                                         >  1:676670/1‑143 (MQ=255)
     gAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCtaaat                                                                                                                                     >  2:774515/1‑143 (MQ=255)
          aTTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATc                                                                                                                                >  1:665165/1‑143 (MQ=255)
            tGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATctg                                                                                                                              >  2:479844/1‑143 (MQ=255)
            ggtcaaCTATGGTACAATCTGCCCTTTCGGCGCGAGCCCACCGGGCAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATctg                                                                                                                              <  1:374302/139‑1 (MQ=255)
                aaCTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATctgctgc                                                                                                                          >  1:826968/1‑143 (MQ=255)
                      ggTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGACTGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGAc                                                                                                                    <  1:884580/143‑1 (MQ=255)
                      ggTACAAACGGCCCTTTCGTGGTGTGCCAATCGGGGATTGAGTACAATAATACAGATGATTTCGCCTCACTCACCTCTACCGGGGGCGATACCATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGAc                                                                                                                    <  2:513108/143‑1 (MQ=255)
                          cAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAACCGGAGAAACTGTAACGCTACATAAATCTGCTGCTGTGACTGAt                                                                                                                >  1:278579/1‑143 (MQ=255)
                              cTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACCGATAAgc                                                                                                            >  2:551726/1‑143 (MQ=255)
                               tGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgct                                                                                                           <  2:676670/143‑1 (MQ=255)
                                    ttttGGCGTGTTCCAATCGTGGAATGAGTACTATAATATTGATGATTTCACCTCACTGATCTATACCGGTGGCGTTACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgctggctg                                                                                                      <  2:479960/143‑1 (MQ=255)
                                            tgtgCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAAtgtt                                                                                              <  1:359658/143‑1 (MQ=255)
                                             gtgACAATCGGTGAATGAGTACAAGAATGCAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGaaatgtt                                                                                             <  2:171576/143‑6 (MQ=255)
                                                         gAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGaa                                                                                 <  1:774515/143‑1 (MQ=255)
                                                                   aataatACAGATGATTTCACCGCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCgg                                                                       <  1:835220/143‑1 (MQ=255)
                                                                     taataCAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCAGGTg                                                                     >  1:548076/1‑143 (MQ=255)
                                                                        taCAGAGGATGTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATc                                                                  <  1:542118/143‑1 (MQ=255)
                                                                                atTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACg                                                          >  1:417500/1‑143 (MQ=255)
                                                                                atTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACg                                                          <  1:149091/143‑1 (MQ=255)
                                                                                atTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCCGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACg                                                          <  1:695776/143‑1 (MQ=255)
                                                                                    cACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGCTCAAATTACGGTAt                                                      >  2:849169/1‑143 (MQ=255)
                                                                                       cTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTGTCac                                                   >  1:279045/1‑142 (MQ=255)
                                                                                          aCTGATCTATACCGGTGGCGATACGATCACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCg                                                <  2:378862/143‑1 (MQ=255)
                                                                                          aCTGATCTACACCGCTCGCGATACGATTACACGAAGCGGAGAAACTGTCACGCTGAATACACCTCCTTCTTTTACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCg                                                <  1:243700/143‑1 (MQ=255)
                                                                                            tGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGt                                              >  1:53566/1‑143 (MQ=255)
                                                                                                      cGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAGCGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGaaa                                    >  2:224685/1‑143 (MQ=255)
                                                                                                      cGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAGCGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGaaa                                    >  2:673643/1‑143 (MQ=255)
                                                                                                           gCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATAcc                               >  1:718277/1‑143 (MQ=255)
                                                                                                            cGATACGATTACACGAAGCGGAGAAACTGTAGCGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCa                              <  1:224685/143‑1 (MQ=255)
                                                                                                            cGATACGATTACACGAAGCGGAGAAACTGTAGCGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCa                              <  1:673643/143‑1 (MQ=255)
                                                                                                             gATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAGTGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCCTCTGGaaaaccccgg                             >  1:555934/1‑137 (MQ=255)
                                                                                                                  gATTACACGAGGCGGAGAAGCTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGAGCAAATGACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTg                        <  2:173627/143‑1 (MQ=255)
                                                                                                                         cGAAGCGGAGAAACTGTAACGCTAAATAAAGCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGcatcat                 <  1:251080/143‑1 (MQ=255)
                                                                                                                           aaGCGGACAAACTGTAGCGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCAtcatca               <  2:56753/143‑1 (MQ=255)
                                                                                                                             gCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAAt             <  2:548076/143‑1 (MQ=255)
                                                                                                                               ggAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCATCTGGAAAATACCAGCGGTGGCATCATCAAtca           >  2:698779/1‑141 (MQ=255)
                                                                                                                                  gaAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCAGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACtt        >  2:721112/1‑143 (MQ=255)
                                                                                                                                   aaaCTGTAACGCTAACTAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACTTa       <  2:649833/143‑1 (MQ=255)
                                                                                                                                    aaCTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACTTAg      <  1:606851/143‑1 (MQ=255)
                                                                                                                                        gTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACTTAGTaaa  <  2:202914/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CACAGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGGGAATGTTGTTAATAGCGGAACGCTTTCCGGTGATCAAATTACGGTATCAAGCGGTCTTCTGGAAAATACCAGCGGTGGCATCATCAATAACTTAGTAAA  >  NC_000913/1471286‑1471564

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