Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,584	T→C	G78G (GGT→GGC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,584	0	T	C	97.1%	103.8 / ‑6.6	35	G78G (GGT→GGC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (9/24);  total (9/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTG‑CTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC  >  NC_000913/71338‑71823                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 tGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGgctgggctgg                                                                                                                                                                                                                                               <  1:178322/250‑1 (MQ=255)                         gCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGACTGGCG‑TTTTAAAAAGCCGTTGCATCGCt                                                                                                                                                                                                                     <  2:65731/251‑1 (MQ=255)                            ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                                                                                  <  2:336469/251‑1 (MQ=255)                             gAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTc                                                                                                                                                                                                                 <  2:247733/251‑1 (MQ=255)                               acacTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCAt                                                                                                                                                                                                               <  1:278090/251‑1 (MQ=255)                                 acTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGAATGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCAttt                                                                                                                                                                                                             <  2:251761/251‑1 (MQ=255)                                    ttattaCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGGACGCCTGGGTGGGCGGGGTTGTCGGGTGTCTTG‑TGGGGGACTGGGTCTCCCTCTTGCTGGGCTGGCGGTTTTAAAAA‑CCGGTGCATCGCTGGGCAtttttt                                                                                                                                                                                                          >  1:95307/1‑248 (MQ=255)                                                  cccGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGc                                                                                                                                                                                            <  1:3238/250‑1 (MQ=255)                                                  aCCGTGTATGCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAg                                                                                                                                                                                             <  2:52009/250‑1 (MQ=255)                                                               tGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGAGCGCTGATTGGCAGCGGCGAGGTAAGGTTTTGGCACGCCTGGGTGGGAGGGGTTGTCGGGTG‑TCTGCTGGGGGACTGGGTCTCCTTTTTGCTGGGCGGGCGGTTTTAAAAA‑CCGTTGCATCGCTGGTCATTTCTGAAGAaaaaaaaaacattaatcgataa                                                                                                                                                                                >  1:50887/1‑236 (MQ=255)                                                                 aTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTTCTCGATAAAAc                                                                                                                                                                             <  2:29931/251‑1 (MQ=255)                                                                         ggggtTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGGTTGAATCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCGGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCAGGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTTCTCGATAAAACTGAACATg                                                                                                                                                                     <  2:189758/248‑1 (MQ=255)                                                                              gggtggCCGTTCTGGGGTCGCGGGTGCTGGTCGGTTTGATTCGACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCGCGTCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGTGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGTTCAGTGCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACAGGCGTTg                                                                                                                                                                <  2:42702/250‑1 (MQ=255)                                                                                  ggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTGTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATc                                                                                                                                                            <  1:4794/251‑1 (MQ=255)                                                                                    ccTTTCTGGAGTCGCTGGCGCGGGTCGGTTGGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCATGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTGCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCaa                                                                                                                                                          <  2:239647/251‑1 (MQ=255)                                                                                                cgctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGTATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACTTGCGTTGCATCAACACAGCATGtt                                                                                                                                               <  2:183737/250‑1 (MQ=255)                                                                                                     gcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAt                                                                                                                                         <  1:13744/251‑1 (MQ=255)                                                                                                      cgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAtt                                                                                                                                        <  1:322979/251‑1 (MQ=255)                                                                                                                tttGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTgggcggtc                                                                                                                              >  1:40386/1‑251 (MQ=255)                                                                                                                         ccccGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTttgttg                                                                                                                     <  2:50887/250‑1 (MQ=255)                                                                                                                         aCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGATGGCG‑TTTTAAAAAGCCGGTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTttg                                                                                                                        <  1:238961/248‑1 (MQ=255)                                                                                                                             ggTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAAc                                                                                                                                                                                                 <  1:291859/171‑1 (MQ=255)                                                                                                                                        gATGGCGGGGCTGGGGGCGCTGATTGGCAGCGGCGATTTAAGTTTCTTGCACGCCTGGCTGGCATGGATTGTCGGTTT‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTTGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCTATAAAACTGACCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGGTGGCCCGActcatct                                                                                                        <  2:231897/249‑7 (MQ=255)                                                                                                                                            gCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGt                                                                                                  <  1:70493/251‑1 (MQ=255)                                                                                                                                               gggCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGcc                                                                                               <  1:186286/251‑1 (MQ=255)                                                                                                                                                                       ggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCCGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCt                                                                       >  2:167354/1‑251 (MQ=255)                                                                                                                                                                              tAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTg                                                                >  1:91042/1‑251 (MQ=255)                                                                                                                                                                                       ggCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAt                                                         <  1:11003/249‑1 (MQ=255)                                                                                                                                                                                         cACGCCTGGCTGGCAGGGATGGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttcat                                                     <  2:330176/251‑4 (MQ=255)                                                                                                                                                                                               tggctggcAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACACGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGcc                                               >  1:324192/1‑251 (MQ=255)                                                                                                                                                                                                ggctggcAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCg                                              <  1:45257/251‑1 (MQ=255)                                                                                                                                                                                                     ggcagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAat                                           <  1:198057/227‑1 (MQ=255)                                                                                                                                                                                                          ggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATc                                     >  1:69227/19‑251 (MQ=255)                                                                                                                                                                                                                          gctGGGCGACTGGATCTCCTTCTGGCTGGTCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTTGTCATTTCTGAAGAAAAACAAAGCATTACTCGAGAAAACTGAACATGCGTTGCATTAACACAGCATGTTCACCATTCTGTTCGGTCGCTTTGTTGGCCCGACTCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGAGCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCtgtg                       <  1:332082/247‑1 (MQ=255)                                                                                                                                                                                                                                 gACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAAAGCTGGATCTGCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGccaccggt              >  2:1704/1‑249 (MQ=255)                                                                                                                                                                                                                                             ttctGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAACCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGACTGCTGTGGCCACCGTTTTACTTCCTGcc  >  1:345378/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 TGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTG‑CTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC  >  NC_000913/71338‑71823

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