Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,189	C→T	intergenic (‑1136/‑1408)	topAI ← / → ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,535,189	0	C	T	98.4%	190.2 / NA	64	intergenic (‑1136/‑1408)	topAI/ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/39);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCG  >  NC_000913/4534942‑4535432                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      actGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTAAATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGGACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAgg                                                                                                                                                                                                                                                    >  2:525354/1‑251 (MQ=255)          cACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCt                                                                                                                                                                                                                                            <  1:64374/250‑1 (MQ=255)          cACTTTGCTATCACGGTAAGTGGCAGACTAAAAGAGGGATGGAAAAAAGACGATTTTCCTCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCAAAATGGAGGTAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCACGGCGGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCt                                                                                                                                                                                                                                            <  1:553528/250‑1 (MQ=255)            cTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATc                                                                                                                                                                                                                                         >  2:500571/1‑251 (MQ=255)                gCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGGCGCCATTTTCGGCGTAGGTTAATGCTATCaaa                                                                                                                                                                                                                                      <  1:60949/250‑1 (MQ=255)                gCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCaaa                                                                                                                                                                                                                                      <  1:581380/250‑1 (MQ=255)                   aTCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑g                                                                                                                                                                                                                                  <  1:438410/250‑1 (MQ=255)                     cACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGTTTAATGCTATCAAAAT‑GTCa                                                                                                                                                                                                                               <  2:503300/251‑1 (MQ=255)                      aCGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCAt                                                                                                                                                                                                                              >  1:479306/1‑251 (MQ=255)                                 cAGCCTCAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGCAAGCTg                                                                                                                                                                                                                   >  2:501141/1‑251 (MQ=255)                                    ccTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTCCAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGa                                                                                                                                                                                                                 >  2:252911/1‑250 (MQ=255)                                      tAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAAt                                                                                                                                                                                                              >  1:575661/1‑251 (MQ=255)                                        aaaGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTa                                                                                                                                                                                                            <  2:197324/251‑1 (MQ=255)                                         aaGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTaa                                                                                                                                                                                                           <  2:473760/251‑1 (MQ=255)                                                agggAAAAAGGACGATTTGCCCCGCAAGGGGATATTTACATAGTGATGGTAGAGGCAGGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGTCGGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGGGTTGTGACGCCATTTTCGGCGTAGGTTAATTCTAACAAAAT‑GGTATGCCGTAAGCTGGAAATTAATGGCaaa                                                                                                                                                                                                    <  2:357984/249‑1 (MQ=255)                                                aTGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGCAGCCAAATGGCGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCaa                                                                                                                                                                                                     >  2:175081/1‑250 (MQ=255)                                                   gAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGAc                                                                                                                                                                                                 <  1:243156/251‑1 (MQ=255)                                                                tgtCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGGCATTACGGGGAACCACTGGACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTa                                                                                                                                                                                     <  1:534136/248‑1 (MQ=255)                                                                     ccGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGATGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATTGCTGTTGGCACATAAGGATGATATAGAGCGCCATTACTGGTAACCTCTGTACAAAGTTGTTGTTGTTGTGACGCCATTTTCGGCGTAGTTTAATGCTATCAAAAT‑GTCTTTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACa                                                                                                                                                                              <  1:326392/251‑1 (MQ=255)                                                                      cGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAg                                                                                                                                                                             <  2:350801/251‑1 (MQ=255)                                                                        cAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTg                                                                                                                                                                           <  1:35731/251‑1 (MQ=255)                                                                        cAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTg                                                                                                                                                                           >  1:442426/1‑251 (MQ=255)                                                                        cAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTg                                                                                                                                                                           >  1:29588/1‑251 (MQ=255)                                                                                            gATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGc                                                                                                                                                        >  2:190961/1‑250 (MQ=255)                                                                                                      gCATGAAGCCAAAGGGAGGGAGCCGTCATAGCCCTCTTGAAGCGGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGaa                                                                                                                                             <  2:458138/251‑1 (MQ=255)                                                                                                                    agagggagCAGTCATAGCCCTCTTTAAGATGAAACACTAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTCAGGGGACCCACTGTACAAAGTTGTTGGGGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCACACT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGGGCAGGGTTTACTG‑ACAGTGGGGACTGTGCCCTCCTTGCAGGGGGAGGAAAGGGCGATTGTtct                                                                                                                               >  2:226550/2‑251 (MQ=255)                                                                                                                        ggagCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCttt                                                                                                                           >  1:561786/1‑251 (MQ=255)                                                                                                                             cGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTcc                                                                                                                      <  1:252911/251‑1 (MQ=255)                                                                                                                                cATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGt                                                                                                                   <  1:500571/251‑1 (MQ=255)                                                                                                                                 atagcccTCTTGAAGCTGAAACCATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGCGCGCCATTCCGGGGAACCACTGTACAAAGATGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGG‑CGGAAAGGGCGGATGTTCCTTTTTCTCCttgct                                                                                                                 >  2:124527/5‑247 (MQ=255)                                                                                                                                   aaCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATACGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGc                                                                                                                                                        >  1:413865/1‑211 (MQ=255)                                                                                                                                   aaCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCAGTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGc                                                                                                                                                        <  2:413865/211‑1 (MQ=255)                                                                                                                                    aCCCTCTTGAAGCATAAACAATAGCAACGGTGCATGTGGCACATAAGGGTGATATATAGCGCCATTTCGGGGAACCTCTGTACAAAATTGTTGGGGTTGTGACGCCAATTTCAGCGTAGGGTTATTCTATCAAAAA‑TTCATTTCGTAAGCTGGGAATTAATGGCAAAGACGAGCGGGGTTTTCTG‑ACAGTGGGGACTGTGCACTCCGTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTAc                                                                                                               <  1:583649/251‑1 (MQ=255)                                                                                                                                     cccTCGTGAAGCTGAAACAATAGCAGCAGGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACt                                                                                                              <  2:598084/251‑1 (MQ=255)                                                                                                                                          ttGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTT‑CGGCGGAGGTTAATGCTATTAAAATTGT‑ATTGCGGAAGCTTGGAATTAAAGGCAAAGACGAAGAGGGGTTACTGGAAAG‑GGGGGCTGTGGACTCCTTGGAGGGGGGGGGAAGGGGGGTTGGTTTCTTTTTTCCTGTCTACTCaaaa                                                                                                         >  2:567794/1‑250 (MQ=255)                                                                                                                                             aaGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACAt                                                                                                      <  1:42687/251‑1 (MQ=255)                                                                                                                                               gCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGACAGGGCGATTGTTCTCttagatcgaagagcgtcgtgtaggg                                                                                                     >  2:315319/1‑227 (MQ=255)                                                                                                                                               gCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCtt                                                                                                                            <  1:315319/227‑1 (MQ=255)                                                                                                                                                   aaaCAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGAt                                                                                                <  1:143422/251‑1 (MQ=255)                                                                                                                                                       aaTAGCAGCATGGCTGTTGCCACATAAGGATAATATAGAGCGCCATTACGGGGAACCCCTGTACAAAGTTGTTGGTGTTGTGACGCC‑TTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAg                                                                                           <  1:567794/251‑1 (MQ=255)                                                                                                                                                        aTAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAg                                                                                           <  1:381506/251‑1 (MQ=255)                                                                                                                                                          agcagcATGGCTGGTGTCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTTCAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGTGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGt                                                                                         <  1:41386/251‑1 (MQ=255)                                                                                                                                                            cagcaTGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCa                                                                                                                                                                                            <  1:339459/151‑1 (MQ=255)                                                                                                                                                            cagcaTGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCa                                                                                                                                                                                            >  2:339459/1‑151 (MQ=255)                                                                                                                                                            cagcaTGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCTATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGttt                                                                                       <  1:210827/251‑1 (MQ=255)                                                                                                                                                                         ggCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAAc                                                                           >  1:470065/1‑250 (MQ=255)                                                                                                                                                                          gCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAA‑T‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATTGCAGGTTTCTTGATAGCAACAGc                                                                        <  1:171977/251‑1 (MQ=255)                                                                                                                                                                              aTAAGGATGATATAGATCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACt                                                                     <  1:117626/251‑1 (MQ=255)                                                                                                                                                                                    atgatATAGAGCGCCATTACGTGGAACCACTGTACAAAGTTGTTGGTGTTGTGACTCCATTTTCGGCGTAGGTTAAGGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTGTACTG‑ACAGTGGGGACTGTGCACTCCTTTCAGGGGGCGGAAAGGGCGATTTTGCTCTTTTCTCCTGTCTACTCAGAACATGAAGTTGGCAGTTGTCTTGATAGCAACAGCACTAACCTc                                                               <  1:600100/251‑1 (MQ=255)                                                                                                                                                                                           agagCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTCGGTTAATGATATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTg                                                        >  2:96208/1‑251 (MQ=255)                                                                                                                                                                                                  cATTACGGGGAACCACTGTACAAAGTTGTTGGTTTTGTGACGCAATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTTCTTGCAGGGGGCTGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCTCTATCCTCAATGTTGCTGtctc                                                 <  1:627776/251‑1 (MQ=255)                                                                                                                                                                                                     tACGGGGAACCTCTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAc                                               <  2:353437/250‑1 (MQ=255)                                                                                                                                                                                                       cGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTGCGGCTTAGTTTAATGCGATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGTACTGTGCACTCCTTGCAGGGGTCGGAAAGGGCGATTTTTCTCTTTTCTGCTGTCTACGCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAcgcg                                            <  1:394346/251‑1 (MQ=255)                                                                                                                                                                                                          ggAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTa                                         <  2:29588/251‑1 (MQ=255)                                                                                                                                                                                                           gAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTa                                         >  1:368757/1‑250 (MQ=255)                                                                                                                                                                                                              ccACTGTACAAAGTTGTTGGAGTTGTGACGCAATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGCCGAGCAGGGATTACTG‑ACAGTGGGGACTGTGAACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTATCTTTTATCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTaa                                        >  1:231266/1‑248 (MQ=255)                                                                                                                                                                                                               cACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTTCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGAt                                    <  1:236628/251‑1 (MQ=255)                                                                                                                                                                                                                             gttgGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTAtt                      >  2:317935/1‑251 (MQ=255)                                                                                                                                                                                                                                     tgtgACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCATCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGat               <  1:92606/250‑1 (MQ=255)                                                                                                                                                                                                                                        gACGCCATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATgg           <  1:501141/251‑1 (MQ=255)                                                                                                                                                                                                                                             cATTTTCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCtt      >  2:395/1‑251 (MQ=255)                                                                                                                                                                                                                                                tttCGGCGTAGGTTAATGCTATCAAAAT‑GTCAGTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTAt    <  1:590992/250‑1 (MQ=255)                                                                                                                                                                                                                                                 ttCGGCGTAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCg  >  1:150262/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ggCGTAGGTTAATGCTATCAAA‑T‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑AAAGTGGGGAATGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATTAAGATGGCAGGTTTCTTGATAGCCACAGCACTATCCTCAATGTTGCTGTCCCACGCGGTAAAAGTAGTTTTTTGGTATTTGGCGATTTGGACCTTATCg  >  1:336303/1‑247 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ACTGTAACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAAT‑GTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTG‑ACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCG  >  NC_000913/4534942‑4535432

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