Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,111,188 G→C 7.8% P108A (CCA→GCA)  PP_5361 ← protein P47K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,111,1880GC 7.8% 151.4 / 3.6 51P108A (CCA→GCA) PP_5361protein P47K
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base C (2/2);  total (23/28)
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.08e-01

CCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCA  >  NC_002947/6111047‑6111330
                                                                                                                                               |                                                                                                                                              
cccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGc                                                                                                                                         >  2:274085/1‑149 (MQ=255)
   ggCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCa                                                                                                                                      <  5:364514/149‑1 (MQ=255)
         ccAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATa                                                                                                                                >  4:319469/1‑149 (MQ=255)
         ccAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCTgagacctt                                                                                                                                <  1:137416/149‑5 (MQ=255)
          cAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                               >  8:44309/1‑149 (MQ=255)
          cAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                               >  4:71815/1‑149 (MQ=255)
                 tCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGa                                                                                                                        <  3:319469/149‑1 (MQ=255)
                     ctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                     <  8:67740/148‑1 (MQ=255)
                        ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATc                                                                                                                 <  5:16408/149‑1 (MQ=255)
                        ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCAAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTAGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGCACTCGAt                                                                                                                  >  5:169680/1‑148 (MQ=255)
                          tGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                               <  8:329767/133‑1 (MQ=255)
                          tGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                               >  7:329767/1‑133 (MQ=255)
                              aGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCaaa                                                                                                           >  1:267062/1‑149 (MQ=255)
                              aGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCaaa                                                                                                           >  5:132396/1‑149 (MQ=255)
                                  gcgcAGGAAGTTCAAGCCATCGACCACGGTTACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaataa                                                                                                        <  6:180041/148‑1 (MQ=255)
                                  gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGAGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAACCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaat                                                                                                          >  1:289378/1‑148 (MQ=255)
                                  gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaataa                                                                                                        <  3:71815/148‑1 (MQ=255)
                                   cgcAGGACTGTCTAGCCATCGACCACGATGACCTTGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCTCTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                      <  5:325052/149‑1 (MQ=255)
                                                aGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCtg                                                                                                                                       >  7:198780/1‑102 (MQ=255)
                                                aGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCtg                                                                                                                                       <  8:198780/103‑2 (MQ=255)
                                                 gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCtt                                                                                        <  8:97057/149‑1 (MQ=255)
                                                   cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCg                                                                                      <  5:23880/149‑1 (MQ=255)
                                                   cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCg                                                                                      <  7:44309/149‑1 (MQ=255)
                                                     tCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGGCTCAGCCACTGGCAGTGGCTCAGAGATACCGGGGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGc                                                                                    <  5:85372/149‑1 (MQ=255)
                                                        aCCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGc                                                                                                              <  8:337769/120‑1 (MQ=255)
                                                            cGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                             <  6:132396/149‑1 (MQ=255)
                                                            cGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                             <  2:267062/149‑1 (MQ=255)
                                                                  cagtggtggCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcc                                                                       >  6:126185/4‑149 (MQ=255)
                                                                                 gcgcCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCATCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTc                                                         <  2:289378/148‑1 (MQ=255)
                                                                                  cgcCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACATCCTCGAGCAGGTCtt                                                       >  8:18209/1‑149 (MQ=255)
                                                                                   gcCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTc                                                      <  1:274085/149‑1 (MQ=255)
                                                                                           ggACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTc                                                         <  8:339212/138‑1 (MQ=255)
                                                                                           ggACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTc                                                         >  7:339212/1‑138 (MQ=255)
                                                                                           ggACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACg                                              >  6:325111/1‑149 (MQ=255)
                                                                                            gACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTc                                                                               >  1:168413/1‑115 (MQ=255)
                                                                                            gACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTc                                                                               <  2:168413/115‑1 (MQ=255)
                                                                                            gACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGt                                             >  2:95304/1‑149 (MQ=255)
                                                                                             aCAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACATCCTCGAGCAGGTCTTCGCGCAACGTg                                            <  7:18209/149‑1 (MQ=255)
                                                                                               aGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGAGCAGCAAATAATCGAACCGGCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTg                                            >  5:368742/1‑147 (MQ=255)
                                                                                                       cGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAAt                                    <  1:358159/149‑1 (MQ=255)
                                                                                                             tGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGcc                              >  6:150490/1‑149 (MQ=255)
                                                                                                                     gcgGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTc                       <  5:325111/148‑1 (MQ=255)
                                                                                                                       ggAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtc                    >  4:182616/1‑149 (MQ=255)
                                                                                                                         aaaGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGc                                                              >  2:215964/1‑105 (MQ=255)
                                                                                                                         aaaGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGc                                                              <  1:215964/105‑1 (MQ=255)
                                                                                                                                   tctgAGCCACTGCCAGTGGCTGAGAGATACCGGTGGACTCGATCAGCAAATAATCGAACCGCCACTCTTCCGCCAGTCGCGCAACCTCCTCGAGCAGCTCATCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGAGCAGTTtc        >  5:255602/1‑149 (MQ=255)
                                                                                                                                   tctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtc        >  2:176505/1‑149 (MQ=255)
                                                                                                                                       aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtct       <  1:95304/146‑1 (MQ=255)
                                                                                                                                         ccACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa  <  5:126185/149‑1 (MQ=255)
                                                                                                                                          cACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCgcgc                                                                           <  6:143288/75‑1 (MQ=255)
                                                                                                                                          cACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCgcgc                                                                           >  5:143288/1‑75 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACT‑‑GCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCA  >  NC_002947/6111047‑6111330

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