Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,516,644 C→T 13.7% intergenic (+641/‑494) 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,6440CT13.7% 100.2 / 7.9 45intergenic (+641/‑494)PP_5576/PP_5578hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (3/3);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GAAGAGACCGATCCGCTGGGCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGA  >  NC_002947/3516495‑3516793
                                                                                                                                                      |                                                                                                                                                     
gtagaaACCGAACTACTGGGCCGCAAGACCACCTACGAATATCATAACCTAACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTCGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGc                                                                                                                                                       >  2:470501/3‑150 (MQ=17)
          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                                                                                                            <  1:508015/123‑1 (MQ=40)
                                                                      tCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGcaca                                                                                                            >  2:508015/1‑123 (MQ=40)
                                                                      tCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCa                                                                                 >  1:1104380/1‑150 (MQ=39)
                                                                       cAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCaa                                                                                <  2:52718/150‑1 (MQ=38)
                                                                              ccGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCtg                                                                         >  1:63735/1‑150 (MQ=39)
                                                                                            cTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGc                                                                                               >  1:148521/1‑114 (MQ=39)
                                                                                            cTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGc                                                                                               <  2:148521/114‑1 (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                                           >  2:72336/1‑150 (MQ=32)
                                                                                                            gacgacAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGc                                           >  1:515150/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                                   <  1:720795/139‑1 (MQ=35)
                                                                                                                               tCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGa                                   >  2:720795/1‑139 (MQ=35)
                                                                                                                                 gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                      >  1:438646/1‑150 (MQ=35)
                                                                                                                                 gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                      >  1:807544/1‑150 (MQ=35)
                                                                                                                                 gCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAag                      <  1:470501/150‑1 (MQ=35)
                                                                                                                                  tcgctgttcaggtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTa                                                                                         <  1:902940‑M1/68‑1 (MQ=255)
                                                                                                                                  tcgctgttcaggtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTa                                                                                         >  2:902940‑M1/15‑82 (MQ=255)
                                                                                                                                         cGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGta              <  1:66891/150‑1 (MQ=39)
                                                                                                                                             gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGa                                                                                                  >  2:704525‑M1/4‑62 (MQ=255)
                                                                                                                                             gtaCTCGGTCATCTGGCCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGa                                                                                                  <  1:704525‑M1/59‑1 (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)
                                                                                                                                                      |                                                                                                                                                     
GAAGAGACCGATCCGCTGGGCCGCAAGACCACCTACGAATATCATCACCTGACCACGCTGGTGAACCAGGTCAGCTACCCGGATGGCAGCACCTGGCGAGCGCGCTACGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGCCTGCCGCACACCATCATCGATGCGACCTACAAGTCCAAGTACCTGTGGTGGAACACCCTGGCCCAGGTCGAACGCTATCAGGACTGCTCGGGCAAGAGCACGTACTACCGCTATGA  >  NC_002947/3516495‑3516793

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: