Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,217,252	G→T	5.9%	L48L (CTG→CTT)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,217,252	0	G	T	5.9%	238.3 / 2.5	68	L48L (CTG→CTT)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (51/13);  new base T (2/2);  total (53/15)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGG  >  NC_000913/4217116‑4217391                                                                                                                                         |                                                                                                                                            cccGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGgc                                                                                                                                          >  2:576796/1‑140 (MQ=255)      acacAACAAATTGAAGAATTACAGAAATAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACTCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAga                                                                                                                                     >  1:62923/1‑139 (MQ=255)            caaATTGAAGACTTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAatat                                                                                                                               <  2:1454490/140‑1 (MQ=255)            caaATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAAt                                                                                                                               <  2:1246696/140‑1 (MQ=255)                   aagaaTTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGt                                                                                                                        <  1:1547761/140‑1 (MQ=255)                                  agTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAg                                                                                                         >  1:822375/1‑140 (MQ=255)                                            aCCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACGGGGCGCAGCGAACATGAGGCGTCTCCTGCACGGGGAGTCGAAAAAAg                                                                                               >  2:740062/1‑140 (MQ=255)                                                    gggAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATc                                                                                       <  1:47741/140‑1 (MQ=255)                                                    gggAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATc                                                                                       <  1:46223/140‑1 (MQ=255)                                                    gggAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATc                                                                                       <  1:48044/140‑1 (MQ=255)                                                           cATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGcc                                                                                >  1:659729/1‑140 (MQ=255)                                                           cATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGcc                                                                                >  1:659612/1‑140 (MQ=255)                                                                    cccATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCac                                                                       >  1:1496095/1‑140 (MQ=255)                                                                    cccATACAGTGCGGAAGATGTGGTGAAATTACGCGGTACAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCac                                                                       >  2:767087/1‑140 (MQ=255)                                                                     ccATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCact                                                                      >  2:1102819/1‑140 (MQ=255)                                                                        tACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCactgac                                                                   >  1:1380541/1‑140 (MQ=255)                                                                        tACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCactgac                                                                   >  1:884200/1‑140 (MQ=255)                                                                         aCAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTTAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCactgact                                                                  <  1:206396/140‑1 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAATATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:196515/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAATTGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:387023/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGTAAAATGGCGACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  2:802582/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:669308/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:67343/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:70190/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:212865/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  2:1347484/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:328212/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  2:1566601/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:148928/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  2:646552/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:1355939/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:1355491/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:1132362/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:1041355/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  1:1070133/1‑140 (MQ=255)                                                                             tGCGGAAGATGTGGAGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCg                                                              >  2:849831/1‑140 (MQ=255)                                                                               cGGAAGATGTGGTTAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTTGGCGCAg                                                                                                                                      <  1:1469413/66‑1 (MQ=255)                                                                               cGGAAGATGTGGTTAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTTGGCGCAg                                                                                                                                      <  1:1469417/66‑1 (MQ=255)                                                                               cGGAAGATGTGGTTAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTTGGCGCAg                                                                                                                                      >  2:1469417/1‑66 (MQ=255)                                                                               cGGAAGATGTGGTTAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTTGGCGCAg                                                                                                                                      >  2:1469413/1‑66 (MQ=255)                                                                               cGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGt                                                            >  1:747575/1‑140 (MQ=255)                                                                               cGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGt                                                            >  2:673005/1‑140 (MQ=255)                                                                               cGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGt                                                            >  2:476371/1‑140 (MQ=255)                                                                               cGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGt                                                            >  2:1622593/1‑140 (MQ=255)                                                                               cGGAAGATGTGGTGAAATTACGCGGTACAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGAAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGTCTACATCAACAGCCTCGGCGCACTGACTGGCGGt                                                            >  2:1079915/1‑140 (MQ=255)                                                                                  aaGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCTGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAg                                                         >  1:1424517/1‑140 (MQ=255)                                                                                  aaGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAg                                                         >  1:1477977/1‑140 (MQ=255)                                                                                     aagtggtgAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGTGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACAACAACAGCCTCGGCGCACTGACTGGCGGTCAGgcg                                                      <  2:1481760/138‑1 (MQ=255)                                                                                           tgAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCaa                                                <  1:405411/140‑1 (MQ=255)                                                                                              aaTTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAg                                             <  1:821317/140‑1 (MQ=255)                                                                                                ttACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGc                                           >  1:199629/1‑140 (MQ=255)                                                                                                        ttCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCgg                                   >  1:672857/1‑140 (MQ=255)                                                                                                                 tCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAg                          <  2:148928/140‑1 (MQ=255)                                                                                                                  ccTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGc                         >  1:460041/1‑140 (MQ=255)                                                                                                                  ccTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGc                         >  1:664745/1‑140 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCATCGAAAATGTGGCGTCTGCTTCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCACGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGCGTATTGAAGCAGTc                     >  2:1068969/1‑140 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGcc                                                                                >  2:1504159/1‑81 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGcc                                                                                <  1:1504159/81‑1 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTc                     >  2:1384773/1‑140 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTc                     >  2:1068874/1‑140 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTc                     >  2:1068892/1‑140 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTAAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAgac                     >  2:18455/1‑138 (MQ=255)                                                                                                                      aaTGCACGCTGGCGCAACTGCGCGCAGCGAAATTGTGGCGTCTGCTGCACGGGGAGGCGAAAAAAGGCTACATCCAAAGCCTCGGCGCCCTGCCTGGCGGTCAGGCGCTGCAACTGGCGAAAGCGGGGATTGAAGCAggc                     >  2:18608/1‑138 (MQ=255)                                                                                                                       aTGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCt                    >  2:1240293/1‑140 (MQ=255)                                                                                                                           aCGCTGGCGCATCTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGCAAGCGGGGATTGAAGCAGTCTATCt                >  2:1453464/1‑140 (MQ=255)                                                                                                                           aCGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTAACt                >  2:1263544/1‑140 (MQ=255)                                                                                                                             gCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGt              <  2:361324/140‑1 (MQ=255)                                                                                                                                  cgcAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGa         >  2:509389/1‑140 (MQ=255)                                                                                                                                         gggCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAgg  <  1:1102819/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                            CCCGTACACAACAAATTGAAGAATTACAGAAAGAGTGGACTCAACCGCGTTGGGAAGGCATTACTCGCCCATACAGTGCGGAAGATGTGGTGAAATTACGCGGTTCAGTCAATCCTGAATGCACGCTGGCGCAACTGGGCGCAGCGAAAATGTGGCGTCTGCTGCACGGTGAGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGG  >  NC_000913/4217116‑4217391

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