Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,321	Δ1 bp	coding (12/1182 nt)	yjgN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,321	0	A	.	98.6%	305.7 / ‑3.2	69	coding (12/1182 nt)	yjgN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (36/32);  minor base C (1/0);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGT  >  NC_000913/4477254‑4477391                                                                    |                                                                       taAGTTATATCAATGTATTTTTTAA‑TTCCGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTccc                                                              <  2:3317196/76‑1 (MQ=255)  aaGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCa                                                             >  1:2984113/1‑76 (MQ=255)       atatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgt                                                        >  2:4365006/1‑76 (MQ=255)       atatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgt                                                        >  2:4407521/1‑76 (MQ=255)        tatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtt                                                       <  1:2341833/76‑1 (MQ=255)          tCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttg                                                     >  2:2092250/1‑76 (MQ=255)           cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgt                                                    >  2:1782507/1‑76 (MQ=255)           cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATGAATGAAATGGAT‑GTTCCGTCCCATTCgtttgt                                                    <  2:2506631/76‑1 (MQ=255)            aaTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt                                                   >  1:549859/1‑76 (MQ=255)            aaTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt                                                   >  1:1750205/1‑76 (MQ=255)              tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt                                                 >  1:3121426/1‑76 (MQ=255)              tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt                                                 <  1:5810702/76‑1 (MQ=255)               ggATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGttttt                                                <  1:64151/76‑1 (MQ=255)                gATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTc                                               >  1:4293487/1‑76 (MQ=255)                  ttttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTTGTTTGTTTTTCAt                                             <  2:2257573/76‑1 (MQ=255)                   tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg                                            <  1:3643951/76‑1 (MQ=255)                    ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCCTgg                                           >  2:4989531/1‑76 (MQ=255)                    ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTCTGTTTTTCATgg                                           >  1:4791917/1‑76 (MQ=255)                     tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt                                          <  1:5695099/76‑1 (MQ=255)                        taattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTAc                                        <  1:4587763/73‑1 (MQ=255)                         aattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                                       <  1:3985182/74‑1 (MQ=255)                         aattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                                       <  1:676782/74‑1 (MQ=255)                          attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                      >  1:3611156/2‑76 (MQ=255)                          attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                      <  2:3808763/75‑1 (MQ=255)                           ttCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAgg                                     >  2:1374736/1‑76 (MQ=255)                           ttCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAgg                                     <  1:520249/76‑1 (MQ=255)                            tCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGt                                    <  2:4791916/76‑1 (MQ=255)                              aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga                                  >  1:4536946/1‑76 (MQ=255)                              aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga                                  <  2:683640/76‑1 (MQ=255)                               gggAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag                                 >  2:4237131/1‑76 (MQ=255)                               gggAATTTTTATGGCTCAAGTTATTAATGAAATGGA‑CGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag                                 >  1:5400432/1‑76 (MQ=255)                                 gAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagag                               <  2:5400431/76‑1 (MQ=255)                                  aaTTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagaga                              <  2:1819838/76‑1 (MQ=255)                                  aaTTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagaga                              <  1:4726971/76‑1 (MQ=255)                                  aaTTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagaga                              >  1:3786425/1‑76 (MQ=255)                                    tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGatt                            >  1:4959610/1‑75 (MQ=255)                                    tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGata                            >  1:4846230/1‑76 (MQ=255)                                    tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGata                            <  1:4730834/76‑1 (MQ=255)                                    gtttgATGGCTCAAGTTATGAATGAAATGGAT‑GTTCCGTCCCATTCGGTGGTGTTTCATGGTACAGGTGAGAGata                            <  1:4756881/75‑1 (MQ=255)                                     ttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGatat                           <  2:483744/76‑1 (MQ=255)                                      tttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtt                          <  1:1205406/76‑1 (MQ=255)                                      tttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtt                          >  1:747557/1‑76 (MQ=255)                                       gtgtGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAttt                         <  2:309864/73‑1 (MQ=255)                                        tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtttt                        <  1:3211826/76‑1 (MQ=255)                                        tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtttt                        >  2:2667652/1‑76 (MQ=255)                                         aTGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTc                       >  1:2865951/1‑76 (MQ=255)                                         aTGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTc                       >  2:88628/1‑76 (MQ=255)                                          tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCt                      >  1:2290311/1‑76 (MQ=255)                                            gCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCttt                    >  2:914399/1‑75 (MQ=255)                                             cTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAt                   <  1:4365007/76‑1 (MQ=255)                                             cTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAt                   <  1:4407522/76‑1 (MQ=255)                                              tCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAtt                  <  2:5528353/76‑1 (MQ=255)                                               cAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAttt                 >  1:5098401/1‑76 (MQ=255)                                               cAAGTTATTAATGAAATGGAT‑GTGCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAttt                 <  2:2655735/76‑1 (MQ=255)                                                 aGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgt               <  1:5707117/76‑1 (MQ=255)                                                  gTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgtg              >  2:3262337/1‑76 (MQ=255)                                                  gTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgtg              <  2:3447524/76‑1 (MQ=255)                                                     attaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtgg           >  2:2261594/1‑76 (MQ=255)                                                      ttaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtggt          >  1:2295420/1‑76 (MQ=255)                                                       taATGAAATGGAT‑GTGCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTgtggtg         <  2:815752/76‑1 (MQ=255)                                                        aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGa        >  2:3617816/1‑76 (MQ=255)                                                        aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGa        >  1:3381248/1‑76 (MQ=255)                                                         aTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa       >  1:668844/1‑76 (MQ=255)                                                         aTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa       >  2:3490714/1‑76 (MQ=255)                                                          tGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAt      >  1:2817275/1‑76 (MQ=255)                                                           gAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg     >  1:1608326/1‑76 (MQ=255)                                                            aaaTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtgt    <  1:2687390/76‑1 (MQ=255)                                                             aaTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtgtg   >  2:5809944/1‑76 (MQ=255)                                                              aTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATgtgt  >  2:29030/1‑76 (MQ=255)                                                                    |                                                                       TAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGTGT  >  NC_000913/4477254‑4477391

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