Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,385,157 T→A 15.4% intergenic (‑121/+40) proS ← / ← oprD proline‑‑tRNA ligase/basic amino acid specific porin OprD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,385,1570TA15.4% 60.3 / 5.3 26intergenic (‑121/+40)proS/oprDproline‑‑tRNA ligase/basic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base A (2/2);  total (12/14)
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.46e-01

AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCA  >  NC_002947/1385017‑1385297
                                                                                                                                            |                                                                                                                                            
aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTgccgggct                                                                                                                                     >  2:91244‑M2/1‑142 (MQ=255)
    atTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTgccgggcttttt                                                                                                                                 <  2:266419‑M2/150‑13 (MQ=255)
        aCTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAgg                                                                                                                             >  1:720533/1‑150 (MQ=255)
                          ccGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGcc                                                                                                           >  2:99017/1‑150 (MQ=255)
                                      ccTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa                                                                                               >  1:202152/1‑150 (MQ=255)
                                         aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtcac                                                                                                              >  2:95397‑M2/1‑101 (MQ=255)
                                         aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtcac                                                                                                              <  1:95397‑M2/132‑32 (MQ=255)
                                                  aaaTAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGt                                                                                   <  1:99017/150‑1 (MQ=255)
                                                                     gggggggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAg                                                                >  1:273175/1‑150 (MQ=255)
                                                                          ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgt                                                           >  2:143218/1‑150 (MQ=255)
                                                                           gAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgtc                                                          <  1:266419/150‑1 (MQ=255)
                                                                              aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtc                                                                                                                <  1:346861‑M2/93‑30 (MQ=255)
                                                                              aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtc                                                                                                                >  2:346861‑M2/1‑64 (MQ=255)
                                                                               cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatc                                                      <  1:5931/150‑1 (MQ=255)
                                                                                    gCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATgg                                                 <  2:202152/150‑1 (MQ=255)
                                                                                         cAGGATGCCCCAGGGCGTGACACTTGTCCGCCAGGAAAAAGCCCGGTGTGGTTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCGCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAAt                                            <  2:278251/150‑1 (MQ=255)
                                                                                             aTGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgta                                        <  2:720533/150‑1 (MQ=255)
                                                                                               gCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtagt                                      <  1:143218/150‑1 (MQ=255)
                                                                                                                ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa                                                                 <  1:841612/106‑1 (MQ=255)
                                                                                                                ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa                                                                 >  2:841612/1‑106 (MQ=255)
                                                                                                                            aaaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc         <  2:273175/150‑1 (MQ=255)
                                                                                                                                aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac     >  2:647615/1‑150 (MQ=255)
                                                                                                                                aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac     <  1:647615/150‑1 (MQ=255)
                                                                                                                                   ccGGCGTGATTGCCGGGCTTTTTCACGTGCGCCCGACTAAAACAGAAAGTGATGGGATGCTCAGCGGGTACTGGACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACCACATCGTAAGAACGGTAGGTGGCCTGacgca  >  1:574932/1‑150 (MQ=255)
                                                                                                                                          gATTgccgggctttttcttggctgacaagtgtcacgccctggggcatcctgtacgccttcgtttcccccccgtacaatgcggggctatttcaacgtcttaggcaatcgacacggctaaccatgc                     >  2:393002‑M2/1‑4 (MQ=255)
                                                                                                                                          gATTgccgggctttttcttggctgacaagtgtcacgccctggggcatcctgtacgccttcgtttcccccccgtacaatgcggggctatttcaacgtcttaggcaatcgacacggctaaccatgc                     <  1:393002‑M2/124‑121 (MQ=255)
                                                                                                                                            |                                                                                                                                            
AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCA  >  NC_002947/1385017‑1385297

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