Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,454 A→T 8.3% intergenic (+45/+27) cspD → / ← icd DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,520,4540AT 8.3% 131.7 / 3.2 48intergenic (+45/+27)cspD/icdDNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (17/27);  new base T (2/2);  total (19/29)
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.85e-01

CCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGA  >  NC_002947/4520307‑4520602
                                                                                                                                                    |                                                                                                                                                    
ccccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAACAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTa                                                                                                                                                      <  1:94281/149‑1 (MQ=255)
      gggCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTAGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                               >  8:177532/1‑149 (MQ=255)
          ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                               <  5:231281/145‑1 (MQ=255)
          ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                               >  6:231281/1‑145 (MQ=255)
                gCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGttt                                                                                                                                     >  3:159897/1‑149 (MQ=255)
                   gTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACCGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAtttttt                                                                                                                                  <  8:294200/149‑3 (MQ=38)
                   gTTGAAATCAAAGATGCCCCCGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCt                                                                                                                                  <  8:284651/149‑1 (MQ=255)
                    ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTg                                                                                                                                 <  7:258713/149‑1 (MQ=255)
                    ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTg                                                                                                                                 <  4:254691/149‑1 (MQ=255)
                       aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                               <  6:294213/148‑1 (MQ=255)
                        aaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTg                                                                                                                              <  4:159897/148‑1 (MQ=255)
                          tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                           >  4:293192/1‑149 (MQ=255)
                            cacgATGCCGCGGCCGCTGCGGCATAGGCTGACAGAGCCTTGGGCGAGCCGCAAAGCGAAACTCCAA‑ACGCCTGACTGCCTGCGCTCTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGACGCGACCGGTTTTCTGTTGCCtt                                                                                                                          <  6:189901/145‑1 (MQ=255)
                            aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTa                                                                                                                         <  2:125256/149‑1 (MQ=255)
                             aaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                        <  1:339970/149‑1 (MQ=255)
                              aGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACa                                                                                                                       <  1:317841/149‑1 (MQ=255)
                                   ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt                                                                                                                  <  5:162860/149‑1 (MQ=255)
                                   ccGCGCCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt                                                                                                                  <  8:337367/149‑1 (MQ=255)
                                     gcgGCCGCTGCGGCATAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGttt                                                                                                                <  6:6123/149‑1 (MQ=255)
                                         ccGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGAt                                                                                                             >  1:213328/1‑148 (MQ=255)
                                         ccGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATc                                                                                                            <  3:293192/149‑1 (MQ=255)
                                            cTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAtt                                                                                                         <  1:85594/149‑1 (MQ=255)
                                                           ccAGAGCCTTTGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                                          <  7:177532/149‑1 (MQ=255)
                                                           ccAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                                          <  2:213328/149‑1 (MQ=255)
                                                              gagCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAg                                                                                       >  7:116128/1‑149 (MQ=255)
                                                                     gggCGAGCCGCAAAGCGAAAATGCAA‑ACGCCTGACTGCCTCCGCTGTGCCTATAGACAAAAAAAGCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAgctgct                                                                                  <  6:209614/147‑1 (MQ=255)
                                                                              gCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                 >  3:32063/1‑107 (MQ=255)
                                                                              gCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                 <  4:32063/107‑1 (MQ=255)
                                                                                          tGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                        >  4:34180/1‑88 (MQ=255)
                                                                                          tGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                        <  3:34180/88‑1 (MQ=255)
                                                                                           gCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGt                                                                                                                                <  4:54623/79‑1 (MQ=255)
                                                                                           gCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGt                                                                                                                                >  3:54623/1‑79 (MQ=255)
                                                                                           gCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                                               >  1:210905/1‑144 (MQ=255)
                                                                                           gCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                                               <  2:210905/144‑1 (MQ=255)
                                                                                                  gCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTc                                                    <  4:248618/149‑1 (MQ=255)
                                                                                                     tGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGa                                                  >  3:246756/1‑148 (MQ=255)
                                                                                                           cctgcGGTGTGCCGATAGAAAAAAAAATCGGCCGATTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTa                                           <  3:185910/149‑1 (MQ=255)
                                                                                                                  tgtgCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACg                                    >  8:149479/1‑149 (MQ=255)
                                                                                                                          tAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGc                            <  8:116128/149‑1 (MQ=255)
                                                                                                                               aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                       >  1:313810/1‑149 (MQ=255)
                                                                                                                                 aaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTCCCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATg                     >  5:199030/1‑149 (MQ=255)
                                                                                                                                        ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt              >  7:208663/4‑149 (MQ=255)
                                                                                                                                         gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAtt             >  1:193391/3‑149 (MQ=255)
                                                                                                                                          ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg            <  7:149479/149‑1 (MQ=255)
                                                                                                                                            gCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAgg                                         >  6:94269/1‑118 (MQ=255)
                                                                                                                                            gCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAgg                                         <  5:94269/118‑1 (MQ=255)
                                                                                                                                              ttcttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGc        >  8:2457/1‑149 (MQ=255)
                                                                                                                                                cttTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGcc       <  8:208663/148‑1 (MQ=255)
                                                                                                                                                 tttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCt     >  2:218708/1‑149 (MQ=255)
                                                                                                                                                    aaaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCTCGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga  >  2:183411/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAA‑AAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGA  >  NC_002947/4520307‑4520602

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