Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,897,126	T→C	T249A (ACG→GCG)	ygcS ←	putative transporter YgcS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,897,126	0	T	C	97.4%	106.7 / ‑6.4	39	T249A (ACG→GCG)	ygcS	putative transporter YgcS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/18);  minor base G (0/1);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCG  >  NC_000913/2896882‑2897367                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  aaaCAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCAGCGcgccg                                                                                                                                                                                                                                             >  2:86551/1‑251 (MQ=255) aaaCAGGCCATGACTACCAGCGTTGACGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAACTCACAAACCATGGGATTACGAGGCAGACAACGAAGCCGCGGTTAAACGCCGCGcgccg                                                                                                                                                                                                                                             >  2:194205/1‑251 (MQ=255)     aGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCgccgcc                                                                                                                                                                                                                                           >  1:89623/1‑249 (MQ=255)           tGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAAcgc                                                                                                                                                                                                                                   >  1:115376/1‑251 (MQ=255)                  cAGCGATGCCGCCAGCAGCCAAAAACATCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCATAACTAATCCCACCAGCGCGCCCACAATTAACAACGCATTAAGCAGCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAACCACAAACCAAGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAgaga                                                                                                                                                                                                                            <  1:125322/251‑1 (MQ=255)                                          aCTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGCTTTGATGTGTTTATggg                                                                                                                                                                                                    >  1:106700/1‑251 (MQ=255)                                               ccAGCAAAAATTGGCCATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCGGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGcc                                                                                                                                                                                               <  2:150868/251‑1 (MQ=255)                                                           tGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTc                                                                                                                                                                                    >  1:197572/1‑250 (MQ=255)                                                                     cagcagGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTGTATGGGTCGCCGTTACCACTTCATCGCCCAGTa                                                                                                                                                                         <  2:95832/251‑1 (MQ=255)                                                                              cGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAACATTAGGCTGACAGTCAGCGCACATACCACCCCAACAGTCTGGGCGATAGTTGACAGCCAGATCTAACTCCCACACCATGGGATTAAAAGGCAGACAACGAAGACGCTGCTAACCGCCGCGCGCCACCAGTAACGCGAGGAGAACAAGGATTTAATGAGCTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATggg                                                                                                                                                                >  1:113109/1‑251 (MQ=255)                                                                                    aaCTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTGACCACTTCATCGACCAGTAAACAATGGGGACCa                                                                                                                                                           >  2:72924/1‑250 (MQ=255)                                                                                                cagcGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCa                                                                                                                                              <  1:251420/251‑1 (MQ=255)                                                                                                 agcGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCAc                                                                                                                                             <  1:211266/251‑1 (MQ=255)                                                                                                     cgcCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGAt                                                                                                                                          >  1:222026/1‑250 (MQ=255)                                                                                                              ttAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGCCCAAAATAGCGATGCACGATAGCGTGAGCt                                                                                                                                >  2:84638/1‑251 (MQ=255)                                                                                                                        aTTAAGCATCAGGCTGGCAGTAAGCGCATCTTCCAGTCCAAAAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCaaacc                                                                                                                      <  1:151718/251‑3 (MQ=255)                                                                                                                                       ggCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTGTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTa                                                                                                       <  1:266452/251‑1 (MQ=255)                                                                                                                                        gCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTTATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCTTGAGCTTCTGCAAAACTCCCCTGGGGCAGTAg                                                                                                      <  2:265383/251‑1 (MQ=255)                                                                                                                                         cAGTCAGCGCATCTTCCAGCCCAATAGTCTGGGAGATAGTTCGCAGCCAGGAATAAATCACAAACCATGGGATTACGAGGCAGACAACGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGCGACCTTCTGCAAAACGCCCCTGGCGCAGTAGc                                                                                                     >  1:212210/1‑251 (MQ=255)                                                                                                                                             aaGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGc                                                                                                 <  2:202517/250‑1 (MQ=255)                                                                                                                                                    tCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCgtggtg                                                                                           >  1:134687/1‑250 (MQ=255)                                                                                                                                                       tCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAACCCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGa                                                                                          <  2:49991/248‑1 (MQ=255)                                                                                                                                                            aCCAATCGTCTGGGCGCTAGATGGCAGCCAGGCATAAATAACAACCCATGGGATTACGAGGCAGACAACGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTACCGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGAGAGCTTCTGCAAAACGACCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCg                                                                                  >  2:232813/1‑251 (MQ=255)                                                                                                                                                                               gTTGGCAGCCAGGTATAAATCACAACCCATGGGATTACGAGGCAGACAAAGAAGTCGCTGTTAAACGCCCGGCGCCGCCAGTAACGCGAAGAGAACAAGGTGTTGATGTGTTTAGGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGAGAGCGTGAGCTGCTGCAAAACGCCCCTGGCGCAGGAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCCTTATAACGt                                                               <  2:135779/251‑1 (MQ=255)                                                                                                                                                                                 tGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTa                                                              >  2:28021/1‑250 (MQ=255)                                                                                                                                                                                  ggCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAAt                                                            >  2:105509/1‑251 (MQ=255)                                                                                                                                                                                               aaaTCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTg                                                                                                                                    >  1:30623/1‑166 (MQ=255)                                                                                                                                                                                               aaaTCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTg                                                                                                                                    <  2:30623/166‑1 (MQ=255)                                                                                                                                                                                                     caAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCt                                         <  1:29963/251‑1 (MQ=255)                                                                                                                                                                                                     aaaaCCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCt                                         <  2:197572/250‑1 (MQ=255)                                                                                                                                                                                                        aCCATGGCATGACGAGGCAGATAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAATAGAACAAGGTGTTGAGGTGTTGATGGGTCGCCGTTCCCACTTCATCGCCCAGTAAAACATGGGGACCAAAATTGCGATGCACGATAGCGTGAGCTTCTGCAAAACTCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAAGAGCGTAATCaataactctcgaagagctgcc                                      <  2:221268/251‑16 (MQ=255)                                                                                                                                                                                                         ccATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCg                                     >  1:1924/1‑251 (MQ=255)                                                                                                                                                                                                             gggATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAAAATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAAAAACGCGGGCAGAGCTGCCGATGc                                 >  1:162171/1‑251 (MQ=255)                                                                                                                                                                                                               gATTACTATGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACTCGAAGAGAACAAGTTTTTGATGTTGTTATGGGTCGCCGTTACCTCTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCTCGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAACCATCAACGCGGGCAGAGCTGCCGATGCCa                               <  2:148393/251‑1 (MQ=255)                                                                                                                                                                                                                     gAGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGcc                         >  2:83534/1‑251 (MQ=255)                                                                                                                                                                                                                      gggCAGACAAAGAAGACGCTGTTAACAGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGAGGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCTATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCa                        <  2:106700/250‑1 (MQ=255)                                                                                                                                                                                                                      aGGCAGACAAAGAAGACGCTGTTAAACGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGCTGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGcc                         >  1:152779/1‑250 (MQ=255)                                                                                                                                                                                                                                          gTTAAACGCCGCGCGCCGCCAGGAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCGGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGttt     <  1:156218/250‑1 (MQ=255)                                                                                                                                                                                                                                             aaaCGCCGCGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCg  <  2:222026/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  AAACAGGCCATGACTACCAGCGTTGCCGCCAGCAGCAAAAAACTTCCCAGCAAAAATTTGCGATGTGCCAGCAGGTGCGTCAGAACTAATCCCAGCAGCGCGCCCACAATTAACAACGCATTAAGCATCAGGCTGGCAGTCAGCGCATCTTCCAGACCAATAGTCTGGGCGATAGTTGGCAGCCAGGTATAAATCACAAACCATGGGATTACGAGGCAGACAAAGAAGACGCTGTTAAACGCCGTGCGCCGCCAGTAACGCGAAGAGAACAAGGTTTTGATGTGTTTATGGGTCGCCGTTACCACTTCATCGCCCAGTAAAACATGGGGACCAAAATAGCGATGCACGATAGCGTGAGCTTCTGCAAAACGCCCCTGGCGCAGTAGCCAGCGTGGTGATTCTGGCGTTCCCCAGCGTAATAACGTAATCAACAACGCGGGCAGAGCTGCCGATGCCAGTAGCCAGCGCCAGGCCTCCGGGTTTTCG  >  NC_000913/2896882‑2897367

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