Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,365,608	T→G	intergenic (+57/+94)	yhcE → / ← insH1	pseudogene fragment; interrupted by IS5/IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,365,608	0	T	G	97.4%	132.7 / NA	38	intergenic (+57/+94)	yhcE/insH1	pseudogene fragment; interrupted by IS5/IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/19);  total (18/20)
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

GCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGT  >  NC_000913/3365461‑3365699                                                                                                                                                    |                                                                                            gCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATgtct                                                                                          >  1:681203/1‑151 (MQ=255)   cAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATgtcttc                                                                                        >  1:139339/1‑151 (MQ=255)          gCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTg                                                                                 >  1:398591/1‑151 (MQ=255)            aTTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAg                                                                               <  1:936092/151‑1 (MQ=255)            aTTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAg                                                                               <  2:328618/151‑1 (MQ=255)                  aGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGAtt                                                                         <  2:723209/151‑1 (MQ=255)                       gtaACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTtc                                                                    >  1:537013/1‑151 (MQ=255)                         aaCAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTtct                                                                   <  1:294512/150‑1 (MQ=255)                         aaCAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTtct                                                                   <  1:334079/150‑1 (MQ=255)                            aaTCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                               >  1:312565/1‑151 (MQ=255)                                tcGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                            >  2:79696/1‑150 (MQ=255)                                 cGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTa                                                             <  2:217893/148‑1 (MQ=255)                                 cGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTa                                                             <  1:158456/148‑1 (MQ=255)                                  gACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATg                                                         <  1:363730/151‑1 (MQ=255)                                   aCCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                        >  2:231161/1‑151 (MQ=255)                                   aCCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                        <  2:359190/151‑1 (MQ=255)                                    ccTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                       >  2:301923/1‑151 (MQ=255)                                    ccTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                       <  1:858387/151‑1 (MQ=255)                                     cTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                      >  1:41138/1‑151 (MQ=255)                                      tGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                     >  2:636626/1‑151 (MQ=255)                                          gTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                 >  1:150729/1‑151 (MQ=255)                                          gTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                 <  1:164439/151‑1 (MQ=255)                                           ttAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                >  1:742460/1‑151 (MQ=255)                                             aGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                               <  2:632917/150‑1 (MQ=255)                                             aGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                              <  2:652700/151‑1 (MQ=255)                                              gCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                             >  1:523572/1‑151 (MQ=255)                                                cAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCt                                           >  1:751756/1‑151 (MQ=255)                                                    ggCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAttt                                       <  2:138368/151‑1 (MQ=255)                                                                   ccGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                         <  1:389545/150‑1 (MQ=255)                                                                   ccGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                        >  1:249923/1‑151 (MQ=255)                                                                       gcagcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAtt                    <  2:103872/151‑1 (MQ=255)                                                                        cagcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                   >  1:863119/1‑151 (MQ=255)                                                                          gcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa                 <  2:737561/151‑1 (MQ=255)                                                                             catgcaCAGCATTCATTAGCGAAGGTTCGAAGAAGCGGTGAAATTCTTCTCGGCTGGCTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTGTTTCTCACATAAATGCCTGGAAGCAGCCTATTTAGACCGTTTCTTCGCCATTTa                  <  1:391682/142‑1 (MQ=11)                                                                              aGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg             <  1:223439/151‑1 (MQ=255)                                                                                    aGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc       >  2:335860/1‑151 (MQ=25)                                                                                      cATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc     >  1:306831/1‑151 (MQ=18)                                                                                       aTCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc     >  1:233039/1‑150 (MQ=18)                                                                                       aTCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc     >  1:175695/1‑150 (MQ=18)                                                                                         cAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGt  <  1:698107/151‑1 (MQ=14)                                                                                                                                                    |                                                                                            GCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGT  >  NC_000913/3365461‑3365699

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