Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	460,060	C→T	I391I (ATC→ATT)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	460,060	0	C	T	95.9%	124.1 / NA	49	I391I (ATC→ATT)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/20);  minor base G (1/1);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATC‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAA  >  NC_000913/459813‑460305                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       tCAAAAAAGCCCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCCGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGt                                                                                                                                                                                                                                                    >  2:176738/1‑251 (MQ=255)  cAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCTTGCGTGATGACGCGGAAATT‑CGTg                                                                                                                                                                                                                                                   >  2:292544/1‑251 (MQ=255)      aaaGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGCCTGTAGCGCGTGAAAGAGCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTTTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCa                                                                                                                                                                                                                                               <  2:92452/251‑1 (MQ=255)                gTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACt                                                                                                                                                                                                                                     >  1:166258/1‑251 (MQ=255)                            aaaTCATTGATACCGAACAATATGGTATGGATCGCGTGAAAGATCGAATCCTTGAGTATATTGGGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAACCCTCTCTTGGGCAGTCCATTGCAAAAGACACCGGGCGTAAATAAGTCCGGATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGATAACCGCCGCACTTACATCGGTTCt                                                                                                                                                                                                                         >  2:394996/1‑251 (MQ=255)                                ccTTGATACCGCCCATTATGGTTCGGCGAGCGTGAAAGATCGCATACTTGCGTATCCTGCAGTTCAAAGCCGTGTAAACAAAACCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACATCTCTTGGTCAGTCCATTGCCAAAGACACCGGGCGTAAATATGACCGTATGGCGCTGGTCGGCGTGCGTGATGAAGCGGAAATT‑AGTGGTAACCGACGTCCTTACATCGGGTCTATGc                                                                                                                                                                                                                     >  2:171863/1‑251 (MQ=255)                                 cTTGATACAGACCATTATGGTCTGGAGCGCGTGCAAGATCGAATCCTTGAGTATCTTGCGGCTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑AGTGGTCACCGCCGTACTTACATCGGTTCTATGcc                                                                                                                                                                                                                    >  2:157605/1‑251 (MQ=255)                                    gATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCgg                                                                                                                                                                                                                  <  1:96393/250‑1 (MQ=255)                                             cATTATGGTCGGGAGCGCGTGAAAGATCGAAGCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGGGTAAATATTTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGAt                                                                                                                                                                                                        <  1:341484/251‑1 (MQ=255)                                                        ggAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATgg                                                                                                                                                                                             >  2:353200/1‑251 (MQ=255)                                                            cgcgAGAAAGATCGAATCCTTGATTATCTTGCGGTTCAAAGCCGTGCCAACAAAACCAAGGGACCGATCCTCGGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCTTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGaa                                                                                                                                                                                         <  1:278413/251‑1 (MQ=255)                                                               gTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGt                                                                                                                                                                                      >  1:361217/1‑251 (MQ=255)                                                                tGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTg                                                                                                                                                                                     >  1:244204/1‑251 (MQ=255)                                                                tGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTg                                                                                                                                                                                     >  2:168501/1‑251 (MQ=255)                                                                   aaGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGCCCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAACCCTCTCTTGGTCCGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTAGGGCGCTGGGCGGCGTGCGTGATGAGGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGAGCCAGAAAATGGCGAAAGTGGGc                                                                                                                                                                                  <  2:267361/251‑1 (MQ=255)                                                                        cGAATCCTTGAGTATCTTGCGGTTCAACGACGTGGCAACAAAAACAAGGGACCGATCCACTGCATGGGCCAACCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGCCCATTGCCCAAGCCACCGGGCGTAAATATGTCCGTAAGGAGCTGGACGGCGTGCGTGATGAAGCGGAAATG‑CGTGATAACCGGCGTACTTACATCCGTTATATGCCGGGTAAACTGATCCAGAAAAAGGCGAAAGGGGGCGTGa                                                                                                                                                                              >  2:288679/1‑250 (MQ=255)                                                                          aaTCCTCGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGGGGGCGTGaaa                                                                                                                                                                            >  1:118438/1‑250 (MQ=255)                                                                                      tCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGGGGGCGTGAAAAACCCGCTGTTcc                                                                                                                                                               >  1:230711/1‑251 (MQ=255)                                                                                                                aaaTCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGtt                                                                                                                                                                                                                           >  1:75234/1‑165 (MQ=255)                                                                                                                aaaTCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGAAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGtt                                                                                                                                                                                                                           <  2:75234/165‑1 (MQ=255)                                                                                                                    aaagggccgtaTCCTCTGCCTGGTAGGGCCGCAGGGGGTAGGTAAATCCTCTCTTGGTCTGTCCTTTGCCAAAGCCACGGGGGTTAAATATGTCCGTATGTCGCTGGGCGGCGTGCGTGAGGAAGCGGAAATG‑CGGGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTg                                                                                                                                 <  1:270491/241‑1 (MQ=255)                                                                                                                         gaccgtcccTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAACCCTCTTTTGGTCAGTCCGTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATg                                                                                                                            <  2:166258/244‑1 (MQ=255)                                                                                                                         gACCGATCCTCGGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTCAACTGATCCAGAAAATGGCGAAAGTGGGCGTGCAAAACCCGCTGTTCCTGCTCGATGCGATCGACAAAATGTCTTCTGACATg                                                                                                                            >  2:321223/1‑251 (MQ=255)                                                                                                                            cgacccgctGCGTGGGAGGGCCGCCGGGGTTAGGAAAAACCTCTTTTGGTCAGGCCATTGCAAAAGCCACCGGGGGTAAATAGGGCGGTATGGCGCGGGGCGGGGGGCGGGAGGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGt                                                                                                                         <  2:224135/244‑1 (MQ=255)                                                                                                                                ccTCTGCCGGTTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCg                                                                                                                     <  1:453025/251‑1 (MQ=255)                                                                                                                                   ctGCCTGGTAGGGCAGCCGGGGGTAGGCAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGACGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAACGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAACATGTCTTATGACATGCGTGGCGATc                                                                                                                  >  2:297135/1‑251 (MQ=255)                                                                                                                                     gCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCGCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATc                                                                                                                  <  1:6320/249‑1 (MQ=255)                                                                                                                                         ggTAGGGCCGCCGGGGGTAGGTATACCCTCTCTGTGCCATTCCATTCCCAAACCCACGGGGTGTAATTATGTCTGTATGGCGTTGGGCGTGGTGCTGGATAAAGCGGAAATT‑GGTGTTAACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCTAAAGTGGGCGTGAAAAGCCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGcc                                                                                                             <  2:436215/250‑1 (MQ=255)                                                                                                                                                     gggggTAGGTAAAACCTCTCTTGGTCTTCCCATTGCCAAAGCCACCGGGGGTAAATATGTCCGTATGGCGCTGGGGGGCGTGCGGGATGAGGCGGAATTT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTg                                                                                                <  2:191416/251‑1 (MQ=255)                                                                                                                                                         gtaggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGt                                                                                            >  1:283995/1‑251 (MQ=255)                                                                                                                                                                aacccTCTCTGGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGGCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCTCCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGAt                                                                                     <  1:176738/248‑1 (MQ=255)                                                                                                                                                                           ggccggTCCATTGCCAAAGCCCCCGGTCTTAAATATTTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGCCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGaa                                                                          <  1:58061/246‑1 (MQ=255)                                                                                                                                                                           ggTCAGTCCATTACCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATT‑CGAGGTGACCGCCGTACGTACATCGGTTCTATGCCGGGTAAACTAATCCAGAAAATGGCGAAAGTGGGCGTGAGAAACGCGCAGTTCCTGCTCGATGAGATCGACAAAAAGTCTTCTGACATGCGTGGCGATCCGGACTCTGCACTGCTTGAAGTGCTGGATCCCGAGCAGa                                                                           >  2:326912/1‑250 (MQ=255)                                                                                                                                                                             tCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACg                                                                        >  1:376668/1‑251 (MQ=255)                                                                                                                                                                              cAGTCCTTTGCCAAAGCCTCCGGGCGGAAAGATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGTTCACCGCCTTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCTATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGACCCGGCCTCTGCACTGCTTGAAGTGCTTGATCCAGAGCAGAACg                                                                        <  1:292296/250‑1 (MQ=255)                                                                                                                                                                               aGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGt                                                                       >  1:322599/1‑250 (MQ=255)                                                                                                                                                                                                   gggCGTAAACATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCt                                                  >  1:320365/1‑251 (MQ=255)                                                                                                                                                                                                      cGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTgg                                                >  1:160317/1‑250 (MQ=255)                                                                                                                                                                                                          aaTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTg                                           <  1:52173/251‑1 (MQ=255)                                                                                                                                                                                                          aaTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTg                                           >  1:238811/1‑251 (MQ=255)                                                                                                                                                                                                              gtcccGTATGGCGCTGGGCGGCGGGCGTGATGATG‑GGAAATTTCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGAt                                        <  2:154889/247‑1 (MQ=255)                                                                                                                                                                                                                 ccGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGGGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGGGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGCAGGGGATTACg                                    >  1:52889/1‑251 (MQ=255)                                                                                                                                                                                                                   gTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGAt                                  >  1:146873/1‑251 (MQ=255)                                                                                                                                                                                                                       ggCGCTGGGCGGCGTGCGTGATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGGGGATTACGATCTCa                              <  1:300453/251‑1 (MQ=255)                                                                                                                                                                                                                                  tcgagcgtgATGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCGATGCCGGGTAAACTGATCCAGAAAATAGCGAAAGAGGGCGTGCAAAACCCGCTGTTCCTGCTAGATGAGATCGACAAACTGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCCCTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGCCCCCTCCCTGGAAGAGGATTACGATCTCAGCgtctagatg                   >  2:338551/5‑243 (MQ=255)                                                                                                                                                                                                                                   cgtgcgtgATGAAGCGGAAATT‑CGTGGTCACCGACGTACTTACATCGGTTCTATGCCGGGTAAACTGATACAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGGCGTGATGt                  >  2:196466/1‑251 (MQ=255)                                                                                                                                                                                                                                      gcgtgAGGAAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGttt                <  1:292544/250‑1 (MQ=255)                                                                                                                                                                                                                                          gatgaAGCGGAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcgc           >  2:133063/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gAAATT‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGaa  <  1:353200/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       TCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATC‑CGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAA  >  NC_000913/459813‑460305

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