Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,472,775 C→T 7.8% N231N (AAC→AAT)  PP_3962 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,472,7750CT 7.8% 146.6 / 2.3 51N231N (AAC→AAT) PP_3962hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (2/2);  total (28/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 = 4.72e-01

CCAAAGCTCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGCTCTGGC  >  NC_002947/4472633‑4472914
                                                                                                                                              |                                                                                                                                           
ccAAAGCTCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGt                                                                                                                                       >  2:236478/1‑149 (MQ=255)
  aaaGCTCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTc                                                                                                                                     <  3:127734/149‑1 (MQ=255)
      cTCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGc                                                                                                                                 >  5:198678/1‑149 (MQ=255)
       tCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGc                                                                                                                                 >  6:258616/1‑148 (MQ=255)
           cGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATa                                                                                                                            <  1:236478/149‑1 (MQ=255)
                    aTTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTgg                                                                                                                    >  1:46696/1‑148 (MQ=255)
                    aTTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTgg                                                                                                                    >  1:9684/1‑148 (MQ=255)
                           tCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTg                                                                                                            >  7:339533/1‑149 (MQ=255)
                             ccTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGTTGTTCATCTCCATCGGCCCGCCACGTCGATATCAGACAAACCGAACAGCCGTTCACGCCAATACCAACTGGGGCAGTTGCa                                                                                                          >  8:60751/1‑149 (MQ=255)
                             ccTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGc                                                                                                           <  3:117407/148‑1 (MQ=255)
                              cTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCAt                                                                                                         >  4:146304/1‑149 (MQ=255)
                                  ggCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAg                                                                                                     <  5:258616/149‑1 (MQ=255)
                                   gCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCAt                                                                                                         >  6:374373/1‑144 (MQ=255)
                                   gCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCAt                                                                                                         <  5:374373/144‑1 (MQ=255)
                                          ccAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAATAGCCGTTCCAGCCAATAGCACCTGGGa                                                                                                                  >  8:185973/1‑128 (MQ=255)
                                          ccAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAATAGCCGTTCCAGCCAATAGCACCTGGGa                                                                                                                  <  7:191354/128‑1 (MQ=255)
                                          ccAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAATAGCCGTTCCAGCCAATAGCACCTGGGa                                                                                                                  <  7:185973/128‑1 (MQ=255)
                                          ccAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAATAGCCGTTCCAGCCAATAGCACCTGGGa                                                                                                                  >  8:191354/1‑128 (MQ=255)
                                                              cacGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTTAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGACCCTGGCGt                                                                                      >  3:354049/1‑136 (MQ=255)
                                                                 ggCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACgg                                                                      <  3:146304/149‑1 (MQ=255)
                                                                           ttCAGATCACCGATCCAACCGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTaaa                                                            <  2:362597/149‑1 (MQ=255)
                                                                              aGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCa                                                                                                                                   <  7:86519/75‑1 (MQ=255)
                                                                              aGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCa                                                                                                                                   >  8:86519/1‑75 (MQ=255)
                                                                              aGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGt                                                                  <  1:115076/140‑1 (MQ=255)
                                                                              aGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGt                                                                  >  2:115076/1‑140 (MQ=255)
                                                                              aGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCg                                                         >  7:188116/1‑149 (MQ=255)
                                                                                aTCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTgg                                                                                                                    >  1:152770/1‑88 (MQ=255)
                                                                                aTCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTgg                                                                                                                    <  2:152770/88‑1 (MQ=255)
                                                                                      gATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGTCAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTgg                                                 <  7:60751/149‑1 (MQ=255)
                                                                                          cAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCaaa                                             <  6:198678/149‑1 (MQ=255)
                                                                                           aaCAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCTTGGctttg                                            >  3:126225/1‑145 (MQ=255)
                                                                                               ggTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCaa                                        <  8:37399/149‑1 (MQ=255)
                                                                                                       tCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGc                                >  3:258433/1‑149 (MQ=255)
                                                                                                       tCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGACTGGCAAAGCCAAGCTACGGc                                >  8:138731/1‑149 (MQ=255)
                                                                                                        cTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGcc                               >  6:94791/1‑149 (MQ=255)
                                                                                                          gCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCt                             >  5:322212/1‑149 (MQ=255)
                                                                                                          gCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACATGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCt                             >  6:276930/1‑149 (MQ=255)
                                                                                                               ggCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTaa                         <  2:9684/148‑1 (MQ=255)
                                                                                                                gCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCa                                                       <  5:101985/117‑1 (MQ=255)
                                                                                                                gCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCa                                                       >  6:101985/1‑117 (MQ=255)
                                                                                                                gCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAgg                       >  4:314430/1‑149 (MQ=255)
                                                                                                                gCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAgg                       >  6:179878/1‑149 (MQ=255)
                                                                                                                  ccGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCt                     <  5:179878/149‑1 (MQ=255)
                                                                                                                      cACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCa                 <  6:322212/149‑1 (MQ=255)
                                                                                                                         gTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATcc              <  8:188116/149‑1 (MQ=255)
                                                                                                                                 cTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCttt                           >  7:49882/1‑128 (MQ=255)
                                                                                                                                 cTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCttt                           <  8:49882/128‑1 (MQ=255)
                                                                                                                                  tGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGctct     >  6:148333/1‑149 (MQ=255)
                                                                                                                                   gagaAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGCTCTg    >  3:332553/1‑149 (MQ=255)
                                                                                                                                    agaAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGCTCTgg   <  5:148333/149‑1 (MQ=255)
                                                                                                                                     gaAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGCTCTGGc  <  8:339533/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CCAAAGCTCAGCGGGCTGACATTCAAATCCCTGCGGCTGATCCCAAGGCGGATGGTTGGCAACACGGCAACGTTATTCAGATCACCGATCCAACAGGTGTTCATCTGCATCGGCCCGCCACGTCGATATCTGAGAAAGCGAACAGCCGTTCCAGCCAATAGCACCTGGGACAGTTGCATGTAGCTATGCCCCTGGCGTGCCTCTCCTGGGACGGCTGTCGTAAACCGCAGCCTGGCAAAGCCAAGCTACGGCCCTTTAAGGCTTCCATCCTCACGCTCTGGC  >  NC_002947/4472633‑4472914

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