Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,769,648 (C)6→5 intergenic (+162/‑55) yfjQ → / → yfjR CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,769,6430C.97.1% 134.4 / ‑2.4 35intergenic (+157/‑60)yfjQ/yfjRCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (17/17);  minor base A (0/1);  total (17/18)
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.98e-01

GTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769507‑2769777
                                                                                                                                        |                                                                                                                                      
gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:987782/139‑1 (MQ=255)
gttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  1:558461/139‑1 (MQ=255)
 ttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACAACTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:202673/139‑1 (MQ=255)
         ggACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  2:365061/1‑130 (MQ=255)
         ggACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  1:365061/130‑1 (MQ=255)
          gACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTaaa                                                                                                                           >  2:1127637/1‑139 (MQ=255)
                 ccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATca                                                                                                                    >  1:819168/1‑139 (MQ=255)
                                     cTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                <  1:1127637/139‑1 (MQ=255)
                                     cTCATGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                <  1:332461/139‑1 (MQ=255)
                                      tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAAc                                                                                               >  1:126654/1‑139 (MQ=255)
                                        gtgtGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  2:1033408/100‑1 (MQ=255)
                                        gtgtGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  1:1033408/1‑100 (MQ=255)
                                         tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCAcat                                                                                                                 <  2:771703/118‑1 (MQ=255)
                                         tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCAcat                                                                                                                 >  1:771703/1‑118 (MQ=255)
                                          gtgCCTCAATCCCCCGGTTATAGCTTTTACCCCCCGCTACATCTGGCTTTTGCAGAGATAAAAAAAAGCTTCTGCGTTATCCTTACCCTGTACACCCCGTCTTTAGAGTAATCACATCATTTTCAATCAGTTACGTTTc                                                                                           <  2:844750/139‑1 (MQ=18)
                                                 aaTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGTGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCTGTCAACTTTCCTGgaga                                                                                    <  2:819168/139‑1 (MQ=18)
                                                  aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGaa                                                                                   >  2:35793/1‑139 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATc                                                                                                                     >  2:401173/1‑100 (MQ=255)
                                                       ccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATc                                                                                                                     <  1:401173/100‑1 (MQ=255)
                                                            tataGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGa                                                                         <  1:35793/139‑1 (MQ=255)
                                                             ataGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        <  2:126654/139‑1 (MQ=255)
                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtca                                                                                                         <  1:577668/99‑1 (MQ=255)
                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTtca                                                                                                         >  2:577668/1‑99 (MQ=255)
                                                                                       gCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTaaa                                                                                                                           <  1:1239393/62‑1 (MQ=255)
                                                                                       gCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTaaa                                                                                                                           >  2:1239393/1‑62 (MQ=255)
                                                                                            tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                <  1:531549/84‑1 (MQ=255)
                                                                                            tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                >  2:531549/1‑84 (MQ=255)
                                                                                            tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        >  1:628731/1‑108 (MQ=255)
                                                                                            tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGac                                                                        <  2:628731/108‑1 (MQ=255)
                                                                                                      aaaaaTAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctgt                               <  1:555105/139‑1 (MQ=255)
                                                                                                              tttCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGtt                                                                                                  <  2:220635/64‑1 (MQ=255)
                                                                                                              tttCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGtt                                                                                                  >  1:220635/1‑64 (MQ=255)
                                                                                                                              aCCCTGTCCGCCCCCTCTTTAAAGCAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggct                                             <  1:235274/101‑1 (MQ=255)
                                                                                                                              aCCCTGTCCGCCCCCTCTTTAAAGCAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggct                                             >  2:235274/1‑101 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAAcg                                                            >  2:145717/1‑81 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAAcg                                                            <  1:145717/81‑1 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGtt        <  2:1388408/133‑1 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGtt        >  1:1388408/1‑133 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGc      <  2:995567/135‑1 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGc      >  1:995567/1‑135 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCgagt  >  2:390240/1‑136 (MQ=255)
                                                                                                                                   gTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  2:103006/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769507‑2769777

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