Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	97.7%	139.1 / NA	43	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/22);  minor base A (1/0);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993285‑3993767                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             gCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTa                                                                                                                                                                                                                                          <  2:158711/251‑1 (MQ=255)     aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTa                                                                                                                                                                                                                                      <  1:240226/251‑1 (MQ=255)       cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACg                                                                                                                                                                                                                                    <  1:247885/251‑1 (MQ=255)       cAAACTACACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTAACATTACCGGCGGCGGAAGACGGCGTTTCCCGCGGAGGGATCATTCAGGTCTTTTTCGAAGAAACACCAAAAGAGACAGGGTTTACTATCTACCTTCTCGACGCAAGCAACCGGGGTTG‑GGGATCACAACACTGGGAAGCCAGCAAAAAAGAGATAGAGACG‑GACGTCAGTCGCTTCTACACGTCTTCGCAAGAACGCTTTACCTAcc                                                                                                                                                                                                                                    >  1:170894/1‑250 (MQ=255)                cGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGGAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGGATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCCGTCGCTTCTACTCGTCCTCGCATGACCGCTTTACCTACGGCTCAAGCt                                                                                                                                                                                                                           >  1:175098/1‑251 (MQ=255)                              aGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAAc                                                                                                                                                                                                             <  2:210308/251‑1 (MQ=255)                                gTTGAAACCAATCACGTCAAATTACCGGCGGTGGGGGACGGCTTTGCCAGCGAAGGGGTCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACCTTCTCGACGACAGCAACCGGG‑TTGCGGGATATCACCACTGCGCAGGCAGCAAAGAGGGGCT‑GGTACGTGACGTCAGTTGCTTCTAATCGTCATCGCATGACCGCTTTTCCTACGGCCCCAGCTTCATCCACTTCacccc                                                                                                                                                                                                           >  1:22468/1‑247 (MQ=255)                                  tGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTg                                                                                                                                                                                                          <  1:192/250‑1 (MQ=255)                                    aaaCCAATCACGTCAAATTACCGGCGGTGGTGGTCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCg                                                                                                                                                                                                       <  2:82502/251‑1 (MQ=255)                                           tCACGTCAAATTACCGGCGGTGGGGGACGGCTTTGCCAGCGCAGGGATCAACCCGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTACTATCTACAATCTCGACGAAAGCCACCGGG‑TTGAGGGATATCAACACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGGCCCTTCCAATCGTCCTCGCATGACCGATTTACCTACGGCCCAAGCTTCATCAACTTCCACCTGCCGcccgtct                                                                                                                                                                                                >  1:270139/1‑245 (MQ=255)                                                         cggcggTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCATCCTGCCGCAGTTCTATCAGATTGTGAAg                                                                                                                                                                                  <  1:267745/251‑1 (MQ=255)                                                                       ggCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTg                                                                                                                                                                     >  1:248536/1‑250 (MQ=255)                                                                                    aaGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGCAGGTTGATGGTCGTGAACAGGTGAttccg                                                                                                                                                       >  2:91798/1‑251 (MQ=255)                                                                                         atcatcCAGTTCTTTTTCGTAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgttccg                                                                                                                                                  >  2:67122/1‑251 (MQ=255)                                                                                            atcCAGTTCTTTTTCGAAGCAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACAGCTTTACCTACGGCTCACGCTTCATCACCTTCAACCTGCCGCAGTTATATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGCTTCCGTTCCGccc                                                                                                                                               >  2:118186/1‑249 (MQ=255)                                                                                             tcCAGTTCTTTTTCGAAGAAACGCAAGTCGCGCATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGAATATCACCACTGCGATGGCAGCAAAGAGGAGCT‑GGTACGCGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTAATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGGGATACCGTTCCGcac                                                                                                                                               >  2:82826/1‑250 (MQ=255)                                                                                                          cGAAGAAACGCAAGACGAGAATGACTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGt                                                                                                                                  >  2:24456/1‑250 (MQ=255)                                                                                                                    cAAGACGAGAATGGCTTTAATATCTACATTCACGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAATGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgccgcc                                                                                                                       <  2:128867/251‑1 (MQ=255)                                                                                                                                    ttAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCa                                                                                                        <  2:35378/250‑1 (MQ=255)                                                                                                                                         atCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATa                                                                                                   <  2:224720/250‑1 (MQ=255)                                                                                                                                             aCATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                                                                                  <  2:175321/247‑1 (MQ=255)                                                                                                                                                     gacgaAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTAc                                                                                       >  2:102371/1‑250 (MQ=255)                                                                                                                                                        gaAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGc                                                                                    >  2:68157/1‑250 (MQ=255)                                                                                                                                                                      ttGAGGTCTATCACCAATGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtg                                                                       >  2:276868/1‑250 (MQ=255)                                                                                                                                                                           gTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAAc                                                                 <  1:130234/251‑1 (MQ=255)                                                                                                                                                                               atCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGc                                                             >  1:186386/1‑251 (MQ=255)                                                                                                                                                                                     acTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGaaa                                                       <  2:47166/251‑1 (MQ=255)                                                                                                                                                                                                agcaAAGAGGAGCT‑GGGACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAAGCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAt                                            <  1:229998/251‑1 (MQ=255)                                                                                                                                                                                                 gcaAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                                           <  1:3842/251‑1 (MQ=255)                                                                                                                                                                                                 gcaAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCCTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                                           >  1:48500/1‑251 (MQ=255)                                                                                                                                                                                                      gaggagCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGt                                                                        <  1:255646/217‑1 (MQ=255)                                                                                                                                                                                                      gaggagCA‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGATCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGt                                                                        >  2:255646/1‑217 (MQ=255)                                                                                                                                                                                                         gagCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGc                                   >  2:153426/1‑251 (MQ=255)                                                                                                                                                                                                          ggCT‑GGTACGTGACGTCATTCGCTCCTACTCGTCATCGCATGACCGCTTTACCGACGGCTCAAGCTTCATCAACTTCACCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCTTTCCGCACAAAATCTAGCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCa                                  <  2:139797/250‑1 (MQ=255)                                                                                                                                                                                                          agCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCa                                  <  1:230682/251‑1 (MQ=255)                                                                                                                                                                                                           gCT‑GTTACGTGACGTCAGTCGCTTCTAAACGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTACATCAATGTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAAGAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCTGCCTGCTAAACAGGATCACAATACGCCGCTATTACAGCAATATTTATCTTGATGAACGTGCCGTAAAGCGCGTCTTATCCGGCATGCat                                 <  2:264470/251‑2 (MQ=255)                                                                                                                                                                                                                gtccGTGACGCCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCtt                             <  2:170256/248‑1 (MQ=255)                                                                                                                                                                                                                 tACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGACTCAACCTTCATCAACTTCAACCTGCCGCAGTTCTCTCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGAAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCtt                             >  1:19671/1‑250 (MQ=255)                                                                                                                                                                                                                               cGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgact               <  2:206670/250‑1 (MQ=255)                                                                                                                                                                                                                                 cTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt            >  1:128400/1‑251 (MQ=255)                                                                                                                                                                                                                                     tACTCGTCATCGCATGACCGCTTTACCTACGGATCAAGCTTCATCACCTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCAGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAc        >  1:28471/1‑251 (MQ=255)                                                                                                                                                                                                                                         cGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGAAATCTTAGCGGAAGCTGACTGTTTCACCCGc    >  1:250970/1‑251 (MQ=255)                                                                                                                                                                                                                                            cATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct  <  1:67122/250‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             GCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGG‑TTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCT‑GGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993285‑3993767

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