Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,200,160 G→A 100% Q132Q (CAG→CAA)  PP_2807 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,200,1600GA93.3% 78.0 / ‑1.6 30Q132Q (CAG→CAA) PP_2807hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (14/14);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CACCCGCTGTTCCTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAGCCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATCCCCGGTAGCCATCAATGGGAT  >  NC_002947/3200017‑3200304
                                                                                                                                               |                                                                                                                                                
cACCCGCTGTTCCTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAGcccc                                                                                                                                              <  8:237297/148‑1 (MQ=255)
cACCCGCTGTTCCTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAAcccc                                                                                                                                              <  2:2612/148‑1 (MQ=255)
  cccGCTGTTCCTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTg                                                                                                                                            >  2:193676/1‑148 (MQ=255)
           ccTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCTCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAGCCCCTGCACCGTGAc                                                                                                                                   <  1:87984/148‑1 (MQ=255)
                        aGGCAGATCCTGCAAACCCCTACCGATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGct                                                                                                                     >  1:150074/1‑149 (MQ=255)
                                    cAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCt                                                                                                                                             <  8:45947/113‑1 (MQ=255)
                                    cAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCt                                                                                                                                             >  7:45947/1‑113 (MQ=255)
                                     aaaCCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCg                                                                                                        <  4:281333/149‑1 (MQ=255)
                                         cccTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATTCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAg                                                                                                     >  5:24360/1‑148 (MQ=255)
                                         cccTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCGCTTCGCTCTGGCGGCACCCGCAg                                                                                                     >  8:218692/1‑148 (MQ=255)
                                              cccATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATTCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATg                                                                                                 <  6:24360/147‑1 (MQ=255)
                                                                 cGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAg                                                                                    >  6:204525/1‑141 (MQ=255)
                                                                 cGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAg                                                                                    <  5:204525/141‑1 (MQ=255)
                                                                  gACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAat                                                                           >  5:125553/1‑149 (MQ=255)
                                                                     cGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAatgat                                                                        >  2:74110/1‑149 (MQ=255)
                                                                         tCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGACCAGGCTGCAAATGATGCTg                                                                    <  2:150074/149‑1 (MQ=255)
                                                                            agaTGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGcc                                                                 <  6:125553/149‑1 (MQ=255)
                                                                                   gcCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCAttt                                                             >  7:213601/1‑146 (MQ=255)
                                                                                    ccTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCATCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCAttt                                                             >  6:80685/1‑145 (MQ=255)
                                                                                           ttCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCAcc                                                  >  8:5858/1‑149 (MQ=255)
                                                                                              aaCTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAg                                               >  1:121166/1‑149 (MQ=255)
                                                                                                   aGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGaaaa                                          <  5:73597/149‑1 (MQ=255)
                                                                                                                  aTCCAGATTGGCCCGGGCCAAGGCCTCCAACCCCTTCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTTCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGt                           <  5:80685/149‑1 (MQ=255)
                                                                                                                        atgtggccGGGCCAAGGCCTCCAACCCCTGCCCCGGGACGCTGCCACTGCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATccc                     <  8:213601/143‑1 (MQ=255)
                                                                                                                         ttGGCCCGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATcccc                    <  2:121166/149‑1 (MQ=255)
                                                                                                                             cccGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATCTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATCCCCGGTa                <  1:74110/149‑1 (MQ=255)
                                                                                                                               cGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATc                       <  2:272332/140‑1 (MQ=255)
                                                                                                                               cGGGCCAAGGCCTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATc                       >  1:272332/1‑140 (MQ=255)
                                                                                                                                          cTCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATCCCCGGTAGCCATCAATGGGa   <  7:5858/149‑1 (MQ=255)
                                                                                                                                           tCCAACCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATCCCCGGTAGCCATCAATGGGAt  >  5:107234/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CACCCGCTGTTCCTTGGCGCGGCCAGGCAGATCCTGCAAACCCCTACCCATGTCTGGGTGGGCAGCGACCGCATCGAGATGGCGCCTGATATTCAACTGAGCATCACCCAGGCGATCCAGATTGGCCCGGGCCAAGGCCTCCAGCCCCTGCACCGTGACGATGCCACTTCGCTCTGGCGGCACCCGCAGTATGGGCGTGAGGCCAGGCTGCAAATGATGCTGGCCATTTCAGATTTCACCGAGGAAAACGGTGCGACGCGCGTCATCCCCGGTAGCCATCAATGGGAT  >  NC_002947/3200017‑3200304

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