Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,892,539 T→G 7.9% L508R (CTG→CGG)  PP_3435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,892,5390TG 7.9% 142.0 / 2.5 48L508R (CTG→CGG) PP_3435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base G (2/2);  total (27/21)
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.73e-01

AGCTGTCGGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTG  >  NC_002947/3892392‑3892685
                                                                                                                                                   |                                                                                                                                                  
aGCTGTCGGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTg                                                                                                                                                   <  2:266208/149‑1 (MQ=255)
   tGTCGGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGc                                                                                                                                                 >  5:103548/1‑148 (MQ=255)
       gggCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACATGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGCCCCATGCCGGCCGCCCAGTTCCTCGAAGCCATGCCCACTc                                                                                                                                            >  2:239259/1‑149 (MQ=255)
          cAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGcc                                                                                                                                         >  7:94245/1‑149 (MQ=255)
                 gACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGgc                                                                                                                                    <  4:79204/147‑1 (MQ=255)
                    cTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCGGGCGAGTGGCCAGGgcttcga                                                                                                                               <  6:352953/149‑6 (MQ=255)
                       ggCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCa                                                                                                                            <  6:103548/149‑1 (MQ=255)
                       ggCATCGTTGCCGAGGGGGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCGGGCGAGTGGCCAGGgcttcgag                                                                                                                              <  5:371687/147‑7 (MQ=255)
                         cATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGt                                                                                                                          <  1:68665/149‑1 (MQ=255)
                                gCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGc                                                                                                                    >  3:282420/1‑148 (MQ=255)
                                  cGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGc                                                                                                                    >  2:47731/1‑146 (MQ=255)
                                      ggTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATgc                                                                                                             >  6:320894/1‑149 (MQ=255)
                                            gAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGccccc                                                                                                       >  2:1648/1‑149 (MQ=255)
                                            gAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGATCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGccccc                                                                                                       >  1:305247/1‑149 (MQ=255)
                                             aaaCCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCt                                                                                                      <  1:47731/149‑1 (MQ=255)
                                                 cGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCt                                                                                                      >  5:177930/1‑145 (MQ=255)
                                                    ccTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCa                                                                                               <  2:352386/149‑1 (MQ=255)
                                                       gCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATgc                                                                                            >  3:192488/1‑149 (MQ=255)
                                                        cAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATg                                                                                             <  8:94245/147‑1 (MQ=255)
                                                                 tACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGAt                                                                                   >  4:15875/1‑148 (MQ=255)
                                                                   ccTGGCCGCGCTCGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCTTGCCGGCCGCCCAGTTCCTCGAAGACCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCTCAAGATGCGCCCCCTGAGATCTTGCGCGGCTgtt                                                                                   <  4:260672/146‑3 (MQ=255)
                                                                         cgcgCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGc                                                                           >  3:15303/1‑148 (MQ=255)
                                                                                 gggTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTccc                                                                  <  7:288178/149‑1 (MQ=255)
                                                                                   gTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTACCg                                                                 >  7:228543/1‑148 (MQ=255)
                                                                                             agcagGGTTACCTGTTCGCCAGACCCATGCCCGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGACCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCa                                                      >  6:58940/1‑149 (MQ=255)
                                                                                                         tGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCaa                                            >  3:282475/1‑147 (MQ=255)
                                                                                                           ttCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAgg                                        >  7:294476/1‑149 (MQ=255)
                                                                                                               ccAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCAtt                                    <  1:239259/149‑1 (MQ=255)
                                                                                                                cAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTg                                   >  2:315539/1‑149 (MQ=255)
                                                                                                                   aCCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGtt                                 >  8:78252/1‑148 (MQ=255)
                                                                                                                     ccATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCa                              >  4:338916/1‑149 (MQ=255)
                                                                                                                      cATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCAt                             >  1:69244/1‑149 (MQ=255)
                                                                                                                       aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGt                                  >  1:342669/1‑143 (MQ=255)
                                                                                                                       aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGt                                  <  2:342669/143‑1 (MQ=255)
                                                                                                                       aTGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATc                            >  5:275785/1‑149 (MQ=255)
                                                                                                                          ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                 <  5:54936/125‑1 (MQ=255)
                                                                                                                          ccggccgCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCgt                                                 >  6:54936/1‑125 (MQ=255)
                                                                                                                              ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                     <  5:105966/149‑1 (MQ=255)
                                                                                                                              ccgccCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATg                     <  3:15875/149‑1 (MQ=255)
                                                                                                                                   cAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTc                <  1:1648/149‑1 (MQ=255)
                                                                                                                                     gTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcct              <  6:177930/149‑1 (MQ=255)
                                                                                                                                      ttCCTCGCAGCCCTGGCCACTCGCCCGGGCGCAACGCAGGTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAACGGCATTGTTTCATCGGGCATGGCGTcctg             <  2:305247/149‑1 (MQ=255)
                                                                                                                                      ttCCTCGAAGCCCTGGCCACGCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctg             >  3:36383/1‑149 (MQ=255)
                                                                                                                                        gcgttgcgCCCGGGCGAGTGGCCAGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCACGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcc           >  6:203564/8‑149 (MQ=255)
                                                                                                                                         cgttgcgCCCGGGCGAGTGGCCAGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcct          >  4:350476/7‑149 (MQ=255)
                                                                                                                                         cTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTcctgcct          >  2:36762/1‑149 (MQ=255)
                                                                                                                                                cccTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCt   <  4:158279/149‑1 (MQ=255)
                                                                                                                                                 ccTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCtg  <  4:282420/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGTCGGGCAAGCTGGACCTGGGCATCGTTGCCGAGGGTGTGGAAACCGACCTGCAGCGCGACTACCTGGCCGCGCACGGGGTGGACTTCCAGCAGGGTTACCTGTTCGCCAGACCCATGCCGGCCGCCCAGTTCCTCGAAGCCCTGGCCACTCGCCCGGGCGCAACGCAGTTGCCGCAAGATGCGCCCCCTGAGATCATGCGCGGCTGATCAACCCGCTCAATTCCCGTATTCGAGGCATTCGTGTCAAAAGGCATTGTTTCATCGGTCATGGCGTCCTGCCTGTTCGCTG  >  NC_002947/3892392‑3892685

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: ‑