Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,564,525 G→A M1M (ATG→ATA) † yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,5250GA97.1% 232.9 / ‑2.1 69M1M (ATG→ATA) †yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (39/28);  total (40/29)
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.97e-01

GATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564430‑2564620
                                                                                               |                                                                                               
gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCATTGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  2:1543250/1‑100 (MQ=255)
gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGtt                                                                                             >  2:917881/1‑100 (MQ=255)
gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  1:760036/1‑100 (MQ=255)
gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  1:79104/1‑100 (MQ=255)
 aTCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTa                                                                                            <  1:242888/100‑1 (MQ=255)
  tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTAc                                                                                           <  1:675524/100‑1 (MQ=255)
   cTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACg                                                                                          >  2:675007/1‑100 (MQ=255)
    tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         >  1:190104/1‑100 (MQ=255)
    tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         <  1:415607/100‑1 (MQ=255)
    tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         <  1:1392320/100‑1 (MQ=255)
     gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                        <  1:1493198/100‑1 (MQ=255)
         cAATATTCCCTTCCGTTACGCTTACCCTCCACCCATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    >  1:512350/1‑100 (MQ=255)
         cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    >  2:448935/1‑100 (MQ=255)
             atTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAt                                                                                <  2:1439594/100‑1 (MQ=255)
               tCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAttt                                                                              <  2:714927/100‑1 (MQ=255)
                cccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTg                                                                             <  2:517927/100‑1 (MQ=255)
                 ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            >  1:1538435/1‑100 (MQ=255)
                    tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                         >  2:554353/1‑100 (MQ=255)
                        ttACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATCAAAATCATAAGTTACGAAATCAAAATTTGCGACATcc                                                                     <  2:848069/100‑1 (MQ=255)
                         tACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCt                                                                    <  2:1284297/100‑1 (MQ=255)
                          aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                   >  2:1248622/1‑100 (MQ=255)
                          aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                   >  2:1547702/1‑100 (MQ=255)
                           cGCTTACCCCCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  2:1207991/100‑1 (MQ=255)
                            gCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaa                                                                 >  1:512755/1‑100 (MQ=255)
                             cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                >  1:730540/1‑100 (MQ=255)
                               cgatcTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAgg                                                              <  1:1285729/96‑1 (MQ=255)
                                aCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgc                                                             >  2:367738/1‑100 (MQ=255)
                                 cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:905241/1‑100 (MQ=255)
                                 cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:712584/1‑100 (MQ=255)
                                 cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:492898/1‑100 (MQ=255)
                                   cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          >  2:1147828/1‑100 (MQ=255)
                                     ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAggagatcg                                                        >  2:1285724/1‑94 (MQ=255)
                                     ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTg                                                         <  2:1456429/99‑1 (MQ=255)
                                     ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                        <  2:164470/100‑1 (MQ=255)
                                        ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     >  2:800349/1‑100 (MQ=255)
                                         caATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                                    <  2:760034/100‑1 (MQ=255)
                                            tACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGaa                                                 >  2:1303698/1‑100 (MQ=255)
                                              cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               <  2:1280818/100‑1 (MQ=255)
                                              cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               >  1:1453437/1‑100 (MQ=255)
                                                cgTGTTAACCTCCCAATGGGTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             >  1:433028/1‑100 (MQ=255)
                                                   gtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAat                                          <  1:930/100‑1 (MQ=255)
                                                      aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       >  1:644044/1‑100 (MQ=255)
                                                      aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       <  2:1357566/100‑1 (MQ=255)
                                                         cTCCCAATGGATTCTATGAGATCGGAGATAAAAATCATAAGTTACGCAATCAAAATTTGCGACATCCTCAAACGCGCTGCGATGGAAGGACATTATAAg                                     >  2:1348213/1‑99 (MQ=255)
                                                             cAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcgc                                >  2:56472/1‑100 (MQ=255)
                                                             cAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcgc                                >  1:1392160/1‑100 (MQ=255)
                                                             cAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcg                                 >  2:367990/1‑99 (MQ=255)
                                                             cAATGGATTCTATGAGATGGCAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcgc                                >  2:965657/1‑100 (MQ=255)
                                                              aaTGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCt                               >  1:1161535/1‑100 (MQ=255)
                                                               aTGGATTCTATTAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAGTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTc                              <  2:997703/100‑1 (MQ=255)
                                                                 ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            <  2:644042/100‑1 (MQ=255)
                                                                        aTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGGCAATATAAGGGCGCTCAACGTgggg                     >  1:1480695/1‑100 (MQ=255)
                                                                         tGAGATGGGAGATAAAAATCATAAGTTACGATATCAAGATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  2:1065696/100‑1 (MQ=255)
                                                                         tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgggg                     >  2:1084223/1‑99 (MQ=255)
                                                                         tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  2:1501272/100‑1 (MQ=255)
                                                                          gagaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGcc                   <  2:427432/100‑1 (MQ=255)
                                                                          gagaTGGGAGATAAAAATCATAAGCTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCCGCGATGGAAGGACAATATAAGGGCGCGCAACGTGGGGcc                   >  2:1473891/1‑100 (MQ=255)
                                                                           agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  <  1:1341103/100‑1 (MQ=255)
                                                                                ggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAACTGc             >  1:120730/1‑100 (MQ=255)
                                                                                  agaTAAAAATCATAAGTTACGAAATCAAAATTTGAGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  2:717759/1‑97 (MQ=255)
                                                                                   gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  1:1309811/100‑1 (MQ=255)
                                                                                    aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  1:3704/100‑1 (MQ=255)
                                                                                    aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGTGCTCAACGTGGGGCCAAATGCgagg         <  1:186938/100‑1 (MQ=255)
                                                                                    aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCACTCAACGTGGGGCCAAATGCgagg         <  2:57406/100‑1 (MQ=255)
                                                                                     tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:873476/1‑100 (MQ=255)
                                                                                      aaaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgaggag       >  2:933430/1‑100 (MQ=255)
                                                                                      aaaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgaggag       >  2:1155707/1‑100 (MQ=255)
                                                                                        aaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGAt     >  1:1286073/1‑100 (MQ=255)
                                                                                           tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  1:448936/100‑1 (MQ=255)
                                                                                               |                                                                                               
GATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564430‑2564620

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: