Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,111,244 G→C 10.2% 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,2440GC10.2% 146.4 / 8.1 49P108A (CCA→GCA) PP_5361protein P47K
Reads supporting (aligned to +/- strand):  ref base G (19/25);  new base C (3/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111097‑6111392
                                                                                                                                                   |                                                                                                                                                    
tttCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCa                                                                                                                                                    <  1:324731/150‑1 (MQ=255)
    gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                <  1:763949/150‑1 (MQ=255)
       ccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTc                                                                                                                                             >  2:538460/1‑150 (MQ=255)
             ggCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGAt                                                                                                                                       <  2:338780/150‑1 (MQ=255)
              gCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATa                                                                                                                                      >  2:57570/1‑150 (MQ=255)
                   gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                 <  1:983507/150‑1 (MQ=255)
                            tgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTgg                                                                                                                                                      >  1:738433/1‑120 (MQ=255)
                            tgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTgg                                                                                                                                                      <  2:738433/120‑1 (MQ=255)
                             gcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGAGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATc                                                                                                                       <  2:384853/150‑1 (MQ=255)
                              ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCgg                                                                                                                                  <  2:147760/138‑1 (MQ=255)
                              ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCgg                                                                                                                                  >  1:147760/1‑138 (MQ=255)
                                    aGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaat                                                                                                                <  2:890435/150‑1 (MQ=255)
                                        gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                            <  1:57570/150‑1 (MQ=255)
                                              gAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCgg                                                                                                                                  <  1:534960/122‑1 (MQ=255)
                                              gAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCgg                                                                                                                                  >  2:534960/1‑122 (MQ=255)
                                              gAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACg                                                                                                      >  2:556371/1‑150 (MQ=255)
                                                       gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTc                                                                                             <  2:219412/150‑1 (MQ=255)
                                                         cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCg                                                                                                                                   <  1:899263/110‑1 (MQ=255)
                                                         cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCg                                                                                                                                   >  2:899263/1‑110 (MQ=255)
                                                                 aCGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                                   <  1:840351/150‑1 (MQ=255)
                                                                  cGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgcg                                                                                  <  2:984354/150‑1 (MQ=255)
                                                                         cATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcctc                                                                           >  1:786363/1‑150 (MQ=255)
                                                                         cATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcctc                                                                           >  1:1127569/1‑150 (MQ=255)
                                                                         cATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCCGTGGCAGCCATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcctc                                                                           <  2:471453/150‑1 (MQ=255)
                                                                         cATGGTGTCCAGTCGCGCCATATCGGAAAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGATCCGATCAGCAAATAATCGAACCGCCCCTCTTCGGCCAGTCGCGCAAcctc                                                                           <  2:471479/150‑1 (MQ=255)
                                                                              tgtCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGa                                                                      >  1:521870/1‑150 (MQ=255)
                                                                                  cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAg                                                                  <  2:786363/150‑1 (MQ=255)
                                                                                  cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCCCTTCGGCCAGTCGCGCAACCTCCTCGAGCAg                                                                  >  2:613505/1‑150 (MQ=255)
                                                                                   aGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAGGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAgg                                                                 <  1:538460/150‑1 (MQ=255)
                                                                                           cATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTcgcg                                                         >  1:4017/1‑150 (MQ=255)
                                                                                                 ggACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGt                                                   >  1:1131678/1‑150 (MQ=255)
                                                                                                     aGACTGCGACCCTGGTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGccc                                                                                                   <  1:41868/98‑1 (MQ=255)
                                                                                                     aGACTGCGACCCTGGTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGccc                                                                                                   >  2:41868/1‑98 (MQ=255)
                                                                                                              cccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgca                                                <  1:109295/140‑1 (MQ=255)
                                                                                                              cccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgca                                                >  2:109295/1‑140 (MQ=255)
                                                                                                               ccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGc                                     >  2:229271/1‑150 (MQ=255)
                                                                                                                       tcgCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTc                             <  2:4017/150‑1 (MQ=255)
                                                                                                                       tcgCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTc                             <  1:229271/150‑1 (MQ=255)
                                                                                                                              aaGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAACCGTTGCTCGTTTCgg                      >  1:323649/1‑149 (MQ=255)
                                                                                                                                 tGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGa                      <  2:610224/147‑1 (MQ=255)
                                                                                                                                 tGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGa                      >  1:610224/1‑147 (MQ=255)
                                                                                                                                 ggAAGGGCTCAGGCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATGGAGCCGCGCCTCGTCGGGCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCa                   <  2:1131678/149‑1 (MQ=255)
                                                                                                                                      gTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtc              >  2:1196237/1‑150 (MQ=255)
                                                                                                                                        ctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                                   >  1:988203/1‑79 (MQ=255)
                                                                                                                                        ctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                                   <  2:988203/79‑1 (MQ=255)
                                                                                                                                        ctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTGCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtctc            <  1:633733/150‑1 (MQ=255)
                                                                                                                                         tgaGCCACTGCCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCATCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTctct           >  1:997670/3‑150 (MQ=255)
                                                                                                                                            gCCACTGCCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCgcgc                                                        <  2:482153/102‑1 (MQ=255)
                                                                                                                                            gCCACTGCCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCgcgc                                                        >  1:482153/1‑102 (MQ=255)
                                                                                                                                             ccACTGCCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtctc                        <  2:187133/133‑1 (MQ=255)
                                                                                                                                             ccACTGCCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtctc                        >  1:187133/1‑133 (MQ=255)
                                                                                                                                                  ggCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCgg  <  2:521870/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111097‑6111392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑