Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 294,916 T→A 6.8% intergenic (‑117/+223) yagM ← / ← yagN CP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913294,9160TA 6.8% 205.4 / 4.6 59intergenic (‑117/+223)yagM/yagNCP4‑6 prophage; uncharacterized protein/CP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base A (2/2);  total (33/26)
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.78e-01

AAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAATGGTTTGAGATTT  >  NC_000913/294776‑295064
                                                                                                                                            |                                                                                                                                                    
aaaaTTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGa                                                                                                                                              >  1:1367820/1‑149 (MQ=255)
      aGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTtct                                                                                                                                        <  2:1367820/149‑1 (MQ=255)
            acTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTTAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACa                                                                                                                                  >  2:606111/1‑149 (MQ=255)
                 tAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGaaaa                                                                                                                             <  4:115223/149‑1 (MQ=255)
                       tttttAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAAt                                                                                                                            <  3:1329726/144‑1 (MQ=255)
                       tttttAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAAt                                                                                                                            >  4:1329726/1‑144 (MQ=255)
                                 aaCCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTAt                                                                                                                        <  1:1298540/138‑1 (MQ=255)
                                 aaCCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTAt                                                                                                                        >  2:1298540/1‑138 (MQ=255)
                                    ctctTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTAt                                                                                                                        >  1:184820/1‑135 (MQ=255)
                                    ctctTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTAt                                                                                                                        <  2:184820/135‑1 (MQ=255)
                                     tctTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGa                                                                                                         >  1:269446/1‑149 (MQ=255)
                                            tATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCttt                                                                                                  <  2:269446/149‑1 (MQ=255)
                                                            tactaTCTTAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATgg                                                                                  <  1:606111/149‑1 (MQ=255)
                                                            tactaTCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGt                                                                                                                          <  3:1432506/109‑1 (MQ=255)
                                                            tactaTCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGt                                                                                                                          >  4:1432506/1‑109 (MQ=255)
                                                                              aGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCg                                                                                                                                      >  2:491423/1‑79 (MQ=255)
                                                                              aGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCg                                                                                                                                      <  1:491423/79‑1 (MQ=255)
                                                                                      cccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTgg                                                              >  1:381773/1‑143 (MQ=255)
                                                                                      cccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTgg                                                              <  2:381773/143‑1 (MQ=255)
                                                                                      cccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAgg                                                        >  3:619831/1‑149 (MQ=255)
                                                                                       ccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATttat                                                                                                                                                    <  3:384588/56‑1 (MQ=255)
                                                                                       ccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATttat                                                                                                                                                    >  4:384588/1‑56 (MQ=255)
                                                                                       ccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGtttt                                                                <  4:139571/140‑1 (MQ=255)
                                                                                       ccATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGtttt                                                                >  3:139571/1‑140 (MQ=255)
                                                                                        cATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTTTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTa                                                      >  2:599963/1‑149 (MQ=255)
                                                                                         aTTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACa                                                                                                                                  <  2:384043/72‑1 (MQ=255)
                                                                                         aTTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACa                                                                                                                                  >  1:384043/1‑72 (MQ=255)
                                                                                          tttttCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAg                                                     >  3:220527/1‑148 (MQ=255)
                                                                                          tttttCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAg                                                     <  4:220527/148‑1 (MQ=255)
                                                                                                   gggCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAAt                                                                                    <  2:151381/108‑1 (MQ=255)
                                                                                                   gggCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAAt                                                                                    >  1:151381/1‑108 (MQ=255)
                                                                                                          aTCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGt                                    >  4:1247097/1‑149 (MQ=255)
                                                                                                           tCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTa                                                                                      <  1:634905/98‑1 (MQ=255)
                                                                                                           tCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTa                                                                                      >  2:634905/1‑98 (MQ=255)
                                                                                                           tCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTc                                   <  3:1247097/149‑1 (MQ=255)
                                                                                                               tataAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGa                                                                        <  1:489356/108‑1 (MQ=255)
                                                                                                               tataAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGa                                                                        >  2:489356/1‑108 (MQ=255)
                                                                                                                     tCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGc                         >  2:517032/1‑149 (MQ=255)
                                                                                                                         tAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGc                                                    >  3:1092158/1‑118 (MQ=255)
                                                                                                                         tAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGc                                                    <  4:1092158/118‑1 (MQ=255)
                                                                                                                         tAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCAtt                     >  4:1118293/1‑149 (MQ=255)
                                                                                                                         tAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCAtt                     >  2:1222193/1‑149 (MQ=255)
                                                                                                                           aaaTTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAg                                                     <  4:1035997/115‑1 (MQ=255)
                                                                                                                           aaaTTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAg                                                     >  3:1035997/1‑115 (MQ=255)
                                                                                                                            aaTTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTAc                  <  1:517032/149‑1 (MQ=255)
                                                                                                                            aaTTAAATCTGTAAGTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTAc                  >  2:1018934/1‑149 (MQ=255)
                                                                                                                              ttAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGa                >  3:111792/1‑149 (MQ=255)
                                                                                                                                aaaTCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAAt              >  3:1236687/1‑149 (MQ=255)
                                                                                                                                  gtCTGTAATAAATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGaa               >  3:188685/2‑146 (MQ=255)
                                                                                                                                  gtCTGTAATAAATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGaa               <  4:188685/145‑1 (MQ=255)
                                                                                                                                   tCTGTAATAAATTACAGTCTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAgg                                            <  2:883627/116‑1 (MQ=255)
                                                                                                                                   tCTGTAATAAATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAgg                                            >  1:883627/1‑116 (MQ=255)
                                                                                                                                    cTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTg                                                               >  3:1350425/1‑96 (MQ=255)
                                                                                                                                    cTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTg                                                               <  4:1350425/96‑1 (MQ=255)
                                                                                                                                     tGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGAttt                                                                     <  1:508199/89‑1 (MQ=255)
                                                                                                                                     tGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGAttt                                                                     >  2:508199/1‑89 (MQ=255)
                                                                                                                                     tGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGa                               <  4:422284/127‑1 (MQ=255)
                                                                                                                                     tGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGa                               >  3:422284/1‑127 (MQ=255)
                                                                                                                                         aTTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGa                >  2:545271/1‑138 (MQ=255)
                                                                                                                                         aTTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGa                <  1:545271/138‑1 (MQ=255)
                                                                                                                                         aTTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAATGGTTTgaga     >  3:1395118/1‑149 (MQ=255)
                                                                                                                                          tttATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGtt                                                                  <  3:538312/87‑1 (MQ=255)
                                                                                                                                          tttATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGtt                                                                  >  4:538312/1‑87 (MQ=255)
                                                                                                                                          tttATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGc                                                    >  4:238903/1‑101 (MQ=255)
                                                                                                                                          tttATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGc                                                    <  3:238903/101‑1 (MQ=255)
                                                                                                                                          tttATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAATGGTTTGAGAt    >  3:1446441/1‑149 (MQ=255)
                                                                                                                                            tattaCAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAATGGTTTGAGAttt  <  3:1118293/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
AAAATTAGTAACACTATTAGACATTTTTAAGGAAACCTCTTTGTTATCTCAATATAAAACTACTATCTGAAATTTCTCAGTGTTGCCCCATTTTTCTCGGGGCAACATCGCTATAAATCAATAAAATTAAATCTGTAATTTATTACAGACTTTCTCGTACAGAAAATGTATCTATAAAATCGAGGACATCTTTACGTAACCAGTAATGGCGGCCATTGATTTGTTTTGGGCGAGGTAGCCCCTCAGGCAATTTGTCCTGACAAAGCCATTTACGAATGGTTTGAGATTT  >  NC_000913/294776‑295064

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