Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,274,742	T→A	11.9%	intergenic (+81/+18)	ssb → / ← yjcB	single‑stranded DNA‑binding protein/putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,274,742	0	T	A	11.9%	122.9 / 5.5	41	intergenic (+81/+18)	ssb/yjcB	single‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base A (2/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGG  >  NC_000913/4274596‑4274890                                                                                                                                                   |                                                                                                                                                     cagTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAt                                                                                                                                                    <  1:360168/149‑1 (MQ=255)         tcaggCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCtt                                                                                                                                            <  4:345837/146‑1 (MQ=255)         tCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCtt                                                                                                                                            >  1:110116/1‑149 (MQ=255)                      tCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGt                                                                                                                               >  4:333744/1‑149 (MQ=255)                         aaCGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaa                                                                                                                            >  2:32622/1‑149 (MQ=255)                            gAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAAt                                                                                                                         >  2:470459/1‑149 (MQ=255)                                                tgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTc                                                                                                                                           <  4:467289/111‑1 (MQ=255)                                                tgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTc                                                                                                                                           >  3:467289/1‑111 (MQ=255)                                                 gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCACTGTGGCGgtgt                                                                                                    >  1:235360/1‑149 (MQ=255)                                                     aCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                <  1:8654/149‑1 (MQ=255)                                                       aTTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                              <  1:470459/149‑1 (MQ=255)                                                        ttCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTc                                                                                             >  3:48841/1‑149 (MQ=255)                                                           cGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtc                                                                                          <  3:333744/149‑1 (MQ=255)                                                            gTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtct                                                                                         >  4:230628/1‑149 (MQ=255)                                                             ttCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                        <  2:201609/149‑1 (MQ=255)                                                               cTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtctt                                                                                      <  2:235360/149‑1 (MQ=255)                                                                  aTTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTa                                                                                   <  2:110116/149‑1 (MQ=255)                                                                         aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                            >  4:382649/1‑149 (MQ=255)                                                                         aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                            >  1:80005/1‑149 (MQ=255)                                                                          ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAcc                                                                           >  3:428162/1‑149 (MQ=255)                                                                              aaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                        >  1:141046/1‑116 (MQ=255)                                                                              aaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                        <  2:141046/116‑1 (MQ=255)                                                                                     ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttct                                                                                    >  4:370254/1‑129 (MQ=255)                                                                                     ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttct                                                                                    <  3:370254/129‑1 (MQ=255)                                                                                     ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                <  4:48841/149‑1 (MQ=255)                                                                                          tatTGTTTTAAGGTGGATGATTAAAGGATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGcc                                                           >  2:170626/1‑149 (MQ=255)                                                                                           atTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTAGTGCCa                                                          >  2:2413/1‑149 (MQ=255)                                                                                                           tgatTAAAGGATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                          <  1:170626/149‑1 (MQ=255)                                                                                                                 aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                               >  3:356006/1‑90 (MQ=255)                                                                                                                 aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                               <  4:356006/90‑1 (MQ=255)                                                                                                                 aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                          <  2:228866/143‑1 (MQ=255)                                                                                                                 aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                          >  1:228866/1‑143 (MQ=255)                                                                                                                     aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGc                                                            >  3:40428/1‑121 (MQ=255)                                                                                                                     aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGc                                                            <  4:40428/121‑1 (MQ=255)                                                                                                                      tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                               <  3:230628/149‑1 (MQ=255)                                                                                                                                 tAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAcc                    >  3:265810/1‑149 (MQ=255)                                                                                                                                  aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                    <  1:305464/132‑1 (MQ=255)                                                                                                                                  aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                    >  2:305464/1‑132 (MQ=255)                                                                                                                                      aagaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAtt               <  1:32622/149‑1 (MQ=255)                                                                                                                                         aagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                               <  3:189460/130‑1 (MQ=255)                                                                                                                                         aagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGGTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                               >  4:189460/1‑130 (MQ=255)                                                                                                                                                   tATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATTATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAgg  >  1:361433/1‑149 (MQ=255)                                                                                                                                                   tATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAgg  >  1:88559/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGG  >  NC_000913/4274596‑4274890

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