Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,774,325 A→G A336A (GCA→GCG)  yfjW → CP4‑57 prophage; uncharacterized protein YfjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,774,3250AG96.7% 205.0 / ‑2.0 61A336A (GCA→GCG) yfjWCP4‑57 prophage; uncharacterized protein YfjW
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (23/36);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

CATACAAGGTGAATCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCACAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTGCACCTCTATCTAG  >  NC_000913/2774196‑2774450
                                                                                                                                 |                                                                                                                             
cATACAAGGTGAATCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCaa                                                                                                                            <  1:87846/133‑1 (MQ=255)
  tACAAGGTGAATCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAAtt                                                                                                                          >  1:345254/1‑133 (MQ=255)
       ggTGAATCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATgg                                                                                                                     >  2:529369/1‑133 (MQ=255)
            aTCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCt                                                                                                                >  2:579718/1‑133 (MQ=255)
              ccGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCttt                                                                                                              <  1:392781/133‑1 (MQ=255)
               cGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTc                                                                                                             <  1:96225/133‑1 (MQ=255)
                  tgatgtgatgCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAAAACGCAATTGATGGAAGCTTTCGCa                                                                                                          >  2:311741/1‑133 (MQ=255)
                   gatgtgatgCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCAt                                                                                                         >  1:50383/1‑133 (MQ=255)
                     tgtgatgCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCAttt                                                                                                       <  1:61321/133‑1 (MQ=255)
                       tgatgCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTAt                                                                                                     <  1:611376/133‑1 (MQ=255)
                          tgCCTAAATCAGCATCCACGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGcc                                                                                                  <  1:264092/133‑1 (MQ=255)
                           gCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCa                                                                                                 <  2:50383/133‑1 (MQ=255)
                                 agcagcATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATaaca                                                                                           >  2:632506/1‑133 (MQ=255)
                                 agcagcATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATaaca                                                                                           <  1:529369/133‑1 (MQ=255)
                                 agcagcATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATaaca                                                                                           <  1:632506/133‑1 (MQ=255)
                                 agcagcATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATaaca                                                                                           <  2:595700/133‑1 (MQ=255)
                                      caTCCAGGTTGAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGaa                                                                                      >  2:497322/1‑133 (MQ=255)
                                      caTCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGaa                                                                                      >  2:621627/1‑133 (MQ=255)
                                      caTCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGaa                                                                                      >  1:552879/1‑133 (MQ=255)
                                                gaGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTa                                                                            <  1:497322/132‑1 (MQ=255)
                                                aaGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTa                                                                            <  2:42588/133‑1 (MQ=255)
                                                  gTTTGAGCTAGCAAGTTGATAAATTTATGCAAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACa                                                                          <  1:311741/133‑1 (MQ=255)
                                                  gTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACa                                                                          <  2:552879/133‑1 (MQ=255)
                                                  gTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCACAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACa                                                                          <  1:364982/133‑1 (MQ=255)
                                                       agcaagcaAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTTCATGAtt                                                                     >  1:184219/1‑133 (MQ=255)
                                                       agcaagcaAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGAtt                                                                     >  1:354249/1‑133 (MQ=255)
                                                           agcaagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatac                                                                 >  1:482489/1‑133 (MQ=255)
                                                            gcaagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTataca                                                                >  1:37471/1‑133 (MQ=255)
                                                            gcaagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTataca                                                                <  1:579718/133‑1 (MQ=255)
                                                            gcaagATGATAAATTCATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTTCATGATTataca                                                                <  2:184219/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  1:621627/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  1:527395/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  2:482489/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  2:473905/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  2:37471/133‑1 (MQ=255)
                                                              aagATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTatacata                                                              <  2:561586/133‑1 (MQ=255)
                                                                  tgatAAATTTATGCTTGCCTCAAGTGTTACCGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGat                                                          <  1:223587/133‑1 (MQ=255)
                                                                          ttATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATataataat                                                  >  1:318889/1‑133 (MQ=255)
                                                                             tGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACAAATGGTATATTAGATGGTTATACATACGATATAATAATaaa                                               >  2:223587/1‑133 (MQ=255)
                                                                                   ccTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATaa                                         >  1:642797/1‑133 (MQ=255)
                                                                                       aaGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATa                                          <  2:588569/128‑1 (MQ=255)
                                                                                       aaGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATa                                          >  1:588569/1‑128 (MQ=255)
                                                                                        aGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAat                                    >  1:488932/1‑133 (MQ=255)
                                                                                            ttACTGACTCTCAAATAAAATTTAAATCAAACAATGCACAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATGCGATATATTAATAACTGATAAAAAATATaa                                >  1:205652/1‑133 (MQ=255)
                                                                                                    tctcAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGaaa                        >  1:405101/1‑133 (MQ=255)
                                                                                                     ctcAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGaaaa                       <  2:354249/133‑1 (MQ=255)
                                                                                                     ctcAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGaaaa                       <  1:38538/133‑1 (MQ=255)
                                                                                                      tcAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAAtata                                 <  1:535536/122‑1 (MQ=255)
                                                                                                      tcAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAAtata                                 >  2:535536/1‑122 (MQ=255)
                                                                                                       cAAATAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATc                     >  2:158788/1‑133 (MQ=255)
                                                                                                        aaaTAAAATTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCa                    <  1:575716/133‑1 (MQ=255)
                                                                                                            aaaaTTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGtt                <  2:405101/133‑1 (MQ=255)
                                                                                                             aaaTTTAAATTAAACTATGCGCAAATGAAGGTTGCTTTCGCTTTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTg               <  2:215294/133‑1 (MQ=255)
                                                                                                             aaaTTTAAATCAAACTATGCGCAAATGACGGTCGCTTTCGCGTCTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTg               <  2:256230/133‑1 (MQ=255)
                                                                                                             aaaTTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTg               <  2:33900/133‑1 (MQ=255)
                                                                                                             aaaTTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTg               <  1:512512/133‑1 (MQ=255)
                                                                                                             aaaTTTAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTg               <  2:642797/133‑1 (MQ=255)
                                                                                                                  tAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTGCACct          <  2:488932/133‑1 (MQ=255)
                                                                                                                  tAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTGCACct          <  1:158788/133‑1 (MQ=255)
                                                                                                                  tAAATCAAACAATGCGCAATTGATGGTTGCTTTCGCATCTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGACAATCATATTGCACct          >  1:256230/1‑133 (MQ=255)
                                                                                                                          aCAATGCGCAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTGCACCTCTATCTAg  >  2:526700/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                             
CATACAAGGTGAATCCGCTGATGTGATGCCTAAAGCAGCATCCAGGTTAAGTTTGAGCAAGCAAGATGATAAATTTATGCTAGCCTCAAGTGTTACTGACTCTCAAATAAAATTTAAATCAAACAATGCACAATTGATGGTTGCTTTCGCATTTATGCCAATAACAACGAATGGTATTTTACATGATTATACATACGATATAATAATAAATGATAAAAAATATAAAATTGAAAATCATGTTGCACCTCTATCTAG  >  NC_000913/2774196‑2774450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑