Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,174,536 G→C 11.4% P136P (CCG→CCC)  PP_1029 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,174,5360GC11.4% 102.4 / 6.1 35P136P (CCG→CCC) PP_1029membrane protein
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base C (2/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GTTGTACGTTGGGCGTGGTGTTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGCTGCAGTACGATAGTGCGTCTGTC  >  NC_002947/1174391‑1174685
                                                                                                                                                 |                                                                                                                                                     
gTTGTACGTTGGGCGTGGTGTTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGGGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATgg                                                                                                                                                   <  1:364211/150‑1 (MQ=255)
 ttGTACGTTGGGCGTGGTGTTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATggg                                                                                                                                                  <  1:159997/150‑1 (MQ=255)
       gTTGGGCGTGGTGTTCCTGCTGCAACTGAAGTCCGACCTGTTGCGGCTGGTGATGGGGGTAGTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGATCCTGTCCGGAGTCTGGGCGGTGCCGATGGGCACGGc                                                                                                                                            >  1:144273/1‑150 (MQ=255)
                ggtgTTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGttt                                                                                                                                   <  2:100487/150‑1 (MQ=255)
                  tgtTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGt                                                                                                                                 >  2:153266/1‑150 (MQ=255)
                                cTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTc                                                                                                                   <  1:228220/150‑1 (MQ=255)
                                              gTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTc                                                                                                                   <  1:180063/136‑1 (MQ=255)
                                              gTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTc                                                                                                                   >  2:180063/1‑136 (MQ=255)
                                                 gctgctGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTGGGGGCTTTGTTCGGCAGCGGCGGGTTTCt                                                                                                  >  2:445760/1‑150 (MQ=255)
                                                         tGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGccc                                                                                          <  1:336748/150‑1 (MQ=255)
                                                                ggTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGTGTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCt                                                                                   <  2:17034/150‑1 (MQ=255)
                                                                    ttCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTGGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAgc                                                                               >  1:408511/1‑150 (MQ=255)
                                                                     tCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCtt                                                                                        <  2:40378/140‑1 (MQ=255)
                                                                     tCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCtt                                                                                        >  1:40378/1‑140 (MQ=255)
                                                                        tGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCAc                                                                           >  1:258125/1‑150 (MQ=255)
                                                                             gCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTg                                                                      <  2:352578/150‑1 (MQ=255)
                                                                               cTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGGGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTggg                                                                    >  1:279130/1‑149 (MQ=255)
                                                                               cTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGa                                                                    <  1:143119/150‑1 (MQ=255)
                                                                                      gttgaCGGGCTGGCAGTAAAAGTGAGGCCAGCGAGCCTGTCCGGCTCCTGGGACGTGCCGATGGGGACGGTAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTATTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGCAg                                                             <  2:144273/146‑1 (MQ=255)
                                                                                                  gCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTAc                                                   >  1:53696/1‑148 (MQ=255)
                                                                                                  gCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTAc                                                   <  2:53696/148‑1 (MQ=255)
                                                                                                        aaaGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCAcc                                           >  1:358710/1‑150 (MQ=255)
                                                                                                            tGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCaga                                       <  2:258125/150‑1 (MQ=255)
                                                                                                                gCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAgcgc                                   >  1:501057/1‑150 (MQ=255)
                                                                                                                      gAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGAt                             <  2:501057/150‑1 (MQ=255)
                                                                                                                               cGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAgcgcgc                                 <  1:45180/137‑1 (MQ=255)
                                                                                                                               cGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAgcgcgc                                 >  2:45180/1‑137 (MQ=255)
                                                                                                                                ggCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGCTGCAGt                   <  1:491228/150‑1 (MQ=255)
                                                                                                                                 gCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGCTGCAGTa                  <  1:153266/150‑1 (MQ=255)
                                                                                                                                        ggCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAg                                                       <  1:493998/106‑1 (MQ=255)
                                                                                                                                        ggCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAg                                                       >  2:493998/1‑106 (MQ=255)
                                                                                                                                         gccgTGCCCATCGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGc                         <  1:559609/132‑1 (MQ=255)
                                                                                                                                         gccgTGCCCATCGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGc                         >  2:559609/4‑135 (MQ=255)
                                                                                                                                           cgTGCCCATCGGCACGGCAGGTGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAAc                                                         >  2:331999/2‑101 (MQ=255)
                                                                                                                                           cgTGCCCATCGGCACGGCAGGTGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAAc                                                         <  1:331999/100‑1 (MQ=255)
                                                                                                                                                 gATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGGTTCTTTATGCCCTTTACCTGAGAGCACGGCTGGAGGTGAAGGAACAGGTACGCGACACCCAGAGCGCGCTGATCAGATGCATTACGATAGTGCGTCTGTc  >  1:262631/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTTGTACGTTGGGCGTGGTGTTCCTGCTGCAACTGAAGTCCGACCTGTTGCTGCTGTTGATGGGGGTATTCGTGACTGCCTATGCACTTTACGGGCTGGCAGTGAAAGTGAGGCCAGCGAGCCTGTCCGGCTTCTGGGCGGTGCCGATGGGCACGGCAGGCGGTTTGTTCGGGGCTTTGTTCGGCAGCGGCGGGTTTCTTTATGCCCTTTACCTGAGCGCACGGCTGGAGGTGAAGGAACAGGTACGCGCCACCCAGAGCGCGCTGATCAGCTGCAGTACGATAGTGCGTCTGTC  >  NC_002947/1174391‑1174685

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