Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 1,195,468 G→T V116F (GTT→TTT)  icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,195,4680GT96.3% 79.3 / NA 27V116F (GTT→TTT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (14/12);  minor base C (0/1);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGAAAAATCCACACAGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACGTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATC  >  NC_000913_3_bme_pgi/1195332‑1195611
                                                                                                                                        |                                                                                                                                               
cGGTGAAAAATCCACACAGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGcc                                                                                                                                     >  1:2109/1‑149 (MQ=255)
  gTGAAAAATCCACACAGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAg                                                                                                                                   <  2:104370/149‑1 (MQ=255)
         aTCCACACAGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTCGTGGCGTTCTCGCTCTGAACTTTGCCCTGCGCCAGGAACTgg                                                                                                                            <  1:570999/149‑1 (MQ=255)
                 aGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTAc                                                                                                                    <  1:713151/149‑1 (MQ=255)
                                aCGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCg                                                                                                     <  2:42339/149‑1 (MQ=255)
                                     tGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTGAAGTTCCGCTGATCACTCCGGTTGGTGGCGGTATTCGCTCCCTGAAGTTTGACCTGCGTCAGAAACTGGATCTCTACAGCTGCCTGCGTCGGATACg                                                                                                >  2:251770/1‑149 (MQ=255)
                                     tGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACAACTCCGGTTGGTGGCGGTATTCGCTCTCTGAAATTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCGTGCGTCGGGTGCg                                                                                                >  1:518793/1‑149 (MQ=255)
                                      ggaggcctgtTCAAACTCGTGATGTGATTGGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGTTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGt                                                                                               <  2:74776/145‑1 (MQ=255)
                                      ggCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCCCTGAACTTTGCCCTGCGCCAGGAACAGGATCTCTACACCCGCCTTCGACCGCTACGt                                                                                               >  2:726890/1‑149 (MQ=255)
                                              gctgAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTCCTCTCa                                                                                       >  1:524373/1‑149 (MQ=255)
                                                           aTCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGc                                                                          >  2:463071/1‑149 (MQ=255)
                                                                 ttCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGtt                                                                    >  1:519490/1‑149 (MQ=255)
                                                                  tCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTa                                                                   >  2:288099/1‑149 (MQ=255)
                                                                      gAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGGATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCGAAGCCCGGTTAAaca                                                               >  1:195223/1‑149 (MQ=255)
                                                                              cgTTGCGATTAAAGGTCCGCTGGCCACTCCGGTTGGTGGCGGTAGTCGCCCTCTGAACCTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAAc                                                       <  1:345492/149‑1 (MQ=255)
                                                                                 tGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGa                                                    <  2:298941/149‑1 (MQ=255)
                                                                                  gCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGAc                                                   >  1:419168/1‑149 (MQ=255)
                                                                                  gCCATTAAAGGACCGCTGACCACTCCGGTTGGTGGCGGGATTCGCCCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGTCTGCGTCCGGTACGTTACTTCGAGCGCGGACTTGGCCCGGTTAAACACCCTGAACTGAc                                                   <  1:378461/149‑1 (MQ=255)
                                                                                           ggTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGt                                          <  1:463071/149‑1 (MQ=255)
                                                                                                  tGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCTAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCAGTCCAAGCCCGGTTAAACACCCTGAACTGGCCGATATGGTTATCTTc                                   >  2:213755/1‑149 (MQ=255)
                                                                                                           cGGTTGGTGGCGGTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAAc                          >  2:337098/1‑149 (MQ=255)
                                                                                                                    gtggTATTCGATCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGAc                 <  1:464508/147‑1 (MQ=255)
                                                                                                                      ggTGTTCGCTCTCTGAACTTTGCCATGCGCCAGCAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACAt               <  2:320766/149‑1 (MQ=255)
                                                                                                                      ggTATTCGCTCTCTGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACAt               >  1:739481/1‑149 (MQ=255)
                                                                                                                              atctctGAACTTTGCCCTGCGCCAGAAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAAATGCCCGATATGGTTACCTTCCGTGAAAACTCGGAAGACATTTATGCgg       <  1:732343/148‑1 (MQ=255)
                                                                                                                                   tGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATc  <  1:337098/149‑1 (MQ=255)
                                                                                                                                   tGAACTTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATc  >  1:667344/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CGGTGAAAAATCCACACAGGTTTATGGTCAGGACGTCTGGCTGCCTGCTGAAACTCTTGATCTGATTCGTGAATATCGCGTTGCCATTAAAGGTCCGCTGACCACTCCGGTTGGTGGCGGTATTCGCTCTCTGAACGTTGCCCTGCGCCAGGAACTGGATCTCTACATCTGCCTGCGTCCGGTACGTTACTATCAGGGCACTCCAAGCCCGGTTAAACACCCTGAACTGACCGATATGGTTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATC  >  NC_000913_3_bme_pgi/1195332‑1195611

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