Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,743	T→C	G53G (GGT→GGC)	tsaE →	N(6)‑L‑threonylcarbamoyladenine synthase, TsaE subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,743	0	T	C	97.3%	128.7 / NA	37	G53G (GGT→GGC)	tsaE	N(6)‑L‑threonylcarbamoyladenine synthase, TsaE subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/25);  minor base A (1/0);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCGTGA  >  NC_000913/4395496‑4395985                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   ggcgCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGACTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGccat                                                                                                                                                                                                                                                 >  1:318010/1‑251 (MQ=255)    gcgTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAg                                                                                                                                                                                                                                              <  1:354871/251‑1 (MQ=255)           ggAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATGCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGTGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCGATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATg                                                                                                                                                                                                                                       <  1:265971/251‑1 (MQ=255)              aCGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCTCAACCGTAATCTATCTGTATTGCGATTTAGGCTCAGGTAAAACCACCTTTAGCCGGGGCGTTTTACAGGCTCTGGGCCATCAGGGTAATGTCa                                                                                                                                                                                                                                    <  2:408165/251‑1 (MQ=255)                               aaCCGATCTCTTTTCCACGCTAGAGCGTATTTTTAACCCGGAAGTGACTGATAAAAACCATAATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTAGCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACg                                                                                                                                                                                                                   <  1:353754/251‑1 (MQ=255)                                        ctTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTc                                                                                                                                                                                                             <  1:65251/248‑1 (MQ=255)                                                cGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATAc                                                                                                                                                                                                  >  2:266381/1‑251 (MQ=255)                                                                tAACCCGGAAGGGACAGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCGGATGAGCAGGCAA‑CATTAGACCTGGGCGATCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTTTGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAATCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTaac                                                                                                                                                                                   <  2:312173/250‑2 (MQ=255)                                                                 aaCCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATgg                                                                                                                                                                                 <  2:314400/251‑1 (MQ=255)                                                                     cGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGTCTGCGATGGCGCAACCGTAAGCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTGTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTa                                                                                                                                                                             <  1:283534/251‑1 (MQ=255)                                                                         aGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCAc                                                                                                                                                                         >  1:93720/1‑251 (MQ=255)                                                                             aCTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTg                                                                                                                                                                     <  1:91326/251‑1 (MQ=255)                                                                                    aaaaCCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGTCAACTTAATGGTCGATCACTTTGATTTGTa                                                                                                                                                              <  1:105583/251‑1 (MQ=255)                                                                                      aaccaataagaATCGAGGAAATCCGCTCCCTGATGGGGAGGCAAACAT‑AGGCCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCTACAACGTAATGGTCTATCACTTTGATTGGTAcc                                                                                                                                                            <  2:375478/242‑1 (MQ=255)                                                                                                 aTCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGAt                                                                                                                                                 <  2:50920/251‑1 (MQ=255)                                                                                                   cGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATcc                                                                                                                                               >  2:106488/1‑251 (MQ=255)                                                                                                    gAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATccc                                                                                                                                              >  2:113987/1‑251 (MQ=255)                                                                                                      gTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCg                                                                                                                                             >  2:170989/1‑250 (MQ=255)                                                                                                                 cccTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGt                                                                                                                                 <  1:468/251‑1 (MQ=255)                                                                                                                             ggCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTTGAGTTTATGGGGATcc                                                                                                                     <  1:303231/251‑1 (MQ=255)                                                                                                                                         aCCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTAtttt                                                                                                          <  1:266381/251‑1 (MQ=255)                                                                                                                                              ggagtgCGGGTAGCTAAAGCCTGCGATGGCGCAACCGGTATCTATGTGTATGGCGATTGTGGCGCAGGTAAAGCCACCGTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCATGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTAGCGCTTTGAGTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGGTTATTTTGccgc                                                                                                     <  2:113947/246‑3 (MQ=255)                                                                                                                                                  ggCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGAt                                                                                                 <  1:240738/250‑1 (MQ=255)                                                                                                                                                  aGCGGGTAGCGAAAGCATGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGCGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGAt                                                                                                 >  2:353558/1‑251 (MQ=255)                                                                                                                                                          gCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCTATTTATGTGCATGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGTCCACCAGGGTAATGTCAAAAGCCCCACGTATACGCTGTTCGAACCCTATACGCTCGACAACTTAAGGATCTATCACTTTTAGTTGTACTGCCTTGCCGATCCCGATGAGCt                                                                                                                                      <  2:378047/206‑1 (MQ=255)                                                                                                                                                          gCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACATTTCGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTACTGTCAAAAACCACACTCATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCCCTTTGATTTGTACCGCCTTGCCGATCCCGAGGTGCt                                                                                                                                      >  1:378047/1‑206 (MQ=255)                                                                                                                                                                              aaCCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaac                                                                     >  1:374570/1‑251 (MQ=255)                                                                                                                                                                                        tctctGTATGGAGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGAACAgg                                                           >  2:421844/3‑251 (MQ=255)                                                                                                                                                                                             gTATGGCGATTTAGTCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGACATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGgtg                                                         <  2:116484/248‑1 (MQ=255)                                                                                                                                                                                                   cgagttAGGCGCAGGTAAAACCACCTTTAGCAGGGGCTTTTTACAGGCTCTGGGCCATCTGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTGTGATTTGTACCGCCTTGCCGATCTCGATGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAATGACGCCAGCTGCCTGGTGGAGTGGCCACAACAAGGGACAGGTGTTCTTCCTg                                                <  1:424065/247‑1 (MQ=255)                                                                                                                                                                                                     aTTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGAc                                              <  1:163562/251‑1 (MQ=255)                                                                                                                                                                                                        tAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGAccc                                            <  1:170989/250‑1 (MQ=255)                                                                                                                                                                                                                     aaCCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGCAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGCCAACTTACTGGTCTATCACTTTGCTTTGTACCGCCTTGCCGCTCCCGAGGAGCTGGCGTTTCTGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATa                               >  2:395189/1‑250 (MQ=255)                                                                                                                                                                                                                      accaccTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATaca                             <  1:175685/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggggCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATc                  >  1:4866/1‑249 (MQ=255)                                                                                                                                                                                                                                         ttttACAGGCTCTGGGACATCAGGGGAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGcaca           >  2:129947/1‑250 (MQ=255)                                                                                                                                                                                                                                             aCAGGCTCTGGTCCATCAGGTTAATGTCAAAACCCCTACTTATACGCGTGTCTAACCCTATACGCTCGACAACTTAATGGGCTATCACTTTTATTTGTACCGCCTTTCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGATGGAGTGTCCAAAACAAGTTACAGGTGTTCTGCCTTACCCGGATGTCGAAATACACACTGATTATCAGGctgaaggc       <  2:180246/250‑8 (MQ=255)                                                                                                                                                                                                                                                 gcgctGGGCCAGCAGGGAAATGTCAAAAGCCCCTCTTATACGCTGGTCGAACCCTATACGCGCGACAACTTAATGGTCTATCTCTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCGGGGGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGAGTACCAGGCACAAGGCCGTGa  <  2:211034/248‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GGCGCGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAA‑CATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCGTGA  >  NC_000913/4395496‑4395985

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