Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 885,786 T→G L81R (CTG→CGG)  ybjM → inner membrane protein YbjM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273885,7860TG94.7% 60.0 / ‑2.4 19coding (242/378 nt)ybjMinner membrane protein YbjM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (7/11);  minor base . (0/1);  total (7/12)
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.11e-01

TCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCA  >  NZ_CP009273/885678‑885906
                                                                                                            |                                                                                                                        
tCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGc                                                                                             <  1:157945/138‑1 (MQ=255)
      ttcCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCt                                                                                      >  2:294489/1‑139 (MQ=255)
         cTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTgg                                                                                    <  1:175544/138‑1 (MQ=255)
                 gTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACt                                                                          <  2:258496/139‑1 (MQ=255)
                  tCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACt                                                                          <  1:21220/139‑1 (MQ=255)
                  tCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACt                                                                          <  1:341132/139‑1 (MQ=255)
                         cttcttTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGcgcg                                                                   <  1:58867/139‑1 (MQ=255)
                                   caGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGgt                                                         <  2:306222/139‑1 (MQ=255)
                                      cGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTGCTAAGCGCGGTGTTCTGGggtg                                                      >  2:51695/1‑139 (MQ=255)
                                            gaAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTgg                                                <  1:86350/139‑1 (MQ=255)
                                               gTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTggggg                                             <  1:316770/139‑1 (MQ=255)
                                                tCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTggggg                                             >  2:6196/1‑138 (MQ=255)
                                                  cTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCTCTGGGGGCAc                                          <  1:188344/139‑1 (MQ=255)
                                                    gAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACtg                                        >  2:243642/1‑139 (MQ=255)
                                                                   cgcAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGGCGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAgcg                                                                     >  2:195951/1‑95 (MQ=255)
                                                                   cgcAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAgcg                                                                     <  1:195951/95‑1 (MQ=255)
                                                                             gcggcACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGTGCTGGCCTGGTTACTAAGCGCGGTGtt                                                              >  1:189603/1‑92 (MQ=255)
                                                                             gcggcACCCTGTTCGATGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGtt                                                              <  2:189603/92‑1 (MQ=255)
                                                                                             tGCTCATTATGCGGCGGTTTTTTTCACCGACGCGCTCATTCTGGCCAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACcacagca  >  1:60096/1‑136 (MQ=255)
                                                                                                            |                                                                                                                        
TCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGAAACCTCTGTTTGGCGCAATGCTGGCGGCACCCTGTTCGATGCTCATTATGCGGCTGTTTTTTTCACCGACGCGCTCATTCTGGCAAGAGCTGGCATGGTTACTAAGCGCGGTGTTCTGGTGTGCGCTGGGGGCACTGTGTTTCTTATTTATCAGTAGTTTGTTTAAACCACAGCA  >  NZ_CP009273/885678‑885906

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