Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 265,418 A→G R42R (CGT→CGC ykfF ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,4180AG98.1% 178.4 / ‑4.2 52R42R (CGT→CGCykfFCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/30);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265328‑265514
                                                                                          |                                                                                                
gATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccacc                                                                                         <  1:740631/100‑1 (MQ=255)
  tGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAg                                                                                       <  2:206782/100‑1 (MQ=255)
   gCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGa                                                                                      <  1:533752/100‑1 (MQ=255)
    ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                     <  1:1031208/100‑1 (MQ=255)
    ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                     <  1:214062/100‑1 (MQ=255)
    ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                     >  2:1280481/1‑100 (MQ=255)
          gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGt                                                                               <  1:573204/100‑1 (MQ=255)
            gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGa                                                                             >  1:903988/1‑100 (MQ=255)
             cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                            <  2:1127946/100‑1 (MQ=255)
                      cGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgt                                                                   >  1:1414526/1‑100 (MQ=255)
                            aggccTTCACCGGCATCCGGCTCAGAATTCCATGCTCGCCAGATCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                             <  2:207734/97‑1 (MQ=21)
                             gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGCCGGAAGTGACTGCCCTGGTCGACTTc                                                            >  1:25611/1‑100 (MQ=255)
                             gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTc                                                            >  1:1060435/1‑100 (MQ=255)
                             gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTCCCCTGGTCGTCTTc                                                            >  2:855042/1‑100 (MQ=255)
                               ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCg                                                           >  2:1305397/1‑99 (MQ=255)
                                  tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                       <  1:4293/100‑1 (MQ=255)
                                  tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                       >  1:701044/1‑100 (MQ=255)
                                   cACCGGCATCCGGCTCAAAGTTCCAGGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                      <  2:1420907/100‑1 (MQ=255)
                                     ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGTCCTGGTCGTCTTCGAGTGTGa                                                    <  1:1111892/100‑1 (MQ=255)
                                      cGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGtgtg                                                     <  2:1231309/98‑1 (MQ=255)
                                       ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                  >  2:100308/1‑100 (MQ=255)
                                        gCATCCGGCTCAAAGTTCCATGCCCGCCAGACCCTCCGGCCTTCAGTATCGCGAACCACCCGACGGGAGTGACTGCCCTGGTCGTCTTCGAGAGTGat                                                   >  2:46039/1‑98 (MQ=37)
                                           tCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                              <  2:623874/100‑1 (MQ=255)
                                           tCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                              <  1:100308/100‑1 (MQ=255)
                                           tCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                              >  2:803756/1‑100 (MQ=255)
                                            ccGGCTCGAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCt                                             >  1:1405569/1‑100 (MQ=255)
                                             cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                            <  1:698508/100‑1 (MQ=255)
                                             cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                            <  1:1244772/100‑1 (MQ=255)
                                             cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                            <  1:218109/100‑1 (MQ=255)
                                                  cgaaGTTCCATGCCCGCCAGACCATCGTGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACg                                        <  2:1132579/97‑1 (MQ=37)
                                                  cAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGt                                       >  1:243618/1‑100 (MQ=255)
                                                   aaaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTg                                      <  1:717946/100‑1 (MQ=255)
                                                   aaaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTg                                      <  1:1281201/100‑1 (MQ=255)
                                                    aaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                     <  2:1284925/100‑1 (MQ=255)
                                                       ttCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                  >  2:1177915/1‑100 (MQ=255)
                                                         ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATGTTGCTGTACGTGGTAGTg                                >  2:644096/1‑100 (MQ=255)
                                                         ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACgtggtggtg                                >  2:1433957/1‑100 (MQ=255)
                                                            tGCCCGCCAGCCCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCGTCCAGTGTGATATTGCTGTACGTGGTAGTGAcc                             <  2:1040002/100‑1 (MQ=255)
                                                            tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                             <  1:256442/100‑1 (MQ=255)
                                                            tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                             <  1:631656/100‑1 (MQ=255)
                                                            tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                             <  2:484897/100‑1 (MQ=255)
                                                             gCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcg                            >  2:740967/1‑100 (MQ=255)
                                                               ccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTTCGTGGTAGTGACcgct                          >  2:926051/1‑100 (MQ=255)
                                                                     gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                    <  1:1503840/100‑1 (MQ=255)
                                                                              gCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCgg           <  1:904099/100‑1 (MQ=255)
                                                                               ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg          <  2:1414518/100‑1 (MQ=255)
                                                                                  tCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGa       >  1:558120/1‑100 (MQ=255)
                                                                                   cAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGaa      >  2:538067/1‑100 (MQ=255)
                                                                                     gTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAg    <  2:20714/100‑1 (MQ=255)
                                                                                     gTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAg    <  2:747254/100‑1 (MQ=255)
                                                                                     gTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAg    >  2:886793/1‑100 (MQ=255)
                                                                                       aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc  >  1:997459/1‑100 (MQ=255)
                                                                                          |                                                                                                
GATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265328‑265514

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: