Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,757,662 (T)8→7 coding (939/957 nt) yfjH ← CP4‑57 prophage; protein YfjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,6550T.97.2% 139.8 / ‑1.8 36coding (946/957 nt)yfjHCP4‑57 prophage; protein YfjH
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (18/17);  total (19/17)
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.61e-01

TATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTAT  >  NC_000913/2757524‑2757790
                                                                                                                                   |                                                                                                                                       
tATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttct                                                                                                                                 >  2:328523/1‑138 (MQ=255)
tATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttct                                                                                                                                 <  2:556913/140‑3 (MQ=255)
        cgcAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAg                                                                                                                        >  2:616126/1‑140 (MQ=255)
          cAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAgcg                                                                                                                      >  1:57645/1‑140 (MQ=255)
                tCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCg                                                                                                                <  2:628329/140‑1 (MQ=255)
                     tCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGc                                                                                                           <  1:328523/140‑1 (MQ=255)
                                  gCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCt                                                                                              >  2:668843/1‑140 (MQ=255)
                                     ccTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCa                                                                                           <  1:637684/140‑1 (MQ=255)
                                      cTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAt                                                                                          <  2:57645/140‑1 (MQ=255)
                                              aCGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                                                                  <  1:604747/140‑1 (MQ=255)
                                                  aTCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGAcc                                                                              >  1:329063/1‑140 (MQ=255)
                                                         gTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                                                                       >  1:136578/1‑140 (MQ=255)
                                                         gTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                                                                       >  1:530327/1‑140 (MQ=255)
                                                         gTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                                                                       <  1:668843/140‑1 (MQ=255)
                                                            aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                                                                    >  2:313856/1‑140 (MQ=255)
                                                              taaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaaaa                                                                  <  1:425873/140‑1 (MQ=255)
                                                                   aGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTa                                                             <  1:479414/140‑1 (MQ=255)
                                                                    gAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAt                                                                            >  1:152163/1‑124 (MQ=255)
                                                                    gAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAt                                                                            <  2:152163/124‑1 (MQ=255)
                                                                    gAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTAt                                                            >  2:547711/1‑140 (MQ=255)
                                                                            aGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCaaa                                                    <  2:599601/140‑1 (MQ=255)
                                                                                  aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCt                                              <  1:313856/140‑1 (MQ=255)
                                                                                      ggCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGAc                                          <  1:682629/140‑1 (MQ=255)
                                                                                       gcgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACt                                         <  2:544734/140‑1 (MQ=255)
                                                                                        cgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTc                                        <  2:136578/140‑1 (MQ=255)
                                                                                             cTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCAttgt                                   <  1:616126/140‑1 (MQ=255)
                                                                                                     tCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTacac                           >  1:644415/1‑140 (MQ=255)
                                                                                                       tataTGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGa                         >  1:551419/1‑140 (MQ=255)
                                                                                                       tataTGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGa                         >  1:298535/1‑140 (MQ=255)
                                                                                                          atGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGatat                      >  2:127538/1‑140 (MQ=255)
                                                                                                                 gtagtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTAt               <  2:230211/140‑1 (MQ=255)
                                                                                                                   agtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTAttt             >  2:557426/1‑140 (MQ=255)
                                                                                                                       attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGat                        >  2:221429/1‑125 (MQ=255)
                                                                                                                       attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGat                        <  1:221429/125‑1 (MQ=255)
                                                                                                                            ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGattatt     >  2:4306/1‑139 (MQ=255)
                                                                                                                            ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGattatt     <  1:4306/139‑1 (MQ=255)
                                                                                                                            ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGAttattt    >  2:599354/1‑140 (MQ=255)
                                                                                                                              aCTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGAttattta   >  2:468747/1‑140 (MQ=255)
                                                                                                                              aCTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGAttatttat  >  2:187844/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTGACTCATTGTTGCTACACGATATTTTTTATTTTGATTATTTAT  >  NC_000913/2757524‑2757790

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