Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,925	G→C	R77P (CGC→CCC)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,925	0	G	C	97.3%	108.7 / ‑6.7	38	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (21/15);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGAT  >  NC_000913/4636685‑4637156                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          aGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTcgc                                                                                                                                                                                                                                >  1:273707/1‑251 (MQ=255)                gTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGGTTATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGTCCGATTTGCTCTCTGGGGACCCAACGCACGGGCAACTGGCGCAGGCGGTTAATCAGCTCCAACCTCCCCCGGTTTC‑CGCCaaaaacgggggg                                                                                                                                                                                                                 >  2:76503/1‑241 (MQ=255)                   tCGGCATGTGTTTGTTTCTGTTCTATTTTTGACTGGCGGCGGTGGCAGCGTTGTTCGTACTGTCCATTTTTGTCAAAGAAGCGAAACCGGGTGCGCGAAGAGCAATGGAGGGGACGGGTATCGATACCGCTGCGTTGTTGGCGGAGCTGGAGAGGCCCGAGTGGTTCCCCGGTGTCGCAACGCATGGCGATCTGGCGCAGGCGTTTAATCAGCTACAACATCCTCCT‑CTGCGCGCCAATAGCGGTggtat                                                                                                                                                                                                               <  2:173961/250‑4 (MQ=255)                    cGGCTTGCGGTTGTTGCTGGGCTATTTTTTCCTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGTGCGTGCGAAGAGCAACGGAGGGGACTTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTTCCGATTTGCTCTATGGGGACCCAACGCATGGGCAACGGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCAtt                                                                                                                                                                                                              <  2:5292/250‑1 (MQ=255)                     ggCAGGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACTGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTaa                                                                                                                                                                                                            <  2:76657/251‑1 (MQ=255)                      gCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAAc                                                                                                                                                                                                           >  1:50674/1‑251 (MQ=255)                                                     ggtggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAGCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGTTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGa                                                                                                                                                                           <  2:34733/251‑1 (MQ=255)                                                      gtggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAATTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTTCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑GTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGAc                                                                                                                                                                          <  1:121351/251‑1 (MQ=255)                                                           ggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGAATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCGCCCG‑TTTCGCGCCAATATCGGAGGCATTAACAAAGTGCGCAATGAATATCCAT‑TCTATATGACCGATg                                                                                                                                                                     >  2:42367/1‑251 (MQ=255)                                                                   ccTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAGCATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACa                                                                                                                                                                                                          <  1:231430/207‑1 (MQ=255)                                                                   ccTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACa                                                                                                                                                                                                          >  2:231430/1‑207 (MQ=255)                                                                     tGGTTCGAACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAAAGCATGGGCAACCGGCGCCGGCGTTTAATCCGCTACAACATCCCCCG‑TTTCACGCCAATATCGGTGACATTAACAAAGGGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCa                                                                                                                                                            >  2:107030/1‑250 (MQ=255)                                                                                  aaaTTTTTGTAAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATAGAGACCGAAACGTTGCTGGCGGAGCTGGAGCGTCCAGATTTGCTCTTTGGGGACCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAAAATCCCCCG‑TTTGGCGCCAATATCGGTGGCATTAAAAAAGTGGGCAATGAATATC‑ATGTCTATATGACAGATGCGAAGGGCAAAGGATTGTTGGAt                                                                                                                                              >  2:108962/3‑251 (MQ=255)                                                                                       tttGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAGCTACCACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATA‑ATGTCTATATGACCGATGCGCCGGGCAAAGTATTGTTCGATTCGGc                                                                                                                                         >  2:165806/1‑251 (MQ=255)                                                                                          gTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCAGTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCaaa                                                                                                                                      <  2:3637/251‑1 (MQ=255)                                                                                           tCAAAGACGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGAATTGTTCGATTCGGCAAAt                                                                                                                                     >  1:176957/1‑251 (MQ=255)                                                                                           tCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACCAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCACAGTATTGTTCGATTCGGCAAAt                                                                                                                                     >  2:75705/1‑251 (MQ=255)                                                                                                   gTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACAACGCAACGTTGCTGGCGGAGCTCGCGCGTCCCGATTTGCTCTCTCGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCCCTATGACCGATGCGCACGGCAAAGTATTGATCGATTCGGCAAATAAAGCCGt                                                                                                                             >  1:28611/1‑251 (MQ=255)                                                                                                      aagccGGGCGTGCGAAGAGCAACGGAGGGGACTGTGATCGACACCGCAACGTTGCTGGCGGGGCTGGCGCGTCCCGATTTGCGCTCTGGGGACCCATCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGGGCAATGAATAAC‑ATGGCTTTATGACCGATGCTCAGGTCAAAGGAGTGTTCGAGTCGGCAAATAAAGTCGTTg                                                                                                                           <  2:94550/247‑1 (MQ=255)                                                                                                       aaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGCTTCGGCAAATAAAGCCGTTggg                                                                                                                         >  2:51752/1‑250 (MQ=255)                                                                                                                  ggAAGAGCAACGGAGGGGACGTTGATCGACTCAGCATCTGTGGTGGCGGAGCTGGCGCGTCCAGATTTGCTCGCTGGGGAGCCACCTCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑GTTCGCGACAAGATCGGTGGCATTAGCAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGTCAAAGTATTGTTCGATTCGGCAAATAAAGCCGGTGGACAGGATTATTc                                                                                                              <  2:142278/250‑1 (MQ=255)                                                                                                                   gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGCTTCGGCAAATAAAGCAGTTGGCCAGGATTATTcg                                                                                                             >  1:16111/1‑251 (MQ=255)                                                                                                                   gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcg                                                                                                             >  1:23111/1‑251 (MQ=255)                                                                                                                     agagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGCTTATTcgcg                                                                                                           >  1:41646/1‑251 (MQ=255)                                                                                                                      gagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCCAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCgcgc                                                                                                          >  2:162459/1‑251 (MQ=255)                                                                                                                             ggAGGGGCCGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGAGTCCCGATTTGCTCTCTGGGGACACAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGCATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATg                                                                                                   >  2:116755/1‑251 (MQ=255)                                                                                                                                         gatagccacCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTa                                                                                        <  1:150163/244‑1 (MQ=255)                                                                                                                                                              ggCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTg                                                                  >  1:58460/1‑251 (MQ=255)                                                                                                                                                                       ggCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAAGCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcagca                                                         <  1:116755/251‑1 (MQ=255)                                                                                                                                                                        gcgcGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCAc                                                        <  1:278788/251‑1 (MQ=255)                                                                                                                                                                               ccGATCTGCTCTCCGGGGACCCAACGCCCGGGCAACTCGCGCAGGCGTTTAATCAGCTACCCCACCCACCG‑TTTCGCCCCAATATCGGAGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTCTTGTTCGATTCGGCAACTAAAGCCGTTGGCCAGGCTTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCa                                                  >  2:131604/1‑250 (MQ=255)                                                                                                                                                                                                   ccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATC‑CCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCCGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGt                            >  2:250717/1‑251 (MQ=255)                                                                                                                                                                                                        gCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtg                        >  1:269936/1‑251 (MQ=255)                                                                                                                                                                                                             ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATc                                    >  1:89430/1‑234 (MQ=255)                                                                                                                                                                                                             ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATc                                    <  2:89430/234‑1 (MQ=255)                                                                                                                                                                                                                    ggCGCAGGCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGTGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTg            <  1:170654/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCGTTTAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCAc      <  1:162459/251‑1 (MQ=255)                                                                                                                                                                                                                              tttAATCAGCTACAACATCCCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGAt  >  1:16762/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          AGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCG‑TTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATC‑ATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGAT  >  NC_000913/4636685‑4637156

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