Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,108,466 C→T 8.9% W55* (TGG→TGA)  PP_2721 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,108,4660CT 8.9% 123.5 / 2.4 45W55* (TGG→TGA) PP_2721membrane protein
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base T (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

TCAGCAGGTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATCCAGCTCT  >  NC_002947/3108318‑3108611
                                                                                                                                                    |                                                                                                                                                   
tCAGCAGGTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGc                                                                                                                                                     <  7:117585/149‑1 (MQ=255)
tCAGCAGGTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGc                                                                                                                                                     >  4:161171/1‑149 (MQ=255)
   gcagGTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCAc                                                                                                                                                  <  2:142795/149‑1 (MQ=255)
       gTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCAccccc                                                                                                                                              <  3:161171/149‑1 (MQ=255)
            aaGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCACCTCAGGTTCATAAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCCCCCCCAggcg                                                                                                                                         >  6:400577/1‑146 (MQ=255)
             aGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACt                                                                                                                                         <  4:178620/148‑1 (MQ=255)
              ggCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGa                                                                                                                                       <  4:385404/149‑1 (MQ=255)
                 gTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACca                                                                                                                                    >  2:228459/1‑149 (MQ=255)
                 gTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACca                                                                                                                                    <  5:294946/149‑1 (MQ=255)
                  tatcATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACcac                                                                                                                                   <  8:77350/146‑1 (MQ=255)
                      aTCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTa                                                                                                                               >  2:69592/1‑149 (MQ=255)
                      aTCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTa                                                                                                                               >  3:384381/1‑149 (MQ=255)
                      aTCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTa                                                                                                                               >  7:91181/1‑149 (MQ=255)
                             ccgccCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCt                                                                                                                        <  2:54118/149‑1 (MQ=255)
                             ccgccCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCt                                                                                                                        <  1:349652/149‑1 (MQ=255)
                                             gagAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAAcagca                                                                                                         >  7:97397/3‑148 (MQ=255)
                                                  ggCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGcc                                                                                                   >  7:195879/1‑149 (MQ=255)
                                                   gcgcTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGccc                                                                                                  >  4:358413/1‑149 (MQ=255)
                                                    cgcTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCcca                                                                                                 <  8:195879/149‑1 (MQ=255)
                                                         tCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CAt                                                                                           <  8:91181/148‑1 (MQ=255)
                                                         tCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CAt                                                                                           <  1:69592/148‑1 (MQ=255)
                                                           cTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATgcg                                                                                        <  5:400577/149‑1 (MQ=255)
                                                              gAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAAcagca                                                                                                         <  6:204624/131‑1 (MQ=255)
                                                                   gcgTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACATTCAAAGCCCAGGCCTGGGCATGCGACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTg                                                                                <  6:110281/149‑1 (MQ=255)
                                                                          ttCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTg                                                                         >  5:143294/1‑149 (MQ=255)
                                                                          ttCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTg                                                                         <  3:358413/149‑1 (MQ=255)
                                                                               gcgcCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGAc                                                                     <  4:384381/148‑1 (MQ=255)
                                                                                  ccACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACGCCCAGACTGCCCACAGGATGTTGCTGGCGAGAAAGCAGCAGCAGCCCGC‑‑CATGCGGCCTCCTGGCCGGTGAGACCCGa                                                                 >  8:332101/1‑149 (MQ=255)
                                                                                            aGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGc                                                       >  2:322898/1‑149 (MQ=255)
                                                                                                   tcaGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGccac                                                 >  7:330281/1‑148 (MQ=255)
                                                                                                      ggCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCaccacc                                             >  4:372904/1‑149 (MQ=255)
                                                                                                             cacaaGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCCC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCCGGCTGCCACCACCGTGACCa                                      >  5:203835/1‑149 (MQ=255)
                                                                                                                aaGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGc                                                       >  5:290581/1‑129 (MQ=255)
                                                                                                                aaGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGc                                                       <  6:290581/129‑1 (MQ=255)
                                                                                                                  gCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATg                                 >  7:354231/1‑149 (MQ=255)
                                                                                                                       agcaCAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATggccgg                            <  1:322898/149‑1 (MQ=255)
                                                                                                                                aaaGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCa                    <  6:143294/148‑1 (MQ=255)
                                                                                                                                    cccAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCACGGGAGTCCCCGTACAGGACGCTGAGACCCGATAAACCACGCTGCCACCACCCAGACCAGCAAGTCCGGCCACTCCAGc                  >  1:319559/1‑148 (MQ=255)
                                                                                                                                     ccAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGt              >  4:176852/1‑149 (MQ=255)
                                                                                                                                          ccTGGGCATGCCACCCCCAGACTGACCACTGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATc         <  8:97397/149‑1 (MQ=255)
                                                                                                                                          ccTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATc         <  1:228459/149‑1 (MQ=255)
                                                                                                                                              ggCATGTCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGt                                           >  6:387496/1‑111 (MQ=255)
                                                                                                                                              ggCATGTCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGt                                           >  6:382215/1‑111 (MQ=255)
                                                                                                                                              ggCATGTCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGt                                           <  5:387496/111‑1 (MQ=255)
                                                                                                                                              ggCATGTCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGt                                           <  5:382215/111‑1 (MQ=255)
                                                                                                                                              ggCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATCCAGc     >  6:54582/1‑149 (MQ=255)
                                                                                                                                                 aTGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATCCAGctct  >  4:69986/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCAGCAGGTGATAAGGCGTAGCATCAGCGCCGCCCTGAACCTTCTCTGAAGGCGCTTTCCTTGAAGCGCGTTTATTCTTGCGCCACCCCCTCAGGTTCATCAGGCACAACACAAGCTGCAGCACAATCAAAGCCCAGGCCTGGGCATGCCACCCCCAGACTGACCACAGTATGTTGCTGGCGATAAAACAGCAGAAGCCCAC‑‑CATGCGCCGTCCTGGACGGTGCGACCCGATACACCAGGCTGCCACCACCGTGACCAGCATGGCCGGCCACTCCAGCAAGTCGATCCAGCTCT  >  NC_002947/3108318‑3108611

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