Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,053,682	A→G	T14A (ACG→GCG)	shiA →	shikimate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,053,682	0	A	G	98.2%	192.0 / ‑4.3	57	T14A (ACG→GCG)	shiA	shikimate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/20);  total (37/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 = 3.14e-01

TGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCC  >  NC_000913/2053546‑2053818                                                                                                                                         |                                                                                                                                         tGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGgcgc                                                                                                                                       <  1:898684/140‑1 (MQ=255)             ctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTc                                                                                                                           >  1:759190/1‑140 (MQ=255)                    ggTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCcg                                                                                                                    <  1:714815/140‑1 (MQ=255)                      ttCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCcgac                                                                                                                  >  2:249282/1‑140 (MQ=255)                         gTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAg                                                                                                               <  2:897833/140‑1 (MQ=255)                         gTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAg                                                                                                               <  2:818842/140‑1 (MQ=255)                         gTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAg                                                                                                               <  2:198573/140‑1 (MQ=255)                         gTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAg                                                                                                               <  2:114529/140‑1 (MQ=255)                         gTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAg                                                                                                               <  2:1105644/140‑1 (MQ=255)                               aGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGt                                                                                                         >  1:325322/1‑140 (MQ=255)                                gATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGtt                                                                                                        >  2:540473/1‑140 (MQ=255)                                 aTCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTa                                                                                                       <  2:759190/140‑1 (MQ=255)                                   cGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAgg                                                                                                     >  2:1035135/1‑140 (MQ=255)                                     aCGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCa                                                                                                   >  1:581628/1‑140 (MQ=255)                                     aCGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCa                                                                                                   >  1:723429/1‑140 (MQ=255)                                                 gTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTg                                                                                       <  1:1035135/140‑1 (MQ=255)                                                  ttACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGc                                                                                      <  1:624471/140‑1 (MQ=255)                                                   tACCTTTTCCATACTCATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTg                                                                                       <  2:1035760/138‑1 (MQ=255)                                                   tACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTg                                                                                       >  1:1035760/1‑138 (MQ=255)                                                    aCCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCcg                                                                                    <  2:730277/140‑1 (MQ=255)                                                          tCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGtcgtcg                                                                              >  2:694349/1‑140 (MQ=255)                                                               aCTAATTATAAAAAGACAGTATAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTgg                                                                         >  2:169233/1‑140 (MQ=255)                                                               aCTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTgg                                                                         >  1:706532/1‑140 (MQ=255)                                                                   aTTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATg                                                                     >  2:766650/1‑140 (MQ=255)                                                                   aTTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATg                                                                     >  1:284193/1‑140 (MQ=255)                                                                   aTTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATg                                                                     >  1:259298/1‑140 (MQ=255)                                                                         aaaaGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGAttttt                                                               <  2:384331/140‑1 (MQ=255)                                                                           aaGACATTACAAACAGGATCATTTTGGATTCCACGCTCTTCTCCTCTCGTCCCGATGAAGGGGCGGTTTCGTTAAGTCGCGCATGACGAGTTTCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTa                                                             <  2:785932/140‑1 (MQ=255)                                                                                  tACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATg                                                      >  1:964202/1‑140 (MQ=255)                                                                                      aaCAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTCCTCTATGGCAt                                                  <  2:325322/140‑1 (MQ=255)                                                                                      aaCAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCAt                                                  <  2:723429/140‑1 (MQ=255)                                                                                            aTCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTa                                                             >  1:211380/1‑123 (MQ=255)                                                                                            aTCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTa                                                             <  2:211380/123‑1 (MQ=255)                                                                                            aTCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCAccgc                                            >  1:161563/1‑140 (MQ=255)                                                                                              cATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCAccgccg                                          <  1:249282/140‑1 (MQ=255)                                                                                               aTTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACcgccgc                                         <  1:694349/140‑1 (MQ=255)                                                                                                 tATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTTTGATTTTTTCCTCTATGGCATCACCGCCGCAc                                       <  2:706532/140‑1 (MQ=255)                                                                                                 tATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTCCTCTATGGCATCACCGCCGCAc                                       <  2:284193/140‑1 (MQ=255)                                                                                                     gACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCAcctgt                                   >  1:1032858/1‑136 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  1:206956/1‑140 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  1:758054/1‑140 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  1:426554/1‑140 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  1:70646/1‑140 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  1:183228/1‑140 (MQ=255)                                                                                                        tCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGtt                                >  2:925933/1‑140 (MQ=255)                                                                                                          cACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTa                              >  1:25243/1‑140 (MQ=255)                                                                                                          cACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTa                              >  1:17674/1‑140 (MQ=255)                                                                                                          cACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTa                              >  2:749112/1‑140 (MQ=255)                                                                                                          cACGCTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCGGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTa                              <  2:758054/140‑1 (MQ=255)                                                                                                          cACGCTCATCTCCACTCGTCCCGATGAAGCGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTa                              >  1:490975/1‑140 (MQ=255)                                                                                                              cTCATCTCCACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATcg                          <  2:964202/140‑1 (MQ=255)                                                                                                                     ccACTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGttt                   >  1:65834/1‑140 (MQ=255)                                                                                                                       aCTCGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGttttt                 >  2:20236/1‑140 (MQ=255)                                                                                                                          cgcccggATGAGGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGTTTTTTTACTTTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTccc              <  2:581628/134‑1 (MQ=255)                                                                                                                          cGTCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTccc              <  2:469396/140‑1 (MQ=255)                                                                                                                            tCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTc                >  2:404119/1‑136 (MQ=255)                                                                                                                            tCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTc                <  1:404119/136‑1 (MQ=255)                                                                                                                            tCCCGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGc            >  2:341137/1‑140 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTcc               <  2:466916/136‑1 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTcc               >  1:466916/1‑136 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCa           >  2:619948/1‑140 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCa           >  2:560519/1‑140 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCa           >  2:525153/1‑140 (MQ=255)                                                                                                                             cccGATGAAGGGGCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCa           >  1:910012/1‑140 (MQ=255)                                                                                                                                      ggggCGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCGGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGccc  <  1:20236/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCC  >  NC_000913/2053546‑2053818

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