Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,780	T→C	T349T (ACT→ACC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,780	0	T	C	91.7%	84.5 / ‑1.2	37	T349T (ACT→ACC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (15/18);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAA  >  NC_000913/4535‑5020                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cGTTTTATTGATGCGACCAACGTGAAGGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGGGTTGTTCCGCCGCAAAATCTGGCACCTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGTCAATGCGTGAGTTACAAGAACTGGGCTACACCTCgg                                                                                                                                                                                                                                                >  1:116409/1‑250 (MQ=255)       aTTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGc                                                                                                                                                                                                                                          >  1:289628/1‑250 (MQ=255)       aTTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCGGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGc                                                                                                                                                                                                                                          <  1:303599/250‑1 (MQ=255)                  caacgtgaagtaagcggTGCCACGTGTCCGGCACGACGGTCAGTGGTCACTAAAAGCGACTAAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCAGCG‑TGTGGAAGAGGTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAg                                                                                                                                                                                                                              <  1:197575/235‑1 (MQ=255)                  cAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATACAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATATGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGt                                                                                                                                                                                                                                >  1:233262/1‑249 (MQ=255)                      gTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGcgt                                                                                                                                                                                                                           >  1:220562/1‑248 (MQ=255)                        aaccgacaCCGTCCCACGTTCCCTGAACGAGGGTCAGGGGTCCCCCAAAGCGCATAAGGAGACGTTATCCAACGCGAGGGAGGGTAGCCAGCCGCACAGCGGTCCGCGATG‑GGAAGATTTGGTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTGAACCTCGGAGCCGCGCGCTGCCGGAGCGTATc                                                                                                                                                                                                                        <  2:188350/244‑1 (MQ=255)                           cGATCCCGTTCCACGTTTCCTGCCCGACGGTCAGTGGTCACCCAAAGCGACGCAGGCGACGTTATCAAACGCGATGGCCTTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTg                                                                                                                                                                                                                     <  1:47416/251‑1 (MQ=255)                             aTACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTgc‑g                                                                                                                                                                                                                  >  1:152739/1‑251 (MQ=255)                                     ccaagTTTCCGGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGGCGTTATCCAACGAGATGGTCGTGAGTCAGTCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGTTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCGACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GATGCGTg                                                                                                                                                                                                           <  1:343431/246‑1 (MQ=255)                                                    gacgGTCAGTGATCACCCAAAACGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTACGCCGCAAAATCTGGAAACTGAAAGAGCTGGGATATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGCACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCGGCGTGACCAGGTGAATCCAgg                                                                                                                                                                                           >  2:370674/1‑251 (MQ=255)                                                     acagcgcgagggaccccAAAACGCCTCGGGCGACTTTATACACCGCTAAGGACGTGAGTCAGCTAAACAACTGGCCGCG‑TGGGAAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCTGCACGCTGCCGTAGCTTATCGTGT‑GCTGCGTGACCAGTTGAATCCAGGc                                                                                                                                                                                          <  2:200407/237‑1 (MQ=255)                                                                 cccccAAAGCGACTCAGGCCACGTTATCCAACGCGTTGGCCGTGATTCGGCCGAACACCTGGCCGGT‑TGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGGGAATATGGCTTg                                                                                                                                                                              <  2:222125/249‑1 (MQ=255)                                                                        aGCGACTCAGGCGGCGTTATCCAACGCGAGGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                       <  1:101714/251‑1 (MQ=255)                                                                                     aCGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATcc                                                                                                                                                         <  1:379567/251‑1 (MQ=255)                                                                                      cGTTATCCAACGCGATGGACGTGAGTCAGCCGAACCACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGAGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGAGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTGCGCACC‑GCGCCTCCg                                                                                                                                                        >  2:186860/1‑251 (MQ=255)                                                                                                cgcgATGGCCGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAAttt                                                                                                                                               >  1:9081/1‑250 (MQ=255)                                                                                                           tGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGGAATTTAAagagggagtg                                                                                                                                   >  1:197684/1‑245 (MQ=255)                                                                                                                      aacaacGGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtct                                                                                                                        <  1:141914/251‑1 (MQ=255)                                                                                                                        caacTGGCAGCG‑TGTGGAAGCGTTGTTCCGCCGCAAAAACTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATACAGGCGCATATGGAGTGTTCCTACGCACCTCCGC‑CCCGGCGAAATTTTAAGAGAGCGGGTAAGCGGTtcgcg                                                                                                                      >  2:35345/1‑248 (MQ=255)                                                                                                                                ccggtgggggAAGAGTGATTCCGCCGCAAAATCTGGCAACGGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACGCCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc                                                                                                               <  2:118377/243‑1 (MQ=255)                                                                                                                                       tggcagagTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGTCGCACGCGGCCGAAGAGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCATCGGCACC‑GCGCATCCGGCGAAATTTAAAAAGGGGGTGGTATCGATTCCCGGTGGAACGtaggtt                                                                                                        >  2:85658/5‑246 (MQ=255)                                                                                                                                       gggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                        <  1:216401/250‑1 (MQ=255)                                                                                                                                        ggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATc                                                                                                       >  1:262021/1‑251 (MQ=255)                                                                                                                                                             aaTCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTTAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATTCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCCAAAAGAGCTGGCAGAAc                                                                                  <  1:184688/251‑1 (MQ=255)                                                                                                                                                                 tGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGc                                                                              >  1:166399/1‑251 (MQ=255)                                                                                                                                                                         gAAAGAGCTGGGTTATGCAGACGTGGGTGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGAttt                                                                        <  1:63037/249‑1 (MQ=255)                                                                                                                                                                           aaGAGATGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAACAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGGGTAACGTTGGATCTGCCAAAAGAGCTGGAAGAACGTGCTGATTTAccc                                                                    >  1:104402/1‑251 (MQ=255)                                                                                                                                                                                       atgcagcagtgtctggtgacaCCATGCAACATAAAATGCGTGAGTTAACAGAACTGGGCTACACATCGAAGGCGCACGATTCCGCAGCATATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGCATATGGCTTGCTCCTCGGCACC‑GCGCAGCAGGCGATATTTACAGAGAGCGTGGAAACGCTTCCCGGTGTATCGTTGGATGTGCCAAACGACCTGGCAGAACGTGCTGATTTACCCTTGCtttacaat                                                        >  2:32161/21‑246 (MQ=255)                                                                                                                                                                                                gggtgtatgaAACCACGCAACAGACAATGCGTGAGCTCAAAGACCTGGCATACACCACAGAGCCGCCCGCTGCTGTAGCGTATCGTGC‑GCTGCGTGACAAGTTGAAATTAGGCGCATATGGCGTGTTCCTCGTCCCC‑GCTCATCCGGCGAAATATAAAGAGAGCGTGTAAGCGAGTCTCGGTGAAACGTTGGAGCTGCCAAAAGAGCTGGCAGAACGTGCTGATGTTCCCTTTCTGTCACATAATCCGccc                                               <  2:295107/245‑1 (MQ=255)                                                                                                                                                                                                      atgaAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccgccggt                                         >  1:331417/1‑249 (MQ=255)                                                                                                                                                                                                                     aGACAATGCGTGAGATAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgcgt                          >  1:303320/1‑251 (MQ=255)                                                                                                                                                                                                                       aCAATGTGTGATTTAAAATAACAGGGCTACACCTCGGAGCCTCGCTCTGCCTTAGTGTGTCGTTC‑GCTGCGTGACTAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACC‑TCGCATCCGGCTAAATGTAAAGAGAGCGTGGAATCGATTCTCGGTGAAACTTTGGATCTGCAAAAAGAGCTGGCAGAACGTGCTTTTTTCCCCTTGTTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTa                         <  1:215049/250‑1 (MQ=255)                                                                                                                                                                                                                         aaTGCGTGAGGAAAAAGAACTGGGCTAAACCGCGGAGCGGCACGCTGCCGTAGCGTATCGGGG‑GATGGGGGACGAGTTGAATCCAGGCGAATAGGGCGTGTTCCTGGGCACG‑GGGCAAACGGAGAAATTTAAAGAGAGCGGGGAAGCGATTCGCGGGGAAACGTGGAATATGCAAAAAGAGAGGGCAGAACGGGCTGATTTCCCCTTCCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAt                      <  1:68793/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaaaGAACTGGGCTACGCCTCGGAGCCGCACGCGGCCGTAGCGTTTCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGAAGAACGTGGTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatc          <  1:34778/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaCTGGGCTACTCTTGGGAGCCGCAGGCTGCGGTAGCGTTTC‑TGCCGCTGCGTCACCAGTTGACTCCAGGCGAATAGGGCTTGTTCCTCGGCCCC‑GCGCAGCCGGAAAAATTAAAAGAAAGTGTGGAAGCGATTCTCGGTGAAACGTTGGATCTCCCAAAAGAGCTGGCAGAACGTGCTGATTTCCCCTTTCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAAtcatca      <  1:307314/250‑1 (MQ=255)                                                                                                                                                                                                                                              gggCTACACCCCGGATCAGCACGGTGCCGTAGCGTACCGTGC‑GCTGCGTAACCAGTTGAAACCAGGCGAATATGGCTTGTTCCCCGGCACC‑GCGCTTCCGGCGAAATTAAAAGATAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGATCTGGCAGACCGGGCTAATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATTATGAATCATCAGTaa  <  1:338652/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAA  >  NC_000913/4535‑5020

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