Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,305 C→T R395C (CGT→TGT)  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_0009131,196,3050CT96.1% 125.5 / NA 53R395C (CGT→TGT) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (31/18);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTT  >  NC_000913/1196167‑1196444
                                                                                                                                          |                                                                                                                                           
ggcaagtgaaaTCCGGGCTCTGTTATGCTCGGCGCTGGGATGATGCTGCGCCACATGGGTGGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCt                                                                                                                                         <  1:581544/135‑1 (MQ=21)
      gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGTCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGCGTGTCTgatgga                                                                                                                                   >  2:812776/1‑143 (MQ=255)
           tCCTGGCTCTATTTTTCTCTCCGCTGTGATGATGCTGCGACACATGGGTTGGACTGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                                                                                                              <  2:697142/143‑1 (MQ=21)
           tCCTGGCTCTATTATTCTCTCCGCTGAGATTATGCTGCGATATATTTGTTCTACCAAAGCGGCTGACTAAAAGGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                                                                                                              <  2:372092/143‑1 (MQ=21)
              tGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                                                                                                           >  2:730585/1‑143 (MQ=255)
              tGGCTCTATTATTCTCTCCACTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                                                                                                           >  2:136532/1‑143 (MQ=255)
                  tctATTATTCTATCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGTAGCGGCTGACTTAATTGTTAAAGGTATGGCAGGCGCAATCTAGGGGGAAACGGTAACCTACGAGTTCGAGCGCACGATGGATAgcgcggacct                                                                                                                       >  1:857460/1‑137 (MQ=255)
                    tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTCATTGTTAAAGGTATGGCAGGCGCAATCAACGCGAAAACCGTAACCTCTGACTTCGAGCGACTGATGGATGGCCCTaaccgcc                                                                                                                     >  1:89410/1‑138 (MQ=255)
                    tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAGTCAACGCGAAAACCGTAACCTATGACTTCGAGTGTGTGATGGATGGCGCTGAActgc                                                                                                                     >  2:433726/1‑143 (MQ=255)
                            tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAGACTGCTGAAATGt                                                                                                             >  1:637037/1‑143 (MQ=255)
                            tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGt                                                                                                             >  1:461512/1‑143 (MQ=255)
                            tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGt                                                                                                             <  2:834317/143‑1 (MQ=255)
                             ctccgcgcggatgctGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                                                                                                            <  1:36890/134‑1 (MQ=18)
                              tcCGCTGAGATGATGCTGCGCCACATGGGCTGGACCGAAGCGGCTGACTTAATTGTTCATGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTc                                                                                                           >  1:851593/1‑143 (MQ=255)
                               ccGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGACATGTTct                                                                                                          >  2:544986/1‑142 (MQ=255)
                                cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTCAACTGCTAAAACGTTCag                                                                                                         >  1:2775/1‑143 (MQ=255)
                                       atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAATGTATGGAAGGCGATATAAACGCGACAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAATGGTTCTGAGgttcg                                                                                                  >  1:498584/1‑141 (MQ=255)
                                       atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATCCTCAGAGCTTgg                                                                                                  >  1:76179/1‑143 (MQ=255)
                                              tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGATAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg                                                                                           >  1:527747/1‑143 (MQ=255)
                                                cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                                         >  1:72020/1‑143 (MQ=255)
                                                  ccACATTGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTACCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGACATGTTCAGAGTTTGGTGACGCGatca                                                                                       >  1:829327/1‑143 (MQ=255)
                                                   cacaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGGAACCTATGACTGCGAGTGGCTGATGTATGCCGCTAAACTGCGGAAATGTTCAGAGTTTGGTGACGCGatcat                                                                                      <  2:588491/143‑1 (MQ=18)
                                                   cacaTGGGGTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcat                                                                                      <  2:31955/143‑1 (MQ=255)
                                                     caTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCg                                                                                    <  2:291912/143‑1 (MQ=255)
                                                        gggCTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAGGGCGCAATCAGCGCGAGAACCTTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                 <  1:670864/143‑1 (MQ=18)
                                                         ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaaa                                                                                >  2:628060/1‑143 (MQ=255)
                                                          gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGGGTCTGATGGATGGCGCTAAACTGCTGCAATGTTAAGAGGTTGGTGgcgccgtgagtgaaagc                                                                               >  1:176864/1‑130 (MQ=18)
                                                                 cGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGACCGCGCTCCCGCTACGGCTGAAATGTTCAGCGTTTGGTGACGCGATCATCGAAAACACGTAAt                                                                        >  2:646804/1‑143 (MQ=18)
                                                                            aCTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTg                                                             >  2:385127/1‑143 (MQ=255)
                                                                              gtAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACGTGTAATGCCGGAGTTTGtt                                                           >  2:701822/2‑143 (MQ=255)
                                                                               tAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTa                                                          >  2:754739/1‑143 (MQ=255)
                                                                                    gTTACAGGTATTGAAGGCGCAATCAACGCGAAAGCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttt                                                     <  1:639821/143‑1 (MQ=255)
                                                                                     ttAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATtta                                                    >  2:781151/1‑143 (MQ=255)
                                                                                     ttAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAAACGGAATGCCGTAGTTTGTTAAATtta                                                    >  1:140618/1‑143 (MQ=255)
                                                                                      taactttggggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCGGTAGTTTGTTAAAtttct                                                   >  1:858338/10‑141 (MQ=21)
                                                                                           gTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACg                                              >  2:280623/1‑143 (MQ=255)
                                                                                             aTGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACggg                                            >  1:502484/1‑143 (MQ=255)
                                                                                               ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAg                                          >  1:586707/1‑143 (MQ=255)
                                                                                                gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGAGTCTGATGGATGACGCTAAACTGCTGAAATGTTCCGAGTCTGGTGACGCGATCATAGCAAACATGTAATCACGTAGGTTGTTAAATCTATTAACGGgcgc                                         >  1:584518/1‑143 (MQ=255)
                                                                                                    gcgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                     <  2:179532/143‑1 (MQ=255)
                                                                                                      gcAATCAACGTCAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                   <  1:616688/143‑1 (MQ=255)
                                                                                                      gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTGTGATTGATCGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                   >  1:189901/1‑143 (MQ=255)
                                                                                                               gcgAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgggg                          >  2:109997/1‑140 (MQ=255)
                                                                                                                 gaaattcGTAGCCTATGGCTGCGAGTGCCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgtt                        <  2:684425/137‑1 (MQ=25)
                                                                                                                   aacccGTAGCCTATGATCTCGAGTGTCTGATGGATGGCGCTAACCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaataa                      <  1:433726/140‑6 (MQ=21)
                                                                                                                    aaCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGCGCGTAACGCTCCCGTTGttttt                     >  2:55451/1‑143 (MQ=255)
                                                                                                                     aCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttttt                    <  2:72020/143‑1 (MQ=255)
                                                                                                                         tAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTa                >  2:147663/1‑143 (MQ=255)
                                                                                                                         tAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTa                <  1:730585/143‑1 (MQ=255)
                                                                                                                           aCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg              <  2:851593/143‑1 (MQ=255)
                                                                                                                                tGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgct         >  1:11616/1‑143 (MQ=255)
                                                                                                                                       gAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGtt  >  2:34727/1‑143 (MQ=255)
                                                                                                                                       gAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGtt  <  2:139375/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTT  >  NC_000913/1196167‑1196444

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