Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,521,151	C→T	G290G (GGG→GGA)	yjhF ←	putative transporter YjhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,521,151	0	C	T	90.7%	115.3 / ‑5.4	43	G290G (GGG→GGA)	yjhF	putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (17/22);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGTAGAGATCGCCGCGACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATAT  >  NC_000913/4520905‑4521392                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ggcggTAGAGATCGCCGCGACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGGGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGACAGGGCTTGTTTTCCCGCGtccc                                                                                                                                                                                                                                                 >  2:443463/1‑250 (MQ=255)  gcggtaaggaTCGCCGCGAAGGTGGCGGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCAATTCCCGAAACCAGGGTGGAAATATACGGACCAACTCCAGAGTCGATGAGAACTTGCTAGAAAGCACCTCCGCCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGACTCCATCAGGTCATGGATTGTTTACCCGCGtcccgg                                                                                                                                                                                                                                               <  2:9917/243‑3 (MQ=255)      tAGAGATCGCCGCGACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGc‑cc                                                                                                                                                                                                                                          >  1:422605/1‑251 (MQ=255)        gagaTCGCCGCTACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCTCCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCAt                                                                                                                                                                                                                                        <  1:27646/251‑1 (MQ=255)           aTCGCCGCGACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGCTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCCGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGt                                                                                                                                                                                                                                     >  2:282932/1‑251 (MQ=255)               ccgcGACGGGGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTACCAGGC‑CCATGGtat                                                                                                                                                                                                                                   >  2:106833/1‑249 (MQ=255)                          tcTGAACGCAGTCAAATACGACTGAATGCAGCCACACCCCGGGCAATTAGGATCAGGATAATATCCATCCGAGAAACCAGGGTGGAATTATACTGACCAACGCCAGAGTCAAAGAGAACTTGTTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATAACCGCGATAGCTATTCCGCACGGATCCACCAGGTCAGTGTTTGTTTTACCGGATCCCAGGC‑CCAAGGTATACATCGCGAATa                                                                                                                                                                                                                       >  2:312859/2‑250 (MQ=255)                             gAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAg                                                                                                                                                                                                                   <  1:79053/251‑1 (MQ=255)                                 ccAGACAAATACGCAGGAATGCGGCCACGCCCAAGGCCATCAGGATCGGGATAATATACATAACCGAAACCAGGGGGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCAACACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTGTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGACTAGCAGTGa                                                                                                                                                                                                                >  2:433000/1‑250 (MQ=255)                                     aCAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTAAAATCa                                                                                                                                                                                                           <  1:125395/251‑1 (MQ=255)                                            cGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGc                                                                                                                                                                                                    <  1:404296/251‑1 (MQ=255)                                            cGCAGGAATGCGGCCACGCCCCAGGCAATCTGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATAGACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTGTCTGCAAGAATCCAACAGGTCAGGGAGTGTGTGCCCGCGTCTCAGGC‑CCATGGTATAGGTCGAGAGGAGCAGTGAAAGCAGCAttgc                                                                                                                                                                                                    <  2:73640/251‑4 (MQ=255)                                             gCAGTAATGCGGCCCCGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCATCATGGCg                                                                                                                                                                                                   <  2:273798/251‑1 (MQ=255)                                             gCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCg                                                                                                                                                                                                   <  1:319814/251‑1 (MQ=255)                                                ggAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATc                                                                                                                                                                                                <  1:28216/251‑1 (MQ=255)                                                ggAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGAACGCGAATAGCAGTGAAATCAGCATGGCGATc                                                                                                                                                                                                <  2:410018/251‑1 (MQ=255)                                                 gaaggCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCTCCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAACAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCg                                                                                                                                                                                               <  2:203424/247‑1 (MQ=255)                                                            cGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCt                                                                                                                                                                                    >  2:107196/1‑251 (MQ=255)                                                                  aGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAAc                                                                                                                                                                              <  1:52160/251‑1 (MQ=255)                                                                            gATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCCCCACCACCACCGATAATCAGCAGCAGACCGGTAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGt                                                                                                                                                                    <  1:308293/251‑1 (MQ=255)                                                                                     aaTATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATa                                                                                                                                                           <  1:446537/251‑1 (MQ=255)                                                                                            aTACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGa                                                                                                                                                    <  2:329391/251‑1 (MQ=255)                                                                                            aTACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGa                                                                                                                                                    <  1:350717/251‑1 (MQ=255)                                                                                                       cAGGGTGGAAATCTACTGGCCCCCGCAAGCGTCTATGAGACCTTGCTTGAATGCACCACCACCACCGATAGTCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGAAATCATCCATAAGGCAAGAGCTTGTTATCACGCGTACCAGGC‑CCAGGGTATAGATCGCGAATAGCAGTGACATCAGCATGGCGATAGTTGAATTACATAAGAACAGTTACACATTGTAAAATAAGCCTGAATCAGCGTCAc                                                                                                                                         >  2:137063/1‑251 (MQ=255)                                                                                                              gtaaTATACTGACCAACGACAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGAACGGCAATACCCGCGATCGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGATTTCCTGCGTCCCAGGA‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGAGATCGTTGAATTACCTAAGAACAGTATCACGTTGTCAAATAAGACTGAATCCGCGGAACTTTTGGc                                                                                                                                  >  2:440367/3‑251 (MQ=255)                                                                                                              gAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCAACCCCACCGATAATAAGCAGATGACCGGCAATAACCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGC‑CCATGGCATAGATCGCGAATAGCAGTGAAATCAGCATGGCGTTCGTGGAATTACCTACGAACACACTCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTggg                                                                                                                                  >  2:440381/1‑250 (MQ=255)                                                                                                                 atatACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGCTAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATg                                                                                                                               >  1:379022/1‑251 (MQ=255)                                                                                                                            aaCGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAAGAATTGTTTTCCCGCGTCCCAGGC‑CCATGGGATAGATCGCGAATAGAAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGAATAACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCTTGAGTCATATgg                                                                                                                     >  2:338736/1‑250 (MQ=255)                                                                                                                                             gaACTTGCTTGAATGCACCCCCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGcc                                                                                                   <  1:443463/251‑1 (MQ=255)                                                                                                                                                    cttgAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAgg                                                                                             >  2:111848/1‑250 (MQ=255)                                                                                                                                                       gAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATc                                                                                          >  1:208720/1‑250 (MQ=255)                                                                                                                                                                cccacccccGATAATCAGCAGCAGCCCGGCAAGTCCCTCGAGAGCTTGTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAAGAGCAGGGAAATCAGCATGGCGAGCGTTGAATTACCTAAAAACAGTACCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCAGGAGTCATCTGGAGGAGGGTGACAACCGCCATCAGGATCACCGgcagca                                                                                <  2:422605/250‑1 (MQ=255)                                                                                                                                                                            aaTCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTgaga                                                                     >  1:453968/1‑250 (MQ=255)                                                                                                                                                                                      agACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATccc                                                          >  1:131586/1‑251 (MQ=255)                                                                                                                                                                                         ccGGAAATCCCCGCGATGGCTTTTCCGCAAGAATCCAACAGGTCAGGGATTGTTTGCCCTCGCCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGAGAAATCAGCATGGCGATCGTTGAATTCCCTAAGACCAGTATCCCGTTGTAAAGTAAGCCTGAATCCGCGCCACTTTTGGCATGAGTCATCTGGATAAGGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCaaa                                                       <  2:1284/251‑1 (MQ=255)                                                                                                                                                                                           ggCAATACCCGCGATAGCTTTTCCGCAAGAATCAATAAGGACAGAGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATCGCAGTGAACTCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAACACAAGCCTGAATCCGCGGCACTTTTGGCATGCGTCATCTGGCTGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGc                                                     >  1:247321/1‑251 (MQ=255)                                                                                                                                                                                           ggCAATACCCGCGATAGCTTTGCAGAAAGAATCCATCAGGGCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGATTCATCTGGATGATGGGGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGc                                                     <  2:302709/251‑1 (MQ=255)                                                                                                                                                                                                  ccagcgATCGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCACAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATACGCGGCACTTTTGGCATGAGTCATATGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGAGCCCAAAGACTGGCAt                                              >  2:42127/4‑251 (MQ=255)                                                                                                                                                                                                     gtgATAGCTTTTCCGCAAGAATCCATTAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGTTATTGATCGCGAATAGCAGTAAAATCAGCATGTCGATCGTTGAATTACCTAAGAACATTATCACGTTGTAAATTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCAGGCAtt                                             <  1:123136/247‑1 (MQ=255)                                                                                                                                                                                                           gCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCAtttcgctttc                                     <  1:282932/251‑1 (MQ=255)                                                                                                                                                                                                                       gAATCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATcggcg                         >  1:346922/1‑251 (MQ=255)                                                                                                                                                                                                                         aTCCATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGACTCTGCggg                       >  1:16966/1‑251 (MQ=255)                                                                                                                                                                                                                            cATCAGGTCAGGGATTGTTTTCCCGCGTCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACTTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGATTCATCTGTATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGATAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTcg                    <  1:433000/251‑1 (MQ=255)                                                                                                                                                                                                                            cATCAGGTCAGGGATTGTGTTCCGGCGTCCC‑GGCACCTTGGGATAGATCGCGATTACCATTGTAATCAGCATGGAGATCGTGGATGTGCACATGACCAGTAGCACGTAGTAAAATAAGACTGAGGCCGCGGCACTGTTGGCACGAGTCATCTGGAAGATGGTGCCAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCGTGCGCTGAATCGGCGGGTcg                    <  2:282751/251‑1 (MQ=255)                                                                                                                                                                                                                                               ttCCCGCGCCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCATCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATTAGTCTTCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTatat  <  1:366984/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GGCGGTAGAGATCGCCGCGACGGTGGCTGAACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGC‑CCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATAT  >  NC_000913/4520905‑4521392

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