Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,394	G→T	L372L (CTG→CTT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,394	0	G	T	97.5%	139.0 / ‑3.7	40	L372L (CTG→CTT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (12/27);  total (12/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

ATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCT  >  NC_000913/109152‑109636                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     atGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGAt                                                                                                                                                                                                                                              <  1:122759/251‑1 (MQ=255) atGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGGCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGAt                                                                                                                                                                                                                                              <  2:41550/251‑1 (MQ=255)       gAAGGGGAGTCTCTGTCCTCTCCGCCCAACATCATGCTGCTGCACCACGTAACGCCGGCGCTGCGCGCTCACGCCCTGTTTACCCGTGGCGTCGACTACCTCGTTAAAGATGGTGAAGTTCTCATCGTTGACGAACACACCGGTCGTAACATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCCGGCTGTGGACGAGCACGAAGGTGTGCAGATCCCGACCGAAAACCACAAGCTTGCTTCGATCACCtt                                                                                                                                                                                                                                        >  2:38491/1‑251 (MQ=255)                       aCTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGACCTACTTCCGt                                                                                                                                                                                                                        >  2:148701/1‑251 (MQ=255)                                    catcatGCTGATGCACCACGTAACAGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCAACAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCAAGAACTACTTCCGTCTGTATGAAAAAc                                                                                                                                                                                                           >  1:451765/1‑251 (MQ=255)                                       catGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGTCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTg                                                                                                                                                                                                         <  1:278714/250‑1 (MQ=255)                                       catGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTg                                                                                                                                                                                                         <  2:375876/250‑1 (MQ=255)                                                         aaCGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTTTTTATGAAAAACTGGCGGGGATG‑ACCGGTAct                                                                                                                                                                                      <  1:415029/250‑1 (MQ=255)                                                           cggcggCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTActgct                                                                                                                                                                                   <  1:173348/251‑1 (MQ=255)                                                              cggcgCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGACGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTGCTGCGGAt                                                                                                                                                                                >  2:457090/1‑251 (MQ=255)                                                                 cgcTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTTCTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATAcc                                                                                                                                                                             <  1:267948/251‑1 (MQ=255)                                                                     gcgcgcTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAg                                                                                                                                                                         <  1:178665/251‑1 (MQ=255)                                                                      cgcgcTCATGCGCTTTTTACCCGTGACGTCGACTACAGCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGc                                                                                                                                                                        <  1:372351/251‑1 (MQ=255)                                                                                        aCCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAAt                                                                                                                                                      <  1:295174/251‑1 (MQ=255)                                                                                         cccGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATc                                                                                                                                                     <  1:356311/251‑1 (MQ=255)                                                                                                  tCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTg                                                                                                                                            >  1:26486/1‑251 (MQ=255)                                                                                                        aCATCGGTAAAGATGGTGAAGTTAGCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATAcc                                                                                                                                      <  2:429553/251‑1 (MQ=255)                                                                                                              ttAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACAGCTGATACCGAAGCTTTTGAATTTAGCTCAATCTAACAGCTGGATACCGACgtt                                                                                                                                >  2:141550/1‑251 (MQ=255)                                                                                                               tAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCgtt                                                                                                                                <  2:136342/250‑1 (MQ=255)                                                                                                                    aTGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTTGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCg                                                                                                                         <  2:145880/251‑1 (MQ=255)                                                                                                                      ggTGAAGTTAGCATCGTTGACGAACACACCGGGCGGACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGAGG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCGCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGAc                                                                                                                       <  2:52524/251‑1 (MQ=255)                                                                                                                         gAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGGTG‑ACCGGGACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACCAGCTGGATAACGTCGGT‑GTTcacacca                                                                                                                     >  1:183770/1‑244 (MQ=255)                                                                                                                                       gtccgaaccccctgtgcggtccgtGCAGGGCCGTCGCGGGTCCGAGGGTCTGCCCCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATTCAGACCGAAAACCAAACGGTTGCTTCGATCCCCTTCCAGAACTACTTCCGTCTGTATGAAAAATTGGCGGGGATG‑ACCGGTTCTGCTGATACCGAAGATTTCGAATTTAGCTCAATCTACAAGCTGGAGACCGTCGTT‑GTTCCGCCCAACCGTCCAATGATTc                                                                                                      <  2:248208/228‑1 (MQ=255)                                                                                                                                                      gcGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCgg                                                                                       <  2:190010/250‑1 (MQ=255)                                                                                                                                                        gTACCATGCAGGGCAGTCGCTGGTCCGATGGTCTGCACCAGGATGCGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATGGACCG‑TACTGCGGTTAACCAAACTTTCGAATTTAGCTCAATACACAAGAGGGATACCGTCCGT‑GTTGCGGCCCACCGTCCAATGATTcagaacaagctgccgga                                                                                      >  2:320252/1‑234 (MQ=255)                                                                                                                                                               gCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCt                                                                              <  1:148701/251‑1 (MQ=255)                                                                                                                                                               gCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCt                                                                              <  1:457090/251‑1 (MQ=255)                                                                                                                                                                         cGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGa                                                                    <  2:26486/251‑1 (MQ=255)                                                                                                                                                                             ggTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGc                                                                 <  1:31608/250‑1 (MQ=255)                                                                                                                                                                                          gCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGACCTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATTCCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCat                                                    <  2:89230/250‑2 (MQ=255)                                                                                                                                                                                          gCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTc                                                      <  1:348765/248‑1 (MQ=255)                                                                                                                                                                                          gCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAgg                                                   >  1:84262/1‑251 (MQ=255)                                                                                                                                                                                              cAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGAt                                               <  1:1573/251‑1 (MQ=255)                                                                                                                                                                                                                  agGTGTGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATAACAAAGAAc                           <  2:124978/251‑1 (MQ=255)                                                                                                                                                                                                                    gtgtGCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTCCTTCCGTCTGTATGAAAAACTGGCGGTGATG‑ACCTGGACTGCTGATACCCAAGCTTTCGAATTTAGCTCCATCTACAAGCTGGGTAACCTCGTTTGTTCC‑ACCAACCGTACAATGATTCGTCAAGATCTGCCGGACCTGTGCTACATGCCTGAAGCGGAAAAAATTCAGGCGATCCTTGACGCTATCacataacg                          >  2:172010/1‑243 (MQ=255)                                                                                                                                                                                                                        gCAGATCCAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGGCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTAccg                     >  2:245894/1‑249 (MQ=255)                                                                                                                                                                                                                        gCAGAGCCAGAACGAAAACCAAACGCTTGCTGCGATCACCTTCCAGAACTACTTCCGGCTGTATGAAAAACGGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGGCCAATGATTCGGAAAGATCGGCCTGACCTGGTCTACATGACTGAAGCGGAAAAAATGCAGGCGATCAGTGAAGATATCAAAGTACGTACt                      <  2:56275/250‑1 (MQ=255)                                                                                                                                                                                                                               cAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAgg              <  1:427128/251‑1 (MQ=255)                                                                                                                                                                                                                               cAGAACGAAAACCAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCAATGAAGATATCAAAGAACGTACTGCGAAAgg              >  2:364138/1‑251 (MQ=255)                                                                                                                                                                                                                                           cAAACGCTTGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCt  <  1:245894/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     ATGGATGAAGGGGAGTCTCTGTACTCTCCGGCCAACATCATGCTGATGCACCACGTAACGGCGGCGCTGCGCGCTCATGCGCTGTTTACCCGTGACGTCGACTACATCGTTAAAGATGGTGAAGTTATCATCGTTGACGAACACACCGGTCGTACCATGCAGGGCCGTCGCTGGTCCGATGGTCTGCACCAGGCTGTGGAAGCGAAAGAAGGTGTGCAGATCCAGAACGAAAACCAAACGCTGGCTTCGATCACCTTCCAGAACTACTTCCGTCTGTATGAAAAACTGGCGGGGATG‑ACCGGTACTGCTGATACCGAAGCTTTCGAATTTAGCTCAATCTACAAGCTGGATACCGTCGTT‑GTTCCGACCAACCGTCCAATGATTCGTAAAGATCTGCCGGACCTGGTCTACATGACTGAAGCGGAAAAAATTCAGGCGATCATTGAAGATATCAAAGAACGTACTGCGAAAGGCCAGCCGGTGCT  >  NC_000913/109152‑109636

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