Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,127,338	G→A	G178D (GGC→GAC)	yceM →	putative oxidoreductase YceM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,127,338	0	G	A	97.4%	100.4 / NA	38	G178D (GGC→GAC)	yceM	putative oxidoreductase YceM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/14);  minor base C (0/1);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCA  >  NC_000913/1127092‑1127578                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gggCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACaaa                                                                                                                                                                                                                                              <  1:155296/250‑1 (MQ=255)    cAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGcc                                                                                                                                                                                                                                           >  2:81699/1‑251 (MQ=255)      gAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCctc                                                                                                                                                                                                                                         >  1:91496/1‑251 (MQ=255)                   aTGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCCCCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGGTGTCGAGCGACACAGCCTCTCTGGAAGGCGGt                                                                                                                                                                                                                            >  2:90641/1‑251 (MQ=255)                      cTGAACGGCTGGTGGAACTGGCGGCGCGGAACAAACTGACGTTGATGGTCGGTTTTAACCGTAGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACg                                                                                                                                                                                                                         >  2:110771/1‑251 (MQ=255)                              cTGGTGGAACTGGCGGCGCGGAAAAACCTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGAc                                                                                                                                                                                                                 >  2:53485/1‑251 (MQ=255)                                     aaCTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACAGCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAacgac                                                                                                                                                                                                          <  2:267403/251‑1 (MQ=255)                                         ggcggcGTGGAAAAAACTGACGTTGATTGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCt                                                                                                                                                                                                       <  1:1151/250‑1 (MQ=255)                                                    aaaaaCTGACGGTGAGGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTATCAATAGCGTCGGGCCACACGAGCTTTATTTCACGTTGCTGGATGATTATCTGCATGGGGTGGATACCGCTCTGTTGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTGTCGATGCTGGCGAAGTGCTg                                                                                                                                                                                           <  2:14969/251‑1 (MQ=255)                                                           gccgataaTGGTCGGGTTTAACCGTCTATTCGCACCACTCTACGGTGGGTTCAAACCTCAACTCGCCCCCGCAGCCTCGCTCAGAATGGATGAACATCGTAGCAATAGCGTCAGGCCACACGATCTTTATTTCACGTTGCTGGGTGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACCAAGCCTCTATTGATGGCGTTACTCTCCTGACTAACGACACTGGCGTACTCCTGTTTGCCg                                                                                                                                                                                    >  2:153676/8‑251 (MQ=255)                                                                    ggTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCAt                                                                                                                                                                            <  2:243063/250‑1 (MQ=255)                                                                    ggTCGGTGTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGGAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATGCCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCAtt                                                                                                                                                                           <  2:91496/251‑1 (MQ=255)                                                                        ggTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTc                                                                                                                                                                       >  1:223263/1‑251 (MQ=255)                                                                        ggTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTc                                                                                                                                                                       >  1:78668/1‑251 (MQ=255)                                                                                  gtcgtTTCGCACCACTCTACGGTGAGTTAAACACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCt                                                                                                                                                             >  1:77250/1‑251 (MQ=255)                                                                                             ccacTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATcac                                                                                                                                                  >  2:4011/1‑251 (MQ=255)                                                                                                     cGGTGAGTTAAAAACGAAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCCATAGCGTCGGGCAACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGAAGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGtc                                                                                                                                          >  2:28171/1‑250 (MQ=255)                                                                                                                aaaCGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGAAATAGCGTCGGGCCACACGATATTTATTTCCAGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGATGTGGTTGACGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGACATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGcggt                                                                                                                               >  1:267643/1‑248 (MQ=255)                                                                                                                      aaCTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCAGTTTGCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGa                                                                                                                         <  1:230490/251‑1 (MQ=255)                                                                                                                          cGCCACCGCAGCCTCGCTAAGAATGGATAAACAGCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTc                                                                                                                     <  2:222561/251‑1 (MQ=255)                                                                                                                            ccACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCCTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTc                                                                                                                     >  1:18393/1‑249 (MQ=255)                                                                                                                             cACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGAGCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGc                                                                                                                  <  2:224700/251‑1 (MQ=255)                                                                                                                             cACCGCAGCCTCGATAAGAATGGCTAAACATCGTAGCAATCGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTACCGACGCTGGCGAAATGCTGTTTGCCGAGCACCCTTTTTCGGCTGGTCCTTTGCAGATCCCCCCCTGAATGCATCGCCGTGCCGGACGTCAGc                                                                                                                  >  2:40230/1‑251 (MQ=255)                                                                                                                                     ccTCGCTAAGAATGGATAAACATCGTAGCAAGAGCGTCGGGCCACACGAGCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGCCAAAGCCTCTCTTGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAAcc                                                                                                          <  2:242316/251‑1 (MQ=255)                                                                                                                                                aTGGATAAACATCGTAGCAAGCGCGTCGGGCCACGCGATCTTTAGTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTTCCGAGCACCATTTTTCGGCCGGTCCTTTTCAGAGCACCACCTGCAGGCAGCGCCGTGCCGGAATTCAGCGTGTACCCGTACAGGCCgc                                                                                               <  2:26900/251‑2 (MQ=255)                                                                                                                                                     tAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGtgactg                                                                                          >  1:189832/1‑251 (MQ=255)                                                                                                                                                       aaCATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGGGCCGGAAGTCAGCGTGCAACCGTGCAGGCCGtgactgtc                                                                                        >  2:8791/1‑249 (MQ=255)                                                                                                                                                         cATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGCTGATTATCTGCATGTGGTGGATACCGCGATGTGGTTATCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCAGGAAACCGTGCAGGCCGTGACTGACgg                                                                                      >  1:130331/1‑251 (MQ=255)                                                                                                                                                         cATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACgg                                                                                      <  1:110771/251‑1 (MQ=255)                                                                                                                                                                 cAATAGCATCGGGGCAAACGATCTTTATTTCACGTGGCTGGATGATTGTCTGCATGTGGTGGATACCGCGCTGTGGTTGACGAGCGACAAAGCCTATCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCGTTGCAGATCAGCAACTGTATGCATCGCCGTGCCGGAAGACAGCGGGAAACCGTGAAGGCAGTGACTGACGGTGCGCTcc                                                                              >  2:22969/1‑250 (MQ=255)                                                                                                                                                                   aTAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGAAGGTGCGCTCATc                                                                            >  1:252811/1‑251 (MQ=255)                                                                                                                                                                       cGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCCGGCCGt                                                                                               >  2:192096/1‑228 (MQ=255)                                                                                                                                                                       cGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGt                                                                                               <  1:192096/228‑1 (MQ=255)                                                                                                                                                                                             ttCACGTTGCTGGATGATTATCGTCATGTGGTGGATACCGCGCTGTGGTTATCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCAGTTTTCGGCTGGTCCTGTTCAGATCACCGCCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGc                                                   <  1:17511/250‑1 (MQ=255)                                                                                                                                                                                                         gatgatTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGAGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACTCTTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTCTGCATCGCCGTGCCGGAAGTCAGCATGAAACCGTGCAGGCCGTGACCGACGGCGCGCTAATCGACATTACGGATATGCGCGAATGGCGTGAGGGGCGCgg                                      >  2:169672/1‑251 (MQ=255)                                                                                                                                                                                                                 tCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCg                              >  1:188957/1‑251 (MQ=255)                                                                                                                                                                                                                  cTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCgg                             >  1:209672/1‑250 (MQ=255)                                                                                                                                                                                                                                              tGTCGAGCGACAAAGCCTCTCGGGATGGCGGTACGCTACTGACTATCGACGCTGGCGAAATGTTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCTGAAGTCAGCGTTAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCATGGCGTAGTGCAGAAACCGAATCCTGTTTGGcc  <  2:55244/250‑2 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 TGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCA  >  NC_000913/1127092‑1127578

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