Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,507 (C)7→8 intergenic (+152/‑203) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5001.C85.7% 143.5 / 15.3 49intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base C (22/20);  total (26/23)
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.00e+00

CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGA  >  NC_000913/2767367‑2767629
                                                                                                                                             |                                                                                                                                    
ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCCCAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccc                                                                                                                               <  2:1158622/139‑1 (MQ=255)
          cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                                                      <  1:821573/139‑1 (MQ=255)
                     tCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAAACCCCCTATAAAACCTCCTTCATAGGGTAAACTTTTTTACATTTGCACCCAACTCCAAATTT‑‑CTCATCCCTGTAAACCTCTCcataa                                                                                                           >  2:4820/1‑139 (MQ=255)
                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAAACCCCCCTAACAACGTCCTGCATAGAGAACACGTTTGCCACTGTGCAACATAGACTAAATTT‑‑CCTCCCCCTTTAAAATTCTCAATAATg                                                                                                         >  2:814009/1‑139 (MQ=255)
                                           aCTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATc                                                                                     <  1:814009/139‑1 (MQ=255)
                                                tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACTTCCTGCATAGGGTACACGTTTTCACCTGTGCACCCACGACTAAATATC‑‑CCCCCCCTTTAAACTTCTCATTATTTTCACCAATTTTTTGGACATCTTatg                                                                                >  2:799778/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATTAACCCACGCCTCACAACGCCCTGCCTAGAGAACACGTGTTCCCCTGCGGACGCACGTCTAAATTTCCCCCCCCTTTCAGCTTCTCAATAATGTCACCAATGTGCTGGAAATCCTatg                                                                                >  1:68180/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATTAACCCCCCGCTCACAACATCCTGCAAAGGGTACACGATCTCCCCTGTGCACCCACGACTAAATTAA‑‑ACACTCCGTTAAACATCTCAATAATGTCACCAAACTTCTggccagcccatg                                                                                >  2:1094166/1‑129 (MQ=255)
                                                         aTCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                       <  1:100072/139‑1 (MQ=255)
                                                          tCAAATAGCATTAACCCCCCCCTCACAACGTCCTTCATAGGGAACACGTTTTCCACTCTGCACCCACGACTCAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGACAACCTATGACGCCttttc                                                                      >  2:781852/1‑137 (MQ=255)
                                                          tCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCCATTTTCTGGAAATCCTATGATGCCTTTcc                                                                      >  2:116170/1‑139 (MQ=255)
                                                            aaaTAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                                                    <  1:4820/139‑1 (MQ=255)
                                                                aGCATTAACCCCCCCCTCACGTCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCTAAGTCTACATTT‑‑CCGCCCCCTTTAATCTTCCTAATAAGGTCAGCAATTTTCTGGAAATCCTCTGATGACTTTCCCGGTaa                                                                <  2:20993/139‑1 (MQ=255)
                                                                 gCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGGTTTCCCCTGTGCACCCACGACTAAATGTC‑‑CCCCCCCTTTAAACCTCTCAATGATGGCCCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                                               >  2:902331/1‑139 (MQ=255)
                                                                   aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCAcc                                                                                                    >  2:1276239/1‑100 (MQ=255)
                                                                   aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCAcc                                                                                                    <  1:1276239/100‑1 (MQ=255)
                                                                   aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             <  1:1067383/139‑1 (MQ=255)
                                                                     ttaaccccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTTAAACTTCTCTATAATGTTACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaagg                                                             >  1:1161682/6‑137 (MQ=255)
                                                                                        aaCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaa                                            <  1:1073953/139‑1 (MQ=255)
                                                                                            tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt                                        >  1:222/1‑139 (MQ=255)
                                                                                             ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt                                       >  2:898013/1‑139 (MQ=255)
                                                                                                 cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                   <  2:222/139‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                                                      >  1:866352/1‑51 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                                                      >  1:914678/1‑51 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                                                      <  2:914678/51‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                                                                                      <  2:866352/51‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCt                              >  1:1029562/1‑139 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTa                             >  2:1010001/1‑139 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                      <  1:226790/78‑1 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                      >  2:226790/1‑78 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             >  1:1277234/1‑103 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             <  2:1277234/103‑1 (MQ=255)
                                                                                                               tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCAtt                     >  1:488713/1‑139 (MQ=255)
                                                                                                                ttCCCCTGTGCCCCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc                                      <  1:667515/121‑1 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                                                          >  1:979533/1‑85 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc                                      >  2:667515/1‑121 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                                                          <  2:979533/85‑1 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                                                       >  1:811782/1‑37 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                                                       <  2:811782/39‑3 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                                                             >  1:1021607/1‑75 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                                                             <  2:1021607/75‑1 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAAt             >  1:881478/1‑139 (MQ=255)
                                                                                                                        tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                              >  2:543771/1‑105 (MQ=255)
                                                                                                                        tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                              <  1:543771/105‑1 (MQ=255)
                                                                                                                        tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATc            <  2:881478/139‑1 (MQ=255)
                                                                                                                         tgcgttgtttaagaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCt           >  1:752447/14‑139 (MQ=255)
                                                                                                                         tgcgttgtttaagaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCt           >  1:357325/14‑139 (MQ=255)
                                                                                                                             ccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                    <  1:876362/110‑1 (MQ=255)
                                                                                                                             ccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                    >  2:876362/1‑110 (MQ=255)
                                                                                                                               aCGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat     <  1:799778/139‑1 (MQ=255)
                                                                                                                                  aCTAAATGTC‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCAGTGAAAAAATCTTCagataga  <  2:68180/139‑1 (MQ=255)
                                                                                                                                  aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  <  1:902331/139‑1 (MQ=255)
                                                                                                                                  aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  >  2:117172/1‑139 (MQ=255)
                                                                                                                                   cTAAATTTCCCCCCCCCTTTAAACTTCGCAATAATGTCACCTATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  <  2:1161682/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                    
CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGA  >  NC_000913/2767367‑2767629

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