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

CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑CCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGA  >  NC_000913/2767367‑2767629
                                                                                                                                             |                                                                                                                                    
ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCCC‑CCC‑CT‑CCCAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑c‑c‑cc                                                                                                                               <  2:1158622/139‑1 (MQ=255)
          cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTaaa                                                                                                                      <  1:821573/139‑1 (MQ=255)
                     tCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑AAC‑CCC‑CTATAAAACCTCCTTCATAGGGTAAACTTTTTTACATTTGCACCCAACTCCAAAT‑T‑T‑‑CTC‑ATCCCTGTAAACCTCTCcataa                                                                                                           >  2:4820/1‑139 (MQ=255)
                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAA‑ACC‑CCC‑CT‑AACAACGTCCTGCATAGAGAACACGTTTGCCACTGTGCAACATAGACTAAAT‑T‑T‑‑C‑CTCCCCCTTTAAAATTCTCAATAATg                                                                                                         >  2:814009/1‑139 (MQ=255)
                                           aCTTGTCCAATATGATCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATc                                                                                     <  1:814009/139‑1 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACTTCCTGCATAGGGTACACGTTTTCACCTGTGCACCCACGACTAAATAT‑C‑‑C‑C‑CCCCCTTTAAACTTCTCATTATTTTCACCAATTTTTTGGACATCTTatg                                                                                >  2:799778/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AA‑CCCACGC‑CT‑CACAACGCCCTGCCTAGAGAACACGTGTTCCCCTGCGGACGCACGTCTAAAT‑T‑TC‑C‑C‑CCCCCTTTCAGCTTCTCAATAATGTCACCAATGTGCTGGAAATCCTatg                                                                                >  1:68180/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AA‑CCC‑CCCGCT‑CACAACATCCTGCAAAGGGTACACGATCTCCCCTGTGCACCCACGACTAAAT‑TAA‑‑A‑C‑ACTCCGTTAAACATCTCAATAATGTCACCAAACTTCTggccagcccatg                                                                                >  2:1094166/1‑129 (MQ=255)
                                                         aTCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                       <  1:100072/139‑1 (MQ=255)
                                                          tCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTTCATAGGGAACACGTTTTCCACTCTGCACCCACGACTCAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGACAACCTATGACGCCttttc                                                                      >  2:781852/1‑137 (MQ=255)
                                                          tCAAATAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCCATTTTCTGGAAATCCTATGATGCCTTTcc                                                                      >  2:116170/1‑139 (MQ=255)
                                                            aaaTAGCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                                                    <  1:4820/139‑1 (MQ=255)
                                                                aGCATT‑AACCCC‑CCC‑CT‑CACGTCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCTAAGTCTACAT‑T‑T‑‑C‑CGCCCCCTTTAATCTTCCTAATAAGGTCAGCAATTTTCTGGAAATCCTCTGATGACTTTCCCGGTaa                                                                <  2:20993/139‑1 (MQ=255)
                                                                 gCATT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGGTTTCCCCTGTGCACCCACGACTAAATGT‑C‑‑C‑C‑CCCCCTTTAAACCTCTCAATGATGGCCCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                                               >  2:902331/1‑139 (MQ=255)
                                                                   aTT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCAcc                                                                                                    >  2:1276239/1‑100 (MQ=255)
                                                                   aTT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCAcc                                                                                                    <  1:1276239/100‑1 (MQ=255)
                                                                   aTT‑AACCCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             <  1:1067383/139‑1 (MQ=255)
                                                                     ttaacccc‑ccc‑cT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCCC‑C‑CCCCCTTTAAACTTCTCTATAATGTTACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaagg                                                             >  1:1161682/6‑137 (MQ=255)
                                                                                        aaCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaa                                            <  1:1073953/139‑1 (MQ=255)
                                                                                            tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt                                        >  1:222/1‑139 (MQ=255)
                                                                                             ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt                                       >  2:898013/1‑139 (MQ=255)
                                                                                                 cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                   <  2:222/139‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTaaa                                                                                                                      >  1:866352/1‑51 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTaaa                                                                                                                      >  1:914678/1‑51 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTaaa                                                                                                                      <  2:914678/51‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTaaa                                                                                                                      <  2:866352/51‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCt                              >  1:1029562/1‑139 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTa                             >  2:1010001/1‑139 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                      <  1:226790/78‑1 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                      >  2:226790/1‑78 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             >  1:1277234/1‑103 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                             <  2:1277234/103‑1 (MQ=255)
                                                                                                               tttCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCAtt                     >  1:488713/1‑139 (MQ=255)
                                                                                                                ttCCCCTGTGCCCCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc                                      <  1:667515/121‑1 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                                                          >  1:979533/1‑85 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc                                      >  2:667515/1‑121 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCCTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                                                          <  2:979533/85‑1 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CCCCCcttta                                                                                                                       >  1:811782/1‑37 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CCCCCcttta                                                                                                                       <  2:811782/39‑3 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                                                             >  1:1021607/1‑75 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                                                             <  2:1021607/75‑1 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAAt             >  1:881478/1‑139 (MQ=255)
                                                                                                                        tgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                              >  2:543771/1‑105 (MQ=255)
                                                                                                                        tgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg                                              <  1:543771/105‑1 (MQ=255)
                                                                                                                        tgCACCCACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATc            <  2:881478/139‑1 (MQ=255)
                                                                                                                         tgcgttgtttaagaaT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCt           >  1:752447/14‑139 (MQ=255)
                                                                                                                         tgcgttgtttaagaaT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCt           >  1:357325/14‑139 (MQ=255)
                                                                                                                             ccACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                    <  1:876362/110‑1 (MQ=255)
                                                                                                                             ccACGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                    >  2:876362/1‑110 (MQ=255)
                                                                                                                               aCGACTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat     <  1:799778/139‑1 (MQ=255)
                                                                                                                                  aCTAAATGT‑C‑‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCAGTGAAAAAATCTTCagataga  <  2:68180/139‑1 (MQ=255)
                                                                                                                                  aCTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  <  1:902331/139‑1 (MQ=255)
                                                                                                                                  aCTAAAT‑T‑TC‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  >  2:117172/1‑139 (MQ=255)
                                                                                                                                   cTAAAT‑T‑TCCC‑C‑CCCCCTTTAAACTTCGCAATAATGTCACCTATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga  <  2:1161682/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                    
CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑CCC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGA  >  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: ‑