Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,075,234	C→T	V1217I (GTA→ATA)	putA ←	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,075,234	0	C	T	91.4%	91.1 / ‑3.8	37	V1217I (GTA→ATA)	putA	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (28/4);  total (30/7)
Fisher's exact test for biased strand distribution p-value = 5.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGG  >  NC_000913/1074987‑1075473                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 atatACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTg                                                                                                                                                                                                                                               >  1:303616/1‑251 (MQ=255)      cAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATTCAGCGCGTCATCCGGCCACAGGATCTGGCTGc                                                                                                                                                                                                                                          >  2:40109/1‑251 (MQ=255)        gCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGccc                                                                                                                                                                                                                                        >  1:284144/1‑251 (MQ=255)                       tatTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACACTTGTGCCATCCCGCGCGGCACCTGCTTCAAACAATGCGCGAAGCTGATCCGACTCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCCGATGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCTCCAGCAc‑ggcg                                                                                                                                                                                                                        >  1:22905/1‑251 (MQ=255)                            cTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCACATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGACTACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGACCTCCGGCCACAGTAGCGGGCTGCCCACCGCCAGCAC‑Ggcggcg                                                                                                                                                                                                                     >  2:169188/1‑249 (MQ=255)                                        gCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGCTGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAgcgc                                                                                                                                                                                                       >  1:363925/1‑251 (MQ=255)                                        gCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGACTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAgcgc                                                                                                                                                                                                       >  1:245053/1‑251 (MQ=255)                                          aaaaCCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCAt                                                                                                                                                                                                     >  1:108398/1‑251 (MQ=255)                                                 tGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCAACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGTCGAGCTGAGTCAGCGAGTCCTGCtc                                                                                                                                                                                              >  1:66389/1‑251 (MQ=255)                                                     ccaaaaCAATTGGGCCCACCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCGCCGTGGAAGATCACCGCATCAAACGGTTGAGCGTTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGCCGATGCAGCGCGTCATCCGGCCACATTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCtcatca                                                                                                                                                                                          <  1:368750/248‑1 (MQ=255)                                                           cAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGc                                                                                                                                                                                     >  1:364655/1‑250 (MQ=255)                                                                gtgCCATCACTCGCGGCAACTGCTTCAACCAATGCGCGCAGCTGATTCGAAGCACCGTGGAAGATCACCGCATCACACGATTGAGCGGTTATATTTTCCGCTTTCGCCAGGTGAATACGTGCGCTGACTGACGATGGCAATGCCTACACTAACTGACGATGCAGCGCGTCATCAGGCCACAGTATCTGGCAGACAACCGCC‑GAACCGGCTGACAGCTGGGTCAGCGCATCATGCTCTTCATCGGccaaac                                                                                                                                                                                >  2:340758/1‑245 (MQ=255)                                                                     aTCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACacaac                                                                                                                                                                           >  1:269575/1‑250 (MQ=255)                                                                         cgcgcgGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACAAAACAcgcg                                                                                                                                                                      >  1:317586/1‑251 (MQ=255)                                                                                       ttCACACAATGCGCGAAGCTGATCAGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGCTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTAACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATCAACAACACGCGCTCACGCGgcagc                                                                                                                                                         >  1:80288/1‑250 (MQ=255)                                                                                        tCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGgcagca                                                                                                                                                        <  1:144720/250‑1 (MQ=255)                                                                                          acacaATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCg                                                                                                                                                     >  1:171827/1‑251 (MQ=255)                                                                                                          gcagcagAATCACCTTGGATGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCGTTCGCCAGTTGAATACGTTCGCGGACGGCCGATGGCAATGTCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCGGGCTGCCCACCGCCAGCAC‑GGCGGCTAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTc                                                                                                                                     <  1:232921/245‑1 (MQ=255)                                                                                                             tCCGAATCCCCGTTGAAGATCCCCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCAATGTCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGACAACGGTAACTGGCTGCCAACTGCCAGCCA‑GGCGGCGCGCCGAGTCATCCCATCCTGCTCCTCCTCGGCAATACGCACCACGCCATCACGCGgcagccgagacaagcggttgagtgcacc                                                                                                                                  >  2:197520/1‑228 (MQ=255)                                                                                                                          ggAAGACCCCCGCATCAAACGGTGGCGGGTTTATATTTCCCGCTTCCGCCAGTATAAGACGTTCGCGGACTGCCGATGTCAAGTCCTTCTCTAACGGCCGATGCAGCGCGTCACCCGGCCACAGGATCTGGCGGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCCCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCgg                                                                                                                      <  1:52398/250‑1 (MQ=255)                                                                                                                            aaGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCTCGgc                                                                                                                     >  1:20082/1‑249 (MQ=255)                                                                                                                                  aCCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGCTGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTACAGGTGTTGCGTTCACCCGTCGGCACCGGCAGCAATCg                                                                                                             >  1:70616/1‑251 (MQ=255)                                                                                                                                         cAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGtt                                                                                                      >  1:39236/1‑251 (MQ=255)                                                                                                                                         cAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGAGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGATGTGtt                                                                                                      >  2:32852/1‑251 (MQ=255)                                                                                                                                                  gAGCGGATATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCCCAGTATCTGGCTGCCCACCGCCAGCAC‑GGAGGCGAGCTGAGACAGCGCCTCATGCTCATAATCGGCACTACACAACACGCGCTCACGCGGCAGCAGCTGCCAGGTGTTGCTATCAACCGCCGGCCCCGCCAGCAATCGTTGTGTTCttgcctgc                                                                                             >  2:170860/1‑251 (MQ=255)                                                                                                                                                       gTTATAATTTCCGCTTTCGCCCGTTGAATAAGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATCCACAAAACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCCACCGTCGGCCCCGGCAGCACTCGCTGTGTACCTGCCTGCGCCac                                                                                        >  2:150125/1‑250 (MQ=255)                                                                                                                                                                                    aCGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGc                                                            >  1:234285/1‑250 (MQ=255)                                                                                                                                                                                           cTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTgg                                                    >  1:158450/1‑251 (MQ=255)                                                                                                                                                                                                 gCCGATGGCAATGCCTTCCCTAACTGACGATGCCGCGCTGCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCACCATATTGCGTACATAACGCCTGCAATTCTGGCCGAtt                                              <  1:319738/251‑1 (MQ=255)                                                                                                                                                                                                           aTGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCAt                                    >  1:137296/1‑251 (MQ=255)                                                                                                                                                                                                                       aCTGACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGcccc                                                                                                                        >  1:199605/1‑155 (MQ=255)                                                                                                                                                                                                                          gACGATGCAGCGCGTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGAAGCAGCGGCCCTGAGGAGCGTTCACCCTTCGACCCCCGCAGCAATCGTTGAGTTCCTGCCTGCGACAACTCGCCATCTTGCGTACATAACGACAGCCATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCCttt                     >  2:271648/1‑251 (MQ=255)                                                                                                                                                                                                                          gACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCAttt                     <  1:250248/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGCGCGTCATCCGGCAACAGTATATGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGCGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCCGGAGTTGCGGTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTAAATAACGCCTGAAATTCTTGACGATTTGCTGCCCATTCCCGGAGGACATTTAACGGCTg             >  2:216057/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTCATCCGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTc         >  1:319204/1‑250 (MQ=255)                                                                                                                                                                                                                                           cccGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑TGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACATCACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGc    <  1:302217/250‑1 (MQ=255)                                                                                                                                                                                                                                             cGGCCACAGTATCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCgg  >  1:132681/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 ATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCTGATCCGAATCACCGTGGAAGATCACCGCATCAAACGGTTGAGCGGTTATATTTTCCGCTTTCGCCAGTTGAATACGTTCGCTGACTGCCGATGGCAATGCCTTCACTAACTGACGATGCAGCGCGTCATCCGGCCACAGTACCTGGCTGCCCACCGCCAGCAC‑GGCGGCGAGCTGAGTCAGCGCATCCTGCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGG  >  NC_000913/1074987‑1075473

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