Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,286,793	C→T	pseudogene (2122/2525 nt)	yejO ←	adhesin‑like autotransporter YejO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,286,793	0	C	T	96.2%	78.6 / ‑5.8	26	pseudogene (2122/2525 nt)	yejO	adhesin‑like autotransporter YejO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/13);  minor base A (0/1);  total (12/14)
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

CGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCA  >  NC_000913/2286557‑2287035                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   cGTCATTACCTTTCACCTGGTTAGAACCGGCTTATTCCTGACGCACGGCCGCTTTCATCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACgg                                                                                                                                                                                                                                      <  1:214562/251‑1 (MQ=255)         aCTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGcc                                                                                                                                                                                                                              >  2:197983/1‑250 (MQ=255)              cACCTGGTTAGACTCGGCGTTTTCCTGACGCGCGTCCGCTTGCAGCCAGGGTTCC‑AGCTTCATTCCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCaa                                                                                                                                                                                                                         <  2:220838/251‑1 (MQ=255)                              gctttgtCCTGACGCTCGGCCTTTTTCA‑CCGGGGGTCCGAGCGCCGTACCGTTT‑GCAAGGTCCTTGGTG‑TAGTTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTTTTTCCCACATCCTCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATTTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAAGTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGcc                                                                                                                                                                                                         <  2:139565/245‑1 (MQ=255)                                 tATTCCTGCCGCACGGCCGCTTTCAGCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAACGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGtt                                                                                                                                                                                                      >  2:115819/1‑251 (MQ=255)                                      cTGACGCACGGCCGCTTTCAGCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTAt                                                                                                                                                                                                 >  2:42008/1‑251 (MQ=255)                                                           gccaggctcccagtgtcttccctttttGCA‑GTTCAATG‑TGTTAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATTAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGtt                                                                                                                                                                             <  2:20435/228‑1 (MQ=255)                                                                   tgaaGCGTAGTCCCGTTTTTCA‑GGTCCGTG‑TAGTTGCTTACCGCCGTTCCCGCTTCA‑CGTCGTAATATTCGGTTATTTCCCACATCCGCGCGCATGCCGTTTGATAACTGGTAGTCCTGACCATCTGTGGTAAGGCCTGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCTTCAACCCAACGGAACCCGCTCTCAACATGAGCCCCCGCGCCGTTGCTATTGTAATCGCCCTACGCTGATGCCCCAt                                                                                                                                                                      <  2:208163/247‑1 (MQ=255)                                                                         gtcgtACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCAt                                                                                                                                                               >  1:2831/1‑251 (MQ=255)                                                                                              cATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGACGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGAGTTgg                                                                                                                                           >  1:255805/1‑251 (MQ=255)                                                                                                                            cAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATc                                                                                                               <  1:197983/250‑1 (MQ=255)                                                                                                                                           cGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGATACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATTGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGtt                                                                                                >  2:245648/1‑251 (MQ=255)                                                                                                                                              gTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTg                                                                                             <  1:87305/251‑1 (MQ=255)                                                                                                                                                    cccACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCt                                                                                        <  2:68556/250‑1 (MQ=255)                                                                                                                                                    cccACATCCGCGCGCATGCCGATTGATACCGTCTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGCTAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGCATGGTCTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTc                                                                                       >  1:106027/1‑251 (MQ=255)                                                                                                                                                                     gCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGcc                                                                       <  1:42008/250‑1 (MQ=255)                                                                                                                                                                          ttGATAACGTGTAGCCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAgg                                                                 >  1:106495/1‑251 (MQ=255)                                                                                                                                                                                   tgaaGTCCAGACCATCTGTGGTACAGCCGTTAACGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACAATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTTCCATGGATGGTGTTGGCACAACGGTCAACTTTCACCACCCCATCAACATAGGCACTGTTCTGATGCTCCCCACCGGCATCAGCCCCCAGGGTATAGCt                                                         >  2:150428/4‑250 (MQ=255)                                                                                                                                                                                                    gtgGTAAAGCCGGTAATGGCAAGATAGGGTCTAACACTACCCAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACACGCGCCGTTACTATTATACTCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAAATTTCACCACCCCATCAACATAGGCACCGTTCTGATGCCCCCAACAGGCATAAGCCCCCAGGGTATACCTATCGATATTACCTTTgcc                                       >  2:231402/1‑251 (MQ=255)                                                                                                                                                                                                               ggTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACTCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGAt                               <  1:246875/248‑1 (MQ=255)                                                                                                                                                                                                               ggTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACTCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGAt                               >  2:246875/1‑248 (MQ=255)                                                                                                                                                                                                                gTAAAGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaaa                           <  2:59124/251‑1 (MQ=255)                                                                                                                                                                                                                   aaGGCCAGATAGGGTCTAACACTCCACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAcc                        >  1:148619/1‑251 (MQ=255)                                                                                                                                                                                                                          aaTNGGGTCTAACACTCCACAATTCGTCATCCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCa                 <  2:258987/250‑1 (MQ=255)                                                                                                                                                                                                                                      cacTCCACACTACGTCATCCCAACGGATCCCTCTCTCAACATGAGCACCCGGTCCGTTACTATTTTAATCGCCAAGCGCTGTGGCCCCATTACTCATCTTGCCATGGATGGTTTTGGCAAAACGGTCAACTTGCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCTGCATAAGCCCCCAGATTATGTCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAAAGAGAGTaa     <  1:169619/251‑1 (MQ=255)                                                                                                                                                                                                                                          ccACAATTCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCa  <  1:203712/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   CGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCC‑AGCGTCGTACCGTTTTGCA‑GGTCCATG‑TGATAGCTTACCGCCGTTCCCGCTTCAGCG‑CGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGCCGGTAAAGGCCAGATAGGGTCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCA  >  NC_000913/2286557‑2287035

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