Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,346 C→G 11.2% Q261H (CAG→CAC)  PP_2941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,344,3460CG11.2% 117.4 / 5.6 44Q261H (CAG→CAC) PP_2941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base G (3/2);  total (23/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 = 9.99e-01

CCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATC  >  NC_002947/3344198‑3344490
                                                                                                                                                    |                                                                                                                                                  
ccATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGTAATGTCAGAAAAGCGGAGCTTGTGCTGCGATGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGGTTCAc                                                                                                                                                    <  1:172703/149‑1 (MQ=255)
  aTTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTCTGCAGCCGCTTCACt                                                                                                                                                   >  5:104796/1‑148 (MQ=255)
     ctcGAAGGTAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAATCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGaa                                                                                                                                                <  4:242725/148‑1 (MQ=255)
     ctcGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAg                                                                                                                                               >  2:33798/1‑149 (MQ=255)
     ctcGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAg                                                                                                                                               >  2:75973/1‑149 (MQ=255)
        gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGc                                                                                                                                              <  6:354211/147‑1 (MQ=255)
        gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGc                                                                                                                                              >  5:354211/1‑147 (MQ=255)
           gggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctg                                                                                                                                         <  1:33798/149‑2 (MQ=255)
                 gTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGATGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtca                                                                                                                                    <  6:104796/148‑1 (MQ=255)
                  ttATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcagc                                                                                                                                  <  5:203851/149‑1 (MQ=255)
                     tGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTTCAGCCGCTTCACTGAAGCGGCTCTAtcagctc                                                                                                                                <  4:250341/148‑1 (MQ=255)
                        cTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                        <  7:105777/137‑3 (MQ=255)
                        cTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                        >  8:105777/1‑135 (MQ=255)
                         tGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTg                                                                                                                                                  <  8:85467/126‑1 (MQ=255)
                         tGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTg                                                                                                                                                  >  7:85467/1‑126 (MQ=255)
                           aaCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCt                                                                                                                         <  5:213162/149‑1 (MQ=255)
                            aCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCt                                                                                                                         <  1:333965/148‑1 (MQ=255)
                                               aaaaGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                         >  1:134510/1‑145 (MQ=255)
                                                     ggAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCATGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtgt                                                                                               >  1:30704/1‑149 (MQ=255)
                                                            gtgCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCa                                                                                          <  8:20671/147‑1 (MQ=255)
                                                             tgctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTc                                                                                       <  4:123210/149‑1 (MQ=255)
                                                               ctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCATGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                     <  2:30704/149‑1 (MQ=255)
                                                                           ttCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTCCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         >  6:10023/1‑149 (MQ=255)
                                                                                         aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTc                                                         >  5:138669/1‑149 (MQ=255)
                                                                                          tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCg                                                        <  1:239139/149‑1 (MQ=255)
                                                                                           cGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGc                                                       >  3:314053/1‑149 (MQ=255)
                                                                                             gACCCACCGCAACCACGGGCAGGTCGACCACATGCGAATCTGTGCAGCCGCTTCACTGAAGAGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGAAATgg                                                                                                   >  7:201474/1‑105 (MQ=255)
                                                                                                   ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAAt                                                                                                     <  2:304168/97‑1 (MQ=255)
                                                                                                   ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAAt                                                                                                     >  1:304168/1‑97 (MQ=255)
                                                                                                       aaCCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAgcg                                           >  1:273231/1‑149 (MQ=255)
                                                                                                         ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         >  3:190040/1‑119 (MQ=255)
                                                                                                         ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         <  4:190040/119‑1 (MQ=255)
                                                                                                         ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTATATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTAGCACAGGTTGAgcgc                                          >  3:19453/1‑148 (MQ=255)
                                                                                                                   gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTAGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCgg                               >  1:284762/1‑149 (MQ=255)
                                                                                                                           ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCg                       >  2:154610/1‑149 (MQ=255)
                                                                                                                            cTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCTCTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTACTCGGCTTGTTCGc                      >  7:294233/1‑149 (MQ=255)
                                                                                                                               ctaaTCTGTGCAGCCGCTTCTCTGAAGCGGGTCTATCAGCTCAGGGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                    <  4:19453/146‑1 (MQ=255)
                                                                                                                               cGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                    <  2:273231/148‑1 (MQ=255)
                                                                                                                                   tCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGt               >  7:301722/1‑149 (MQ=255)
                                                                                                                                       tagaGCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgcg           >  8:210617/4‑149 (MQ=255)
                                                                                                                                        gCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                          <  7:112475/133‑1 (MQ=255)
                                                                                                                                        gCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                          >  8:112475/1‑133 (MQ=255)
                                                                                                                                           gCCGCTTCACTGAAGCGGCTCTATCAGATCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGtt       >  7:36784/1‑149 (MQ=255)
                                                                                                                                                  cACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATc  <  6:138669/147‑1 (MQ=255)
                                                                                                                                                   aCTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCCCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATc  <  1:94781/148‑1 (MQ=255)
                                                                                                                                                    cTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg          <  6:297807/137‑1 (MQ=255)
                                                                                                                                                    cTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg          >  5:297807/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGG‑‑CTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATC  >  NC_002947/3344198‑3344490

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