Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	272,330:1	+G	8.2%	intergenic (+141/‑250)	insN → / → eyeA	pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	272,330	1	.	G	8.2%	196.5 / 13.6	49	intergenic (+141/‑250)	insN/eyeA	pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge
Reads supporting (aligned to +/- strand):  ref base . (29/16);  new base G (2/2);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-02

GACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACGG  >  NC_000913/272183‑272476                                                                                                                                                     |                                                                                                                                                   gACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGc                                                                                                                                                     <  4:356194/148‑1 (MQ=255)     tataGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCa                                                                                                                                                >  3:273254/1‑148 (MQ=255)         gtgtGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGccc                                                                                                                                           <  1:271990/149‑1 (MQ=255)            tgCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAGCCTGAGCGGCAGTGACCTGCGGGCATGccgc                                                                                                                                          <  3:174730/148‑3 (MQ=255)                cGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑TGCATGCCCGCAGGGt                                                                                                                                    <  2:19568/149‑1 (MQ=255)                               gcATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTg                                                                                                                      >  4:130068/1‑148 (MQ=255)                                   ggATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAAc                                                                                                                 >  1:193500/1‑149 (MQ=255)                                     aTGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                               >  2:48738/1‑149 (MQ=255)                                      tGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                                                               >  3:153508/1‑148 (MQ=255)                                         ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAggg                                                                                                                                     <  2:362919/123‑1 (MQ=255)                                         ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAggg                                                                                                                                     >  1:362919/1‑123 (MQ=255)                                          tAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAt                                                                                                           >  3:174496/1‑148 (MQ=255)                                            aaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACggg                                                                                                              >  1:101124/1‑143 (MQ=255)                                            aaTGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACggg                                                                                                              <  2:101124/143‑1 (MQ=255)                                                        gtctgtctGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgca                                                                                                                                         <  3:21049/105‑4 (MQ=255)                                                        gtctgtctGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgca                                                                                                                                         >  4:21049/1‑102 (MQ=255)                                                           tgtctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                          >  3:193120/1‑148 (MQ=255)                                                              ctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGAGCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGcc                                                                                       >  4:229201/1‑148 (MQ=255)                                                                 gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgag                                                                                    <  4:153508/148‑1 (MQ=255)                                                                          aCCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCa                                                                          <  1:261922/149‑1 (MQ=255)                                                                            cGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAg                                                                        <  2:293867/149‑1 (MQ=255)                                                                                   tCACNATGAGAAACCCAACAATCCCACCGGGTNTGACGGTGGAGAACNTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  1:225797/1‑149 (MQ=255)                                                                                   tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  1:22801/1‑149 (MQ=255)                                                                                    cACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  4:362100/1‑148 (MQ=255)                                                                                    cACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  3:194852/1‑148 (MQ=255)                                                                                    cACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  4:107997/1‑148 (MQ=255)                                                                                        aTGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTg                                                             >  3:10294/1‑148 (MQ=255)                                                                                                 ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                   >  1:39863/1‑149 (MQ=255)                                                                                                  caacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaa                                                   >  4:120114/1‑148 (MQ=255)                                                                                                   aacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaaga                                                 >  2:257276/1‑149 (MQ=255)                                                                                                           cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAc                                         >  1:163059/1‑149 (MQ=255)                                                                                                            aCCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                  <  4:215928/91‑1 (MQ=255)                                                                                                            aCCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgag                                                                                                  >  3:215928/1‑91 (MQ=255)                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCg                                      <  4:10294/148‑1 (MQ=255)                                                                                                               gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGGAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGt                                     <  2:264290/149‑1 (MQ=255)                                                                                                                 gtgtGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 <  4:309879/119‑1 (MQ=255)                                                                                                                 gtgtGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGt                                                                 >  3:309879/1‑119 (MQ=255)                                                                                                                         gTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGcc                            >  4:244030/1‑148 (MQ=255)                                                                                                                           ggAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCctct                         >  2:190904/1‑149 (MQ=255)                                                                                                                               aaCCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTgg                     <  1:190904/149‑1 (MQ=255)                                                                                                                                 ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGt                                                                                                                               >  2:395682/1‑41 (MQ=255)                                                                                                                                 ccTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGt                                                                                                                               <  1:395682/41‑1 (MQ=255)                                                                                                                                  cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGtt                  >  2:244568/1‑149 (MQ=255)                                                                                                                                  cTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGt                   <  3:130068/148‑1 (MQ=255)                                                                                                                                         gCAGTGACCTGC‑GGCATGCCCTCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTACGCCGAGCTGATCCTCAAGGTTAAGTACTGAAAGGTTGTAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAAct           >  1:292149/1‑148 (MQ=255)                                                                                                                                          cAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATcc                               >  2:155492/1‑128 (MQ=255)                                                                                                                                          cAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATcc                               <  1:155492/128‑1 (MQ=255)                                                                                                                                              acccTGCGGGCATG‑CCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacc      >  1:344624/3‑149 (MQ=255)                                                                                                                                               aCCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca     >  1:325776/1‑149 (MQ=255)                                                                                                                                                  tGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACgg  <  2:163059/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGC‑GGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCACGG  >  NC_000913/272183‑272476

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: ‑