Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	595,271	T→C	N58S (AAC→AGC)	cusR ←	DNA‑binding transcriptional activator CusR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	595,271	0	T	C	97.4%	125.9 / NA	39	N58S (AAC→AGC)	cusR	DNA‑binding transcriptional activator CusR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/24);  minor base A (1/0);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATG  >  NC_000913/595039‑595514                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    gACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCata                                                                                                                                                                                                                                    >  1:205534/1‑250 (MQ=255)            cGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCTCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCCCTCACGTCCGGCAGCATCATATCCAGGAtt                                                                                                                                                                                                                        >  2:120706/1‑251 (MQ=255)              gCAACCTGAAACGGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTAt                                                                                                                                                                                                                      <  2:290314/251‑1 (MQ=255)                aaCCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATc                                                                                                                                                                                                                     >  2:284458/1‑250 (MQ=255)                       accTGACTTTCGATAATCACCGCCGCGCCACGTCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCTAATGGCTTCACCAGGTAGTCATCTGCCCCGACCTCCAGCAGCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCTTCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGACGATTCTCAGATCATa                                                                                                                                                                                                             <  2:273616/249‑1 (MQ=255)                                 cGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTc                                                                                                                                                                                                   <  2:341253/251‑1 (MQ=255)                                            gccgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCgcg                                                                                                                                                                                                                            >  1:220864/1‑213 (MQ=255)                                                 cccGCGCAGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAGCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCAGGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCAAAAGCACCGTTCATCGCCAGATGGTAGc                                                                                                                                                                                   <  2:225080/251‑1 (MQ=255)                                                   cgcgCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGcc                                                                                                                                                                                  >  2:6842/1‑250 (MQ=255)                                                         gcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCg                                                                                                                                                                                               >  1:67655/1‑231 (MQ=255)                                                         gcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCg                                                                                                                                                                                               <  2:67655/231‑1 (MQ=255)                                                              aaTGTGCGCACCCGTGCCAGCAGTTCAGCAAAAGCGAATGTCTTTCCCAGGTGATCATCTGCCCCCAACCCCAGCCCCTTGACTCCATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCTGCCTCCCTTTATTGGCGGAGCGTAACAGGCGCACGATATTCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTAGCAGATCATTATCGCCGGTCATCGTCTGGTTGTAGCTATTGAGCCCGt                                                                                                                                                                       <  2:242821/250‑1 (MQ=255)                                                                               cagaagTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGCCGCGATGTTCACTGGTGCCCAGCGCGGCAAGCAACAGAATCGGCATCCCTTTATCGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTCTCAGACCATAATCACCGGCCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCCAATCGAcc                                                                                                                                                      >  2:361745/5‑250 (MQ=255)                                                                                       cAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCGGTTGTCGGCCAAATCGACCACAAAACCg                                                                                                                                             <  1:318996/251‑1 (MQ=255)                                                                                        aGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCgg                                                                                                                                            >  1:207589/1‑251 (MQ=255)                                                                                               gCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGt                                                                                                                                     <  1:261397/251‑1 (MQ=255)                                                                                               gCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGt                                                                                                                                     <  2:54219/251‑1 (MQ=255)                                                                                                     ggCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAAccc                                                                                                                               <  1:84341/251‑1 (MQ=255)                                                                                                          cACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTgg                                                                                                                          <  1:375495/251‑1 (MQ=255)                                                                                                            ccAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCGTTGGTc                                                                                                                        <  1:253072/251‑1 (MQ=255)                                                                                                             cAGGGAGTCATCGGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCATTTGGTCa                                                                                                                       <  2:65766/251‑1 (MQ=255)                                                                                                                  aGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTAt                                                                                                                  <  1:284458/251‑1 (MQ=255)                                                                                                                     cATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATtct                                                                                                               >  1:130203/1‑251 (MQ=255)                                                                                                                                                cgcgAGGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCAGCTTCGACAATc                                                                                    <  2:129265/251‑1 (MQ=255)                                                                                                                                                          aaTGGTGCCAAGCGCGGTAAGCAACAGAATCGGCAGCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCa                                                                          <  2:35546/251‑1 (MQ=255)                                                                                                                                                                   aaGCGCGGTAAGCAACAGACTCGGCATCACTGTATTGGCGGAGCGTAACATGCGCACGATATACCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCAACGGTCATCGCCAGATGGTAGCCATTAAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGATAACCCTTTGGTCAAGTATTCTACGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTtcctc                                                                 >  2:268184/1‑251 (MQ=255)                                                                                                                                                                     gcgcGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCg                                                               <  2:76788/251‑1 (MQ=255)                                                                                                                                                                      cgcgGTAAGCAAAAGCATCCGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCCGCCGATAACGTCCGGAAGCATAATATCGAGGCTTATCAGATCATAATCACCGGTCATCGCCAGAAGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGCTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGc                                                              >  2:229601/1‑251 (MQ=255)                                                                                                                                                                            aaGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTg                                                        <  2:130203/251‑1 (MQ=255)                                                                                                                                                                              gaaaCAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATGATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCAAGTTGcc                                                      <  2:111946/249‑1 (MQ=255)                                                                                                                                                                                                        gCGGAGCGTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCa                            <  1:93682/251‑1 (MQ=255)                                                                                                                                                                                                               gTAACATGCGCACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGc                      >  1:99248/1‑250 (MQ=255)                                                                                                                                                                                                                       cgcACGATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTa             <  2:187435/251‑1 (MQ=255)                                                                                                                                                                                                                            gATATCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGc        >  1:224614/1‑251 (MQ=255)                                                                                                                                                                                                                              tatCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGttt                                                                                                        <  2:2058/153‑1 (MQ=255)                                                                                                                                                                                                                              tatCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGttt                                                                                                        >  1:2058/1‑153 (MQ=255)                                                                                                                                                                                                                               atCCCAGCCGCTCACGTCCTGCAGCACAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCaaa     <  2:323856/251‑1 (MQ=255)                                                                                                                                                                                                                               atCCCAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCaa      <  2:48175/250‑1 (MQ=255)                                                                                                                                                                                                                                   cAGCCGCTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATg  >  1:67978/1‑250 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCACGATATCCCAGCCGTTCACGTCCGGCAGCATAATATCGAGGATTATCAGATCATAATCACCGGTCATCGCCAGATGGTAGCCATTCAGCCCGTTGTCGGCCAAATCGACCACAAAACCGGCTTCGGTTAACCCTTTGGTCAAGTATTCTCCGGTTTTCTTTTCATCTTCGACAATCAACAGTTTCATATTTCCTCCGCATGTTGCCCGGGCAATTCTAGAGTAGCGGGATCAGATGGCAATCGCTTATTGGCAAAATG  >  NC_000913/595039‑595514

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