Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,590,225 G→A 10.8% G235S (GGC→AGC)  nicP‑I → porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,590,2250GA10.8% 95.5 / 3.5 37G235S (GGC→AGC) nicP‑Iporin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base A (2/2);  total (21/16)
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.43e-01

TGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCA  >  NC_002947/3590077‑3590370
                                                                                                                                                    |                                                                                                                                                 
tGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAgg                                                                                                                                                  <  1:868024/150‑1 (MQ=255)
   cAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCtt                                                                                                                                               >  2:1383985/1‑150 (MQ=255)
    aaGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCttt                                                                                                                                              <  2:1389020/150‑1 (MQ=255)
          aaCGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGc                                                                                                                                        <  2:1042844/150‑1 (MQ=255)
                  gCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACACGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTGGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGca                                                                                                                                >  1:1152569/1‑150 (MQ=255)
                   cTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGcac                                                                                                                               >  1:582687/1‑150 (MQ=255)
                    tCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGcaca                                                                                                                              <  2:90653/150‑1 (MQ=255)
                                cGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCg                                                                                                                  >  2:1107671/1‑150 (MQ=255)
                                      gACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTAGGCGTTCAATGCCAACCAAAGCAGCGTCAGCCGTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACt                                                                                                            <  1:1383985/150‑1 (MQ=255)
                                                aaGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAg                                                                                                  <  1:233020/150‑1 (MQ=255)
                                                          ttCAACTATGTGGGTGGGGACTCTGCGCTCAATGCCAACCAAGGCATCTTCAGCGTTTGGGGCGGTCAGCTGGGGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTg                                                                                        >  2:1231864/1‑150 (MQ=255)
                                                           tCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGTCGACTGGGTACTCAGTGCCAACCTgg                                                                                       <  1:1107671/150‑1 (MQ=255)
                                                                    tgggtgggGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCa                                                                              >  2:525839/1‑150 (MQ=255)
                                                                               tATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCg                                                                                                                  <  1:82909/103‑1 (MQ=255)
                                                                               tATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCg                                                                                                                  >  2:82909/1‑103 (MQ=255)
                                                                                          tGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCAcc                                                                            >  2:931462/1‑130 (MQ=255)
                                                                                          tGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCAcc                                                                            <  1:931462/130‑1 (MQ=255)
                                                                                                ccaaAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCgctgct                                                  >  2:1113358/1‑150 (MQ=255)
                                                                                                                cTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAAc                                  <  1:1113358/150‑1 (MQ=255)
                                                                                                                 tttGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAAcc                                 >  2:1201278/1‑150 (MQ=255)
                                                                                                                   tGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAg                               <  1:525839/150‑1 (MQ=255)
                                                                                                                     gTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGc                             <  2:1386225/150‑1 (MQ=255)
                                                                                                                     gTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGAAAACATGGGTTCTTTAACGGCAAGAGCTCAAGGCTATGCGCGGCTTTGCAACATCGACAAGCAggg                             >  2:1085120/1‑149 (MQ=255)
                                                                                                                       aCGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGt                           <  2:430421/150‑1 (MQ=255)
                                                                                                                         ggTCAGCTGAAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGGTATTTCACCGCACGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGttc                         >  2:663797/1‑150 (MQ=255)
                                                                                                                                tGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGc                  <  2:1152569/150‑1 (MQ=255)
                                                                                                                                tGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGc                  >  2:466043/1‑150 (MQ=255)
                                                                                                                                 ggaggaCATCTACCAACAAGGCTTTGTTGGCCTGAAGCATTGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCt                 <  1:1085120/150‑1 (MQ=255)
                                                                                                                                    ggaCATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGtt              >  1:717969/1‑150 (MQ=255)
                                                                                                                                    ggaCATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGtt              >  2:980488/1‑150 (MQ=255)
                                                                                                                                       ttcaggccAACAAAGCCTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCa            <  1:1180046/143‑1 (MQ=255)
                                                                                                                                       ttcaggccAACAAAGCCTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCa            >  2:1180046/7‑149 (MQ=255)
                                                                                                                                       cATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCAc           >  1:350263/1‑150 (MQ=255)
                                                                                                                                           tACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGcc       >  1:1194166/1‑150 (MQ=255)
                                                                                                                                            gccAACAAAGCCTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTgg                                                <  2:645304/107‑1 (MQ=255)
                                                                                                                                            gccAACAAAGCCTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTgg                                                >  1:645304/2‑108 (MQ=255)
                                                                                                                                            aCCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGccc      >  2:362509/1‑150 (MQ=255)
                                                                                                                                                acaaGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCa  >  1:898092/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCTGGAGGACATCTACCAACAAGGCTTTGTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCA  >  NC_002947/3590077‑3590370

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: