Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,143,420 G→T 19.9% Q189H (CAG→CAT PP_2754 → OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,143,4200GT19.9% 51.6 / 2.7 20Q189H (CAG→CATPP_2754OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base T (2/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

CTTCGACACCGAGCTGCGCTATGGCGACGTCAGCCCGACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGATTACCGGTTCAACGACCA  >  NC_002947/3143285‑3143554
                                                                                                                                       |                                                                                                                                      
cTTCGACACCGAGCTGCGCTATGTCGACGTCAGCCCGACCAGCAAGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGAGCGGCTTCATCATCAACAACAACACGCTGGATGGCCTGATCCTGCAg                                                                                                                                        >  1:425250/1‑136 (MQ=255)
         cgtgcggcgcTAGGGCGACGTCAGCCCGACCTGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAg                                                                                                                               <  2:475322/130‑1 (MQ=255)
             gtGCGCTATGGCGACGTCAGCCCGACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTgg                                                                                                                           <  2:425250/135‑1 (MQ=255)
              tGCGCTATGGCGACGTCAGCCCGACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGa                                                                                                                          >  2:501314/1‑136 (MQ=255)
                             tCAGCCCGACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCATGGCGGCAAGCTGGAGTCGAGCAGTGACCg                                                                                                           >  2:287936/1‑136 (MQ=255)
                                    gACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGcc                                                                                                    <  1:323021/136‑1 (MQ=255)
                                         gCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCg                                                                                                           >  2:514253/1‑124 (MQ=255)
                                         gCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCg                                                                                                           <  1:514253/124‑1 (MQ=255)
                                         gCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCaa                                                                                               <  1:308658/136‑1 (MQ=255)
                                                                            ccGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCATGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGc                                                            >  1:363759/1‑136 (MQ=255)
                                                                                      tGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAggg                                                  <  2:206154/136‑1 (MQ=255)
                                                                                          cgGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGcc                                              >  2:249605/1‑136 (MQ=255)
                                                                                              ttCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACg                                          <  1:249605/136‑1 (MQ=255)
                                                                                                cATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGAt                                        <  2:193914/136‑1 (MQ=255)
                                                                                                          aacaacTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTAcc                              >  1:356974/1‑136 (MQ=255)
                                                                                                             aacTCGCTGGATGGCCTGAGCCTGCAGGGAGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCTTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTTGTGTACCTgg                           <  2:334231/136‑1 (MQ=255)
                                                                                                                     ggATGGCCTGAGCCTGCATGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGAt                   <  2:363759/136‑1 (MQ=255)
                                                                                                                     ggATGGCCTGAGCCTGCATGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGAt                   <  1:287936/136‑1 (MQ=255)
                                                                                                                             tGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGATTACCGGtt           >  2:460445/1‑136 (MQ=255)
                                                                                                                                      aGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGATTACCGGTTCAACGACCa  <  1:446389/136‑1 (MQ=255)
                                                                                                                                      aGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGATTACCGGTTCAACGACCa  <  1:100014/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CTTCGACACCGAGCTGCGCTATGGCGACGTCAGCCCGACCAGCCCGGTGCTGCACGCCAGCGACATTCGCCTGCTGCCGCAGACGCTGCGCGGCTTCATCTTCAACAACAACTCGCTGGATGGCCTGAGCCTGCAGGGCGGCAAGCTGGAGTCGAGCAGTGACCGCAACGCCAGCAACCACCGTGGTGATCTCGGCACCGTCTATGCCGGGCGTATCAAGGGCGCCAACGATGTGGTGTACCTGGGTGGCGATTACCGGTTCAACGACCA  >  NC_002947/3143285‑3143554

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