Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,318,420 T→A 11.1% L613Q (CTG→CAG)  PP_1150 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,318,4200TA11.1% 93.7 / 3.7 36L613Q (CTG→CAG) PP_1150membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/15);  new base A (2/2);  total (19/17)
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.73e-01

CTTCACGGGCTGGATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAGGAGGACACCCATGGG  >  NC_002947/1318272‑1318562
                                                                                                                                                    |                                                                                                                                              
cTTCACGGGCTGGATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAgct                                                                                                                                                <  1:75511/149‑1 (MQ=255)
    aCGGGCTGGATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGctgctg                                                                                                                                            <  8:171559/149‑1 (MQ=255)
        gCTGGATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGTCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGAtg                                                                                                                                        >  6:159272/1‑149 (MQ=255)
            gATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGAtgctgc                                                                                                                                    <  8:113830/149‑1 (MQ=255)
                tCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGt                                                                                                                                <  6:139486/149‑1 (MQ=255)
                tCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGt                                                                                                                                >  3:65252/1‑149 (MQ=255)
                              gTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGc                                                                                                                  >  8:83266/1‑149 (MQ=255)
                                 cACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCCGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCgg                                                                                                                <  7:118259/148‑1 (MQ=255)
                                      gcgcATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGt                                                                                                                         <  3:107861/134‑1 (MQ=255)
                                      gcgcATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGt                                                                                                                         >  4:107861/1‑134 (MQ=255)
                                          aTGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCg                                                                                                      >  3:250873/1‑149 (MQ=255)
                                               ccTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTAcc                                                                                                >  2:92846/1‑149 (MQ=255)
                                                 tGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCAGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTAccc                                                                                               >  7:121403/1‑149 (MQ=255)
                                                 tGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCAGCCGGGCGTGGGCAGCAGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTAccc                                                                                               >  7:124894/1‑149 (MQ=255)
                                                              ggCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGtt                                                                                                                        >  6:245870/1‑111 (MQ=255)
                                                              ggCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGtt                                                                                                                        <  5:245870/111‑1 (MQ=255)
                                                              ggCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCAGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGc                                                                                  <  8:121403/149‑1 (MQ=255)
                                                              ggCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCAGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGc                                                                                  <  8:124894/149‑1 (MQ=255)
                                                               gCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCt                                                                                 >  1:163142/1‑149 (MQ=255)
                                                                 cccTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTc                                                                               <  1:175499/149‑1 (MQ=255)
                                                                       ggggCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTgg                                                                         >  6:95341/1‑149 (MQ=255)
                                                                             gcctgAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAgg                                                                   <  5:64314/145‑1 (MQ=255)
                                                                               ccGAGCGCGAGTACCTGCAACACCTGCAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGCATGCTGCCACTGAGCAGTTCATTGCGGCCTTTCGTCGGCTACCCGCCAGCGACCCCCATCGCCTGGCCGacgag                                                                 >  8:75917/1‑145 (MQ=255)
                                                                                cGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAggggg                                                                <  5:159272/149‑1 (MQ=255)
                                                                                  aGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTc                                                                               >  5:94714/1‑132 (MQ=255)
                                                                                  aGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTc                                                                               <  6:94714/132‑1 (MQ=255)
                                                                                  aGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCa                                                              <  2:163142/149‑1 (MQ=255)
                                                                                            ccTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCaa                                                     <  4:65252/148‑1 (MQ=255)
                                                                                             cTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTAc                                                                                                 <  4:60564/103‑1 (MQ=255)
                                                                                             cTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTAc                                                                                                 >  3:60564/1‑103 (MQ=255)
                                                                                                       tGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGc                                         >  5:116826/1‑149 (MQ=255)
                                                                                                        ggAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTTCTGCAATTGCAGAAAAGCt                                        >  7:72935/1‑149 (MQ=255)
                                                                                                            gCAGGGCTGGCCAATGGGCCTGCTGCCGGGCGTGGGCAGCTGCTTGATGCTGCGAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTggg                                     <  8:234479/148‑1 (MQ=255)
                                                                                                                       cAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGcc                         >  8:71998/1‑149 (MQ=255)
                                                                                                                                     ccGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAggagga            >  2:127804/1‑148 (MQ=255)
                                                                                                                                     ccGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAggagga            <  5:95341/148‑1 (MQ=255)
                                                                                                                                       gggCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAGGAGGacac         >  3:60944/1‑149 (MQ=255)
                                                                                                                                              ggCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAGGAGGACACCCATggg  <  1:127804/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTTCACGGGCTGGATATCGGCGACGAACTGGTGCACCTGCGCATGTGCCTGGCGGTCGCTCAGGCCCCTTTGGGGCGGGCCGAGCGCGAGTACCTGCAACACGTGGAGGCAGCGCTGGCCAAAGGGCCTGCTGCCGGGCGTGGGCAGCTGCTGGATGCTGCCAGTGAGCAGTTCATTGCGGCCTTGCGTCGGCTACCCGCCAGCGACCCGCTTCGCCTGGCCGAGGGGGCAGTACTGCAATTGCAGAAAAGCTGGGGCAAATGGTGCCGCTGGCAGGAGGACACCCATGGG  >  NC_002947/1318272‑1318562

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