Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,835,842 Δ1 bp 8.9% coding (147/693 nt) PP_4249 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,835,8420C. 8.9% 129.8 / 8.8 45coding (147/693 nt)PP_4249hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base . (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

TTCAAGCAACACCTTCTGCGGGTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCATTAATCATGCT  >  NC_002947/4835698‑4835991
                                                                                                                                                |                                                                                                                                                       
ttCAAGCAACACCTTCTGCGGGTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGgcgc                                                                                                                                                    <  2:524413/150‑1 (MQ=255)
     gCAACACCTTCTGCGGGTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGcgcgcccg                                                                                                                                               >  1:323429/1‑150 (MQ=255)
                tGCGGGTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGcc                                                                                                                                    <  1:175663/150‑1 (MQ=255)
                   gggTCCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGa                                                                                                                                 >  2:933240/1‑150 (MQ=255)
                     gTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTCCGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATGTCGTCCATAAACCCGCTGTCCAGCTGCGCGCGCGCGCGATTGTCCAGaaa                                                                                                                               >  2:1164306/1‑149 (MQ=255)
                       tCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCa                                                                                                                             <  2:323429/150‑1 (MQ=255)
                        cggcggCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAg                                                                                                                            <  1:933240/150‑1 (MQ=255)
                         ggcggcCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCaag                                                                                                                           <  2:1027353/150‑3 (MQ=255)
                               cGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCCCGCCATGGGCCAGAATCAGGCAcccc                                                                                                                     >  2:880979/1‑150 (MQ=255)
                                       aCCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTc                                                                                                             <  2:86558/150‑1 (MQ=255)
                                               cccTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTACGTCCATGAACCCGCTGTCCATGGcgcgggcgcgccgatg                                                                                                                                       >  2:362839/1‑111 (MQ=255)
                                               cccTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTACGTCCATGAACCCGCTGTCCATGGcgcgggcgcgccgatg                                                                                                                                       <  1:362839/115‑5 (MQ=255)
                                                 cTGCGCTCGGCAATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCt                                                                                                   <  1:849601/150‑1 (MQ=255)
                                                         ggCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCcg                                                                                           >  1:333287/1‑150 (MQ=255)
                                                           ccATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCcgcg                                                                                         <  2:333287/150‑1 (MQ=255)
                                                             aTGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCgg‑‑a                                                                                     <  1:996016/150‑1 (MQ=255)
                                                                tACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCgg‑‑agga                                                                                  <  1:498163/150‑1 (MQ=255)
                                                                         tCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATGGcgcgggcgcgccgatg                                                                                                                                       <  1:52968/89‑5 (MQ=255)
                                                                         tCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATGGcgcgggcgcgccgatg                                                                                                                                       >  2:52968/1‑85 (MQ=255)
                                                                           gAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCaa                                                                       >  1:356724/1‑150 (MQ=255)
                                                                               cgcACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCt                                                                   <  1:108545/150‑1 (MQ=255)
                                                                                  aCTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTggg                                                                <  1:5871/150‑1 (MQ=255)
                                                                                    tACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGAc                                                              >  2:940333/1‑150 (MQ=255)
                                                                                     aCCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAgc                                                                           >  1:157084/1‑136 (MQ=255)
                                                                                     aCCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAgc                                                                           <  2:157084/136‑1 (MQ=255)
                                                                                              ccAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGAcc                                                    >  1:123334/1‑150 (MQ=255)
                                                                                               cAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTTCGACCt                                                   >  1:505657/1‑150 (MQ=255)
                                                                                                   agCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGc                                               <  2:505657/150‑1 (MQ=255)
                                                                                                    gCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGcc                                              <  1:632980/150‑1 (MQ=255)
                                                                                                     ccGCAAGCCTTTGCGCCATTTCGGCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGccc                                             <  1:1164306/150‑1 (MQ=255)
                                                                                                      cGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCa                                            >  1:903519/1‑150 (MQ=255)
                                                                                                       gCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATAGGCGCGCCCGCACGAATGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCAt                                           >  1:751434/1‑150 (MQ=255)
                                                                                                         aaGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCAt                                                                                                           >  1:652829/1‑86 (MQ=255)
                                                                                                         aaGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCAt                                                                                                           <  2:652829/86‑1 (MQ=255)
                                                                                                            ccTTTGCCCCATTTCGTCCATGAGCCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCCACCCGCGG‑‑AGGATCGCAGCGCCAGCGTGCGACCTTTGCCACCTCCGCCCagtttcc                                      <  2:620702/150‑7 (MQ=255)
                                                                                                              tttGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATccc                                    >  1:518497/1‑150 (MQ=255)
                                                                                                                            tCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCa                      >  2:120981/1‑150 (MQ=255)
                                                                                                                             ccATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCAc                     >  1:773941/1‑150 (MQ=255)
                                                                                                                             ccATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCAc                     <  2:751434/150‑1 (MQ=255)
                                                                                                                             ccATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCAc                     <  1:120981/150‑1 (MQ=255)
                                                                                                                                     ccGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGGTCTGCATGTGCCCAGGCCGCGGACATTTTCGACATTGGTCGATTGAGCCCCGGCAATACCGGCACTTTGCCCa             >  1:928817/1‑150 (MQ=255)
                                                                                                                                      cGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCAt            <  1:940333/150‑1 (MQ=255)
                                                                                                                                      cGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCAt            >  1:138471/1‑150 (MQ=255)
                                                                                                                                             cATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCATTAATCAt     >  2:1057446/1‑150 (MQ=255)
                                                                                                                                              aTCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCATTAATCATg    <  1:619135/150‑1 (MQ=255)
                                                                                                                                                cGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACAGGCACTTTGCCCATTAATCATGCt  >  1:733259/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
TTCAAGCAACACCTTCTGCGGGTTCGGCGGCCGCTTGCCACCGGTAACCCTGCGCTCGGCCATGTACGGGAACTCGAAGCGCACTACCGCTACCCCAAGAGCCGCAAGCCTTTGCGCCATTTCGTCCATGAACCCGCTGTCCATCGGCGCGCCCGCGCCATGGGCCAGAATCAGGCACCCCTTGTAGTCATCCTGACCTGCAACCCGCGG‑‑AGGATCGCAGCGCAAGCCTGGGACATTTGCGACCTTCGCCCATTGATCCCCGTCAATACCGGCACTTTGCCCATTAATCATGCT  >  NC_002947/4835698‑4835991

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