Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 286,445 A→G 7.7% S351G (AGC→GGC)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,4450AG 7.7% 143.1 / 3.7 52S351G (AGC→GGC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

CTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGG  >  NC_000913/286297‑286592
                                                                                                                                                    |                                                                                                                                                   
cTACTCGCTGATCAGCCTGCTGATTTT‑GTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGagag                                                                                                                                                    >  2:167300/1‑149 (MQ=255)
    tCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCg                                                                                                                                                 >  8:53686/1‑149 (MQ=255)
    tCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCCGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTAAAAACCACCACGCCGCAGCAGTGGCTCATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCg                                                                                                                                                 >  2:225775/1‑149 (MQ=255)
      gCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGc                                                                                                                                                <  2:75339/148‑1 (MQ=255)
      gCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGc                                                                                                                                                >  1:75339/1‑148 (MQ=255)
             aGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTc                                                                                                                                         <  8:86222/148‑1 (MQ=255)
               ccTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGTAGAGCCGCCGCGGTcgcc                                                                                                                                      >  5:182902/1‑149 (MQ=255)
                            gTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGc                                                                                                                          <  7:53686/148‑1 (MQ=255)
                             tCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTg                                                                                                                        <  1:228367/149‑1 (MQ=255)
                                ccccGGCGGAGCACAACGCGATCATTTTTGCCCTCAACATCCTGTTCCAGTTCGTCTATAATACCACAACGCCGCTGCAGTTGCGGATGTCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGACGCCGCCTCGACGGGCTgggt                                                                                                                     >  8:250893/1‑147 (MQ=255)
                                 cccGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCCCCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGgtgt                                                                                                                    >  5:84696/1‑149 (MQ=255)
                                 cccGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGCCGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGgtgt                                                                                                                    >  5:26014/1‑149 (MQ=255)
                                       ggAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGgtg                                                                                                                     <  8:42331/142‑1 (MQ=255)
                                       ggAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGgtg                                                                                                                     >  7:42331/1‑142 (MQ=255)
                                            acaTCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCac                                                                                                             <  2:186248/145‑1 (MQ=255)
                                            acaTCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCac                                                                                                             >  1:186248/1‑145 (MQ=255)
                                              aTCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCacctacct                                                                                                       <  1:167300/149‑1 (MQ=255)
                                               tCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTg                                                                                                      <  2:244588/149‑1 (MQ=255)
                                                 gcgcTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATCCCACCACGCCGCTGCAGTGGCTGTTGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGtt                                                                                                    <  1:225775/149‑1 (MQ=255)
                                                 gcgcTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGtt                                                                                                    >  1:69369/1‑149 (MQ=255)
                                                 gcgcTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGtt                                                                                                    <  7:44594/149‑1 (MQ=255)
                                                       aTTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCt                                                                                              <  4:189204/149‑1 (MQ=255)
                                                              cccTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCa                                                                                                              <  3:136413/126‑1 (MQ=255)
                                                              cccTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCa                                                                                                              >  4:136413/1‑126 (MQ=255)
                                                                cTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTgg                                                                                     >  2:182801/1‑149 (MQ=255)
                                                                     cATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTgg                                                                                <  1:182801/149‑1 (MQ=255)
                                                                          tgttcctgttcGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCCGCCTGAAGATTGGCCTGGCGAt                                                                            >  2:161654/1‑148 (MQ=255)
                                                                          tgttcctgttcGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGAtt                                                                           >  3:132341/1‑149 (MQ=255)
                                                                                  ttcGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTggcggggc                                                                   <  5:11711/149‑1 (MQ=255)
                                                                                      tCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCggtg                                                               >  6:45690/1‑149 (MQ=255)
                                                                                        tttAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTtctc                                                                                                                <  7:137172/98‑1 (MQ=255)
                                                                                        tttAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTtctc                                                                                                                >  8:137172/1‑98 (MQ=255)
                                                                                        tttAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCggtggt                                                             >  7:210308/1‑149 (MQ=255)
                                                                                               ccaccaCGCCGCTGCAGGGGCTGATGGCTGCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTg                                                       <  6:26014/148‑1 (MQ=255)
                                                                                                caccacGCCGCTGCAGTGGCTGATTGGTTCTGACGTGGTCGACTCCGAGGATAGCCGCTGCGGGCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGa                                                     <  2:117745/149‑1 (MQ=255)
                                                                                                     cgccgcTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTg                                                <  6:84696/149‑1 (MQ=255)
                                                                                                           tGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTAc                                          >  2:194177/1‑149 (MQ=255)
                                                                                                             ctgTGGCTGATGGCTTCTGACTTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCGTGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTAc                                          <  7:250893/145‑1 (MQ=255)
                                                                                                             cAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTAc                                          <  8:210308/147‑1 (MQ=255)
                                                                                                             cAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACg                                         <  3:32280/148‑1 (MQ=255)
                                                                                                                tGGCTGATGGCTTGTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGGTGGTGTTCTCCACCTACGTGGTCAGCCTGAAGATTGGCCTGGCGATGGGCGGGGCGGGGGTGGGCTGGATCCTGGCGTACGTCac                                     >  2:62989/1‑148 (MQ=255)
                                                                                                                tGGCTGATGGCTTCGGACGTGGGGGACTACGAGGAGAGCCGCCGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCGGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGGACGTCaa                                     <  6:182902/149‑1 (MQ=255)
                                                                                                                      aTGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGa                                                                                                                               <  1:71097/53‑1 (MQ=255)
                                                                                                                      aTGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGa                                                                                                                               >  2:71097/1‑53 (MQ=255)
                                                                                                                           ttCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGc                            >  2:91140/1‑147 (MQ=255)
                                                                                                                                    ggtggACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGCAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCg                 >  1:46134/1‑149 (MQ=255)
                                                                                                                                        gACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGTTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGc              <  2:69369/148‑1 (MQ=255)
                                                                                                                                           tACGAGGAGGGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGcc                                                                                   <  5:113890/76‑1 (MQ=255)
                                                                                                                                           tACGAGGAGGGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGcc                                                                                   >  6:86905/1‑76 (MQ=255)
                                                                                                                                           tACGAGGAGGGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGcc                                                                                   >  6:113890/1‑76 (MQ=255)
                                                                                                                                           tACGAGGAGGGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGcc                                                                                   <  5:86905/76‑1 (MQ=255)
                                                                                                                                              gaggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGt       <  6:18131/149‑1 (MQ=255)
                                                                                                                                               aggagAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGtt      <  4:132341/149‑1 (MQ=255)
                                                                                                                                                  agagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAg   <  8:12793/149‑1 (MQ=255)
                                                                                                                                                  agagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAg   >  6:129441/1‑149 (MQ=255)
                                                                                                                                                   gagCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAgg  >  3:114593/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTACTCGCTGATCAGCCTGCTGATTTTCGTCACCCCGGCGGAGCACATCGCGCTCATTTTTGCCCTCAACATCCTGTTCCTGTTCGTCTTTAATACCACCACGCCGCTGCAGTGGCTGATGGCTTCTGACGTGGTGGACTACGAGGAGAGCCGCAGCGGTCGCCGCCTCGACGGGCTGGTGTTCTCCACCTACCTGTTCAGCCTGAAGATTGGCCTGGCGATTGGCGGGGCGGTGGTGGGCTGGATCCTGGCGTACGTCAACTATTCCGCCAGCAGCAGCGTGCAGCCGGTTGAGG  >  NC_000913/286297‑286592

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: ‑