Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,126,657	C→T	A52T (GCC→ACC)	priA ←	primosomal protein N'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,126,657	0	C	T	97.5%	134.9 / NA	40	A52T (GCC→ACC)	priA	primosomal protein N'
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/19);  minor base A (1/0);  total (21/19)
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.95e-01

CAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGCATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCGGATCAGTTTGCATACGCTGGTT  >  NC_000913/4126415‑4126891                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           cAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTa                                                                                                                                                                                                                                     >  1:73735/1‑250 (MQ=255) cAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCCCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCATACGGAGTGAGTAAACTCCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCCTCATTTCGCTGGTCTCGCTa                                                                                                                                                                                                                                     <  2:162923/250‑1 (MQ=255)         ccAGATCAACCGCCTGGCTTTGTTCCGTGGCAAACCAATACCACATCAGCGCGTCAGCAGCAGGCCGCCCCTGGCGTAACAACATCGGCAAGGCATGACCCACCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACAACGGCTCACTATCGAGCAATTAGACTACCGATTTTAGCTTATTGCGCGGAAGTTCGCTGGTATCGCTAACTGANAcc                                                                                                                                                                                                                            >  2:258736/2‑251 (MQ=255)            gATCCACCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACcac                                                                                                                                                                                                                          >  2:171942/1‑250 (MQ=255)                 aCCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATccc                                                                                                                                                                                                                    <  1:171942/251‑1 (MQ=255)                   cGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGa                                                                                                                                                                                                                  <  1:125968/251‑1 (MQ=255)                    gCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGAt                                                                                                                                                                                                                 <  1:255219/251‑1 (MQ=255)                       tGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATgc                                                                                                                                                                                                               >  1:277770/1‑250 (MQ=255)                        ggCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATgc                                                                                                                                                                                                               <  1:67593/249‑1 (MQ=255)                          ccTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTc                                                                                                                                                                                                           >  2:235698/1‑251 (MQ=255)                                tCAGTCGCAAACCAGTACCAAATCGGCGCGCTAGCCGCCGGCCGCCCATGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAACCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCCGACGGCGCGAGCAAACACCGGCTCCCTAACCAGCCCTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGGTCGCTGGTATCGCTAACTGTTACCACAATCCCGATGCGCTCctgctg                                                                                                                                                                                                     >  2:178688/1‑251 (MQ=255)                                                 ccACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATTCCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGgcacgc                                                                                                                                                                                    <  1:85221/251‑1 (MQ=255)                                                       cGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGgcacgcgcacg                                                                                                                                                                               <  1:246740/250‑1 (MQ=255)                                                                cgccgcAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAg                                                                                                                                                                     <  1:144980/251‑1 (MQ=255)                                                                           gccCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTc                                                                                                                                                          <  1:211583/251‑1 (MQ=255)                                                                              ccTGGTGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATg                                                                                                                                                       <  1:213774/251‑1 (MQ=255)                                                                                     tagtaAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTg                                                                                                                                                >  2:226693/1‑251 (MQ=255)                                                                                       gtaAAATCGGCAAGGCATGCACCAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCCGCTTTAACCGTCCTGCCTTCTGgc                                                                                                                                              >  2:178036/1‑251 (MQ=255)                                                                                                  aaGGCATGAAACAGCACATCGCCAAGCGGATGATGATAGCAATCTGCCGCCCATAGCAGCCATCGCCAGACGGAGTGAGTACACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATTCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGaccgtc                                                                                                                                   >  2:223470/1‑246 (MQ=255)                                                                                                                     cGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCgg                                                                                                                 <  2:3485/250‑1 (MQ=255)                                                                                                                           tCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGggc                                                                                                          >  2:101242/1‑251 (MQ=255)                                                                                                                              gatgatgatAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGggca                                                                                                         >  2:103952/1‑249 (MQ=255)                                                                                                                                     atAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGt                                                                                                <  1:179548/251‑1 (MQ=255)                                                                                                                                                 ctcccATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACATCGACTACCGCTTTAAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGATTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCATCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGc                                                                                     >  1:263322/3‑250 (MQ=255)                                                                                                                                                        agcagcAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGcatc                                                                             >  2:282334/1‑251 (MQ=255)                                                                                                                                                           agcaATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCatcatc                                                                          <  1:226693/251‑1 (MQ=255)                                                                                                                                                                 cGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGTCAACGTGGGCAACGGGCATAGCATCATCCTGACt                                                                    <  1:82116/251‑1 (MQ=255)                                                                                                                                                                         ggAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTc                                                            >  2:83951/1‑251 (MQ=255)                                                                                                                                                                              gagtAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTgg                                                       <  1:143971/251‑1 (MQ=255)                                                                                                                                                                                          gCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCACTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGAACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATacac                                           >  1:167665/1‑251 (MQ=255)                                                                                                                                                                                             cACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACAtt                                        <  1:263681/251‑1 (MQ=255)                                                                                                                                                                                             cACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACAtt                                        >  1:225447/1‑251 (MQ=255)                                                                                                                                                                                                              aCTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCAAGCGCACGCGACACCAAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGACATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCgg                       >  2:175893/1‑251 (MQ=255)                                                                                                                                                                                                               cTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATacac                                           <  2:77118/230‑1 (MQ=255)                                                                                                                                                                                                               cTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATacac                                           >  1:77118/1‑230 (MQ=255)                                                                                                                                                                                                               cTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGAGGGTTAGTATACACATTGCCGTAGAAGAGTGCgg                       >  2:105297/1‑250 (MQ=255)                                                                                                                                                                                                                 aCCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCGGAt                     <  1:282250/250‑1 (MQ=255)                                                                                                                                                                                                                         tAGCTCATTGAGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCGGATCAGTTTGc             <  1:101242/250‑1 (MQ=255)                                                                                                                                                                                                                                 tGAGCGGGAGTTCGCTGGAATCGCTAACTGATACCACAATCCCGATGGGCTCCTGCTGTTTACCAAACGGCACGCGCACGCGAAACCCAGCTTTAACAGTCATGACTTTTGGCAGCAGATAGTCAAAGGTACGAGAAAGCGGAACGGGAAAGGCAACGTGGGCAAAGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATAAACATTGCCGGAGAAGAGTGAGGAGCAGTTTGCATACGCTgg    >  2:92299/1‑251 (MQ=255)                                                                                                                                                                                                                                   aGCGGCAGTTCGCTGGTATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCGGATCAGTTTGCATACGCTGGtt  >  1:274836/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CAGGCTGTTCAGATCCACCGCCTGGCCTTGTTCAGTGGCAAACCAGTACCACATCGGCGCGTTCGCCGCAGGCCGCCCCTGGCGTAGTAAAATCGGCAAGGCATGAAACAGCACATCGCCAATCGGATGATGATAGTAATCTGCCGCCCATAGCAGCAATCGCCAGACGGAGTGAGTAAACACCGGCTCACTATCCAGCACTTCGACTACCGCTTTTAGCTCATTGAGCGGCAGTTCGCTGGCATCGCTAACTGATACCACAATCCCGATGCGCTCCTGCTGTTTGCCAAACGGCACGCGCACGCGACACCCAGCTTTAACCGTCATGCCTTCTGGCAGCAGATAGTCAAAGGTACGAGGAAGCGGAACGGGCAAGGCAACGTGGGCAACGGGCATAGCATCATCCTGACTTGAAATTCGGTGGGTTAGTATACACATTGCCGTAGAAGAGTGCGGATCAGTTTGCATACGCTGGTT  >  NC_000913/4126415‑4126891

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