Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,239,343	C→T	T155T (ACC→ACT)	dadX →	alanine racemase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,239,343	0	C	T	96.4%	169.4 / NA	56	T155T (ACC→ACT)	dadX	alanine racemase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/24);  minor base G (0/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGC  >  NC_000913/1239096‑1239586                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     aaGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGTGCAATGGCGAATGTGGGCGAAATGACTCTg                                                                                                                                                                                                                                                   <  1:435132/251‑1 (MQ=255)        gATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGc                                                                                                                                                                                                                                            >  1:361237/1‑251 (MQ=255)             tGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCAtttt                                                                                                                                                                                                                                       >  1:400669/1‑251 (MQ=255)                       aGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCgg                                                                                                                                                                                                                             <  1:503919/251‑1 (MQ=255)                        ggATTTTTAAATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGCCCACCTGCGTACACAGCAACTGACAACTAAAAGCAATGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATTGCGTGATTAATGTCTGGCAGAAGTTGCGGGCAATGGCGAATGTTGGCGAATTGACTCTGATGttggtctttgacgatgcgg                                                                                                                                                                                                                             >  2:113214/1‑232 (MQ=255)                          aTTTTTCCATGCTCAGGATCTGGAGATTTATGACCCGCCCCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAACTGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGTATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCCGCAGTTGCGGGCAATGGCGCATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAAc                                                                                                                                                                                                                          >  2:160335/1‑251 (MQ=255)                                 cATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGa                                                                                                                                                                                                                  <  1:376251/251‑1 (MQ=255)                                    gCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATgg                                                                                                                                                                                                               <  1:273214/251‑1 (MQ=255)                                         ggATCTGGAGATTTATGACCAGCAAGGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGATTCCAGCCCGATCGAGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGTATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAAttt                                                                                                                                                                                                          >  2:422671/1‑251 (MQ=255)                                                           ccAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTa                                                                                                                                                                                        <  2:443865/251‑1 (MQ=255)                                                                     cTGACCACCTGCGTACACAACAACTGGCAGCTCAAAGCCCTGCAAAATGCGCGGCTCAACGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGACTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGGTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGCACATCC‑TGATGGAATTTCAGGAGCGATGGCGCGTATTGCGCAggc                                                                                                                                                                              >  2:354621/1‑251 (MQ=255)                                                                       gACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGCATTGTGCAggcgg                                                                                                                                                                            >  2:184217/1‑251 (MQ=255)                                                                          caccTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGcggcgg                                                                                                                                                                         >  1:464620/1‑251 (MQ=255)                                                                                    cacaGCGACTGGCAGCTCACAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGCGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATAGCGAATGTTGGCGAAATGACTCTGATGGCGCATTTTGCCGAAGCAGAAAAACC‑TGATGGAAATTCAGTCGAGAGGACGCGTATTGGGAAGGCGGCGGCGGGGATggg                                                                                                                                                               >  1:504801/1‑250 (MQ=255)                                                                                            cTGGCAGCTCAAAGCACTGCCAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGc                                                                                                                                                       >  1:464645/1‑251 (MQ=255)                                                                                                  gCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcgttcgt                                                                                                                                                 >  1:299047/1‑251 (MQ=255)                                                                                                              gCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTcgg                                                                                                                                     <  1:466489/251‑1 (MQ=255)                                                                                                                aaaaTGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTcggcg                                                                                                                                   <  1:144735/251‑1 (MQ=255)                                                                                                                aaaaTGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTcggcg                                                                                                                                   >  1:153850/1‑251 (MQ=255)                                                                                                                            tAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGGCTCtgggg                                                                                                                       >  1:515058/1‑248 (MQ=255)                                                                                                                                 gCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCAccc                                                                                                                  >  1:485034/1‑251 (MQ=255)                                                                                                                                       ttGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAgc                                                                                                            >  1:184666/1‑251 (MQ=255)                                                                                                                                       ttGGATATCTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGCCTGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGCGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGGTCGTTGTCCAATTCGGCGGCGACTCTGGGGCACGCGGAAgg                                                                                                            >  2:283139/1‑250 (MQ=255)                                                                                                                                                 agcTTAAAGTGAACAGTTGGAAGAATCGGTTGGGCTTCCAGCCCGAGCGCGGGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGGCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACACGGAAGCGCATTTTGAc                                                                                                  <  2:29957/249‑1 (MQ=255)                                                                                                                                                     taacgTGAACAGTGGGATGAATCGGTTGGTCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACGCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCATGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCCCCCGGAAGCGCATTTTGACt                                                                                                 <  2:212158/244‑1 (MQ=255)                                                                                                                                                     tAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTg                                                                                                <  2:153850/249‑1 (MQ=255)                                                                                                                                                      aaaGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTggg                                                                                              <  1:31847/250‑1 (MQ=255)                                                                                                                                                            aaCAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGACTGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGCGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGGGGCACCCGGGAGCGCATTTTGTCTGGGTTCGGcc                                                                                       >  2:125641/1‑251 (MQ=255)                                                                                                                                                             aCAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCt                                                                                      >  1:83786/1‑251 (MQ=255)                                                                                                                                                                       gAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTAttt                                                                            >  2:567860/1‑251 (MQ=255)                                                                                                                                                                        aaTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAAttt                                                                                                                                                                                                          >  1:234610/1‑124 (MQ=255)                                                                                                                                                                        aaTCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAAttt                                                                                                                                                                                                          <  2:234610/124‑1 (MQ=255)                                                                                                                                                                                         gcccggtggcgTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACACGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCgg                                                          <  2:4649/243‑2 (MQ=255)                                                                                                                                                                                                cgcgTGCTTACCGGCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATGCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTc                                                    <  2:413996/250‑1 (MQ=255)                                                                                                                                                                                                 gcgTGCTTACCGTCTGGCAGCAGGTGCGGGCAATGGCGAATGGTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACGGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTTCGGTCAg                                                  <  1:66655/251‑1 (MQ=255)                                                                                                                                                                                                        tACCGTCTCGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCCTTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTg                                           >  1:221216/1‑251 (MQ=255)                                                                                                                                                                                                            gTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGGGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGACGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGGGGCGTGctat                                       >  2:92808/1‑251 (MQ=255)                                                                                                                                                                                                              cgggCAGCAGTTGCGGGTAATGGCGAATGGTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGATGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCTTCCGGACAGGGGCGTTATATc                                      <  2:586619/248‑1 (MQ=255)                                                                                                                                                                                                                ggCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATAC‑TGATGGAATTTCCGGCGAGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTAGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGcc                                   >  1:135823/1‑251 (MQ=255)                                                                                                                                                                                                                ggCAGCAGGTGCGGGCAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCg                                     <  2:275547/249‑1 (MQ=255)                                                                                                                                                                                                                 gaagcagTTACGGGCAATGGCGAATGTTGTCGAAATGACTCTGATGTCGCATTTTGCCGAAACGGAACATCC‑TGATGGAATTTCCGGCGCGTTGGCGCGTATTGAGCCGGCGGCGGAGGGGCTGGCATGTCGGCGTTCGTTGTCCAGTTCGGCGGCGACTCTGTGGCACCCGGACGCGAATTTTGCCTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCGGGTCAGTGGCGTGatacctccc                                  >  2:94524/3‑245 (MQ=255)                                                                                                                                                                                                                    gcagTTGCGGGCAATGGCGAATGTTGGCGAAATGGCTCTGATGGCGCATTTTGCCGAAGCGGATCATCC‑TGATGGAATTTCCGGGGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATa                               <  2:473974/251‑1 (MQ=255)                                                                                                                                                                                                                       gTTGCGGGCAATGGCGAATGTTGGCGAAATGGCTCTGATTGCGCAATTTGCCGAAGCGGAACGTCCATG‑TGGAATTTCCGGCGCGATGGCGCGTATGGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGATCGTTGTCCAATTCGGCGGCGACTCTGTGGCTCCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATATTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATAGCGCCAATGCCg                            <  2:391218/251‑1 (MQ=255)                                                                                                                                                                                                                         tGCGGGCAATGGCGAATGGTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCTGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGTCAGCCGGATGCGCATTTTGACTGGGTTCGGCCTGGCATTAGTGTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCgg                           <  1:354621/250‑1 (MQ=255)                                                                                                                                                                                                                               cAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGt                    <  1:567860/251‑1 (MQ=255)                                                                                                                                                                                                                               cAATGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGt                    <  1:72405/251‑1 (MQ=255)                                                                                                                                                                                                                                  tGGCGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCg                 >  1:470829/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCGAATGTTGGCGAAATGACTCTGATGTCGCCTTTTGCCGAAGCGGAAAATCC‑TGATGGAATTTCCGGCGAGATAGCGCGTATTGAGCAGGCGGCGGAGGGGATGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGAGATTCGCCGTCCGGTCCATGGCGGGATATCGCCAATAACGGCTTACGTACggg               >  2:370069/1‑250 (MQ=255)                                                                                                                                                                                                                                     cGAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGTGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGGCGGCGTTCGTTGTCCAATTCGGCGGCGACGCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGtg              <  1:335915/251‑1 (MQ=255)                                                                                                                                                                                                                                      gAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGGTGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCt                                                                                      >  2:576018/1‑178 (MQ=255)                                                                                                                                                                                                                                      gAATGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCt                                                                                      <  1:576018/178‑1 (MQ=255)                                                                                                                                                                                                                                         tGTTGGCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCCGTGGCGTGATATCGCCAATACCGGATTACGTCCGGtgatg           >  1:539037/1‑250 (MQ=255)                                                                                                                                                                                                                                          gTTGGCGAAATGCCGCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGGGCGATGTCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCCCCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGTATTACGTCCGGTGATGAc         <  2:529039/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑GGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCt      <  2:63932/251‑1 (MQ=255)                                                                                                                                                                                                                                              gCGAAATGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCt      >  1:42780/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aaTGACTCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGc  >  1:428294/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     AAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCGGTTGGGCTTCCAGCCCGATCGCGTGCTTACCGTCTGGCAGCAGTTGCGGGCAATGGCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCC‑TGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTCGGCGGCGACTCTGTGGCACCCGGAAGCGCATTTTGACTGGGTTCGGCCTGGCATTATTTTGTATGGCGCTTCGCCGTCCGGTCAGTGGCGTGATATCGCCAATACCGGATTACGTCCGGTGATGACGCTAAGC  >  NC_000913/1239096‑1239586

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