Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,516,648 A→C 12.2% intergenic (+645/‑490) PP_5576 → / → PP_5578 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,516,6480AC12.2% 112.6 / 11.0 50intergenic (+645/‑490)PP_5576/PP_5578hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/17);  new base C (3/3);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

ATCCGCTGGGCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGACGAG  >  NC_002947/3516505‑3516797
                                                                                                                                                |                                                                                                                                                     
aTCCGCTGGGCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCg                                                                                                                                                 >  2:1093624/1‑150 (MQ=255)
  gcggtagtacgtgctcttgcccgagcagtcctgatagcgttcgacctgggccagggtgttccaccacaggtacttggacttgtaggtcgcatcgatgatggtgtgcggcaggccgtcatcgctgttcaggtaCTCGGCCAGATGACCGAg                                                                                                                                               >  2:409339‑M1/133‑150 (MQ=255)
         gCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGa                                                                                                                                        >  2:377977/1‑150 (MQ=255)
                gACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGa                                                                                                                                 >  2:218399/1‑150 (MQ=255)
                           aaTATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACAGAGCACCTGAACAGCGATGACGACCTGc                                                                                                                      >  2:1106332/1‑150 (MQ=255)
                              atcatcACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCATCGGCCAGATGACCGAGTACCTGAACAGCGATGACGACCTGCCGc                                                                                                                   <  1:1106332/150‑1 (MQ=255)
                               tcatcaCCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGca                                                                                                                  <  1:377977/150‑1 (MQ=255)
                                  tcaCCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATCCATAAGTTGCGCGGAGGTGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCacac                                                                                                               >  2:870142/1‑150 (MQ=14)
                                         aCCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCatcatc                                                                                                        <  1:218399/150‑1 (MQ=255)
                                                    gAACCAGGTCAGCTACCCGGATGGTAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTa                                                                                             >  1:327287/1‑150 (MQ=255)
                                                          ggTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTc                                                                                       >  2:536310/1‑150 (MQ=38)
                                                            tCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGcaca                                                                                                                >  2:508015/1‑123 (MQ=40)
                                                            tCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGcaca                                                                                                                <  1:508015/123‑1 (MQ=40)
                                                            tCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCa                                                                                     >  1:1104380/1‑150 (MQ=39)
                                                             cAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCaa                                                                                    <  2:52718/150‑1 (MQ=38)
                                                                    ccGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCtg                                                                             >  1:63735/1‑150 (MQ=39)
                                                                                  cTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGc                                                                                                   <  2:148521/114‑1 (MQ=39)
                                                                                  cTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGc                                                                                                   >  1:148521/1‑114 (MQ=39)
                                                                                   tGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTg                                                              >  1:317667/1‑150 (MQ=35)
                                                                                    ggCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTgg                                                             <  1:1093624/150‑1 (MQ=34)
                                                                                        aGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCa                                                         <  2:63735/150‑1 (MQ=34)
                                                                                         gcgcgcTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAg                                                        <  1:536310/150‑1 (MQ=34)
                                                                                            cgcTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTc                                                     >  2:462728/1‑150 (MQ=33)
                                                                                                  gacgacAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGc                                               >  1:515150/1‑150 (MQ=32)
                                                                                                  gacgacAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGc                                               >  2:72336/1‑150 (MQ=32)
                                                                                                        aaGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGaa                                                  >  2:420222/1‑141 (MQ=11)
                                                                                                        aaGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGaa                                                  <  1:420222/141‑1 (MQ=11)
                                                                                                        aaGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAg                                         <  2:974416/150‑1 (MQ=35)
                                                                                                             cAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTg                                    <  2:1104380/150‑1 (MQ=35)
                                                                                                                ccTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTc                                 >  2:519524/1‑150 (MQ=35)
                                                                                                                ccTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTc                                 >  2:66891/1‑150 (MQ=35)
                                                                                                                 cTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCg                                >  1:721471/1‑150 (MQ=35)
                                                                                                                     tCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGa                                       >  2:720795/1‑139 (MQ=35)
                                                                                                                     tCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGa                                       <  1:720795/139‑1 (MQ=35)
                                                                                                                       gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                          >  1:438646/1‑150 (MQ=35)
                                                                                                                       gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                          <  1:470501/150‑1 (MQ=35)
                                                                                                                       gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                          >  1:807544/1‑150 (MQ=35)
                                                                                                                        tcgctgttcaggtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTa                                                                                             >  2:902940‑M1/15‑82 (MQ=255)
                                                                                                                        tcgctgttcaggtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTa                                                                                             <  1:902940‑M1/68‑1 (MQ=255)
                                                                                                                               cGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGta                  <  1:66891/150‑1 (MQ=39)
                                                                                                                                   gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGa                                                                                                      <  1:704525‑M1/59‑1 (MQ=255)
                                                                                                                                   gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGa                                                                                                      >  2:704525‑M1/4‑62 (MQ=255)
                                                                                                                                   gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGaa                                                  >  2:865137‑M1/4‑114 (MQ=255)
                                                                                                                                   gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGaa                                                  <  1:865137‑M1/111‑1 (MQ=255)
                                                                                                                                         CGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTa         <  1:661719‑M1/150‑1 (MQ=255)
                                                                                                                                            ccAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGa      >  1:553078/1‑150 (MQ=39)
                                                                                                                                             cAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACGAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGAc     >  1:27927/1‑150 (MQ=38)
                                                                                                                                             cAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGAc     <  2:327287/150‑1 (MQ=39)
                                                                                                                                             cAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTCGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGAc     >  2:152993/1‑150 (MQ=38)
                                                                                                                                               gATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGGACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGACGa   >  2:244646/1‑150 (MQ=38)
                                                                                                                                                aTGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGACGAg  <  2:661719/150‑1 (MQ=40)
                                                                                                                                                |                                                                                                                                                     
ATCCGCTGGGCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGACGAG  >  NC_002947/3516505‑3516797

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