Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,250	0	A	G	32.7%	112.7 / 19.4	55	E161G (GAG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base A (1/36);  new base G (18/0);  total (19/36)
Fisher's exact test for biased strand distribution p-value = 1.32e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGC  >  NC_000913/4622015‑4622473                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    tGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑Cggggagg                                                                                                                                                                                                                                   <  2:894690/238‑1 (MQ=255)   gCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑Cggggagggga                                                                                                                                                                                                                                >  1:1058725/1‑239 (MQ=255)           gAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTAt                                                                                                                                                                                                                          <  2:1127391/236‑1 (MQ=255)            aaGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAg                                                                                                                                                                                                                       <  2:890784/238‑1 (MQ=255)             aGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCCGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGTCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGT‑GGTGGGGGGGGGGGAG‑TTATCacc                                                                                                                                                                                                                      >  2:1413999/1‑237 (MQ=255)                tGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATg                                                                                                                                                                                                                  <  1:938109/239‑1 (MQ=255)                  tCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑TGT‑CGGGGGGGGGGA‑TTATCAGCGAAggg                                                                                                                                                                                                                >  1:150438/1‑238 (MQ=255)                       tAACAGGCGCAGATGGCGACGAACAGCGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCCACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGT‑GGG‑GGTGGGGGGGGG‑TGATCAGGGGAGGATTatt                                                                                                                                                                                                           >  2:976941/1‑237 (MQ=255)                              cgcAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACt                                                                                                                                                                                                           <  2:1058725/232‑1 (MQ=255)                                                  gggttGATGGTATGCCCTGGTATTGGACCGGTCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCctct                                                                                                                                                                                 <  1:546913/235‑1 (MQ=255)                                                  ggTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGGAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTg                                                                                                                                                                                <  1:423032/239‑1 (MQ=255)                                                           gTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGGT‑GGT‑GGGGGGGGGGGG‑ATTTCAACGGATGATTTCTTCACAACAACCCTCGGGAAATCctcctggtcacaaa                                                                                                                                                                       >  2:603097/1‑228 (MQ=255)                                                            tttGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaaaaa                                                                                                                                                                       <  1:1450957/238‑1 (MQ=255)                                                             ttGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGGGGGGGA‑TTATCAGCGATGGATTACTGCACAAAACCCCGCGGGAAATCCTCTGGATCAAAAAc                                                                                                                                                                      >  1:1601507/1‑238 (MQ=255)                                                             ttGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTGTGTT‑GGG‑GGGGGGGGGGGA‑TTTTGAGCGCTTGGATACTATAGAAAACCCTCCGGGAAATTCTCTCGCGCaaaaaac                                                                                                                                                                     >  2:637756/1‑237 (MQ=255)                                                                 ccTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAc                                                                                                                                                                  <  2:144995/238‑1 (MQ=255)                                                                 ccTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAc                                                                                                                                                                  <  2:486315/238‑1 (MQ=255)                                                                  cTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGTCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAccc                                                                                                                                                                <  1:1413999/239‑1 (MQ=255)                                                                                           cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑TGG‑CGGGGGGGGGAAGTTAAAAAGGA‑GGAGTTATGGGCAAAAACCCGCGGGAAATTCTCTGGTGCAAAAACTCACCCATTACGAACAGACGGCAAGTGaaga                                                                                                                                       >  2:843573/1‑236 (MQ=255)                                                                                           cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGG‑CGGGGGGGGGGA‑TTATTAGCGATGGATTACTGCACAAAACCCCGCGGAAATTCCCCTTGGTCAAAAAATCACCCATTACGAAAAGCTGGCCAGGGATg                                                                                                                                        >  1:147615/1‑238 (MQ=255)                                                                                           cAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGG‑CGGGGGGGGGGA‑ATATCAAGGATGGAATATTGCACAACACCCCGCGGGAAATACTCCGGGTCAAAAAAATCCCCATTACGAACAGACGGCAAATGagac                                                                                                                                       >  2:906553/1‑236 (MQ=255)                                                                                               gCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAc                                                                                                                                                                  <  1:746882/208‑1 (MQ=255)                                                                                                         gcgcgcGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGGGGGGGA‑TTATCAGCCGTGGGGTTCTACACAACAACCCTCGGGAATTCCCCTCGGACAAAAAAACACCCACTaacaa                                                                                                                                                        >  2:746882/1‑204 (MQ=255)                                                                                                                 ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTggg                                                                                                                 <  1:603097/238‑1 (MQ=255)                                                                                                                                                    gatgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCACTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACcta                                                                            <  1:910473/239‑1 (MQ=255)                                                                                                                                                      tgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTAc                                                                               <  1:843573/234‑1 (MQ=255)                                                                                                                                                                   tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACt                                                                                                                                                                                                           <  2:1447980/99‑1 (MQ=255)                                                                                                                                                                   tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACt                                                                                                                                                                                                           >  1:1447980/1‑99 (MQ=255)                                                                                                                                                                                aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCctct                                                                                                                                                                                 <  1:139047/112‑1 (MQ=255)                                                                                                                                                                                aTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACaaa                                                                <  1:795562/222‑1 (MQ=255)                                                                                                                                                                                     ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑Cgggga                                                                                                                                                                                                                                     <  1:761916/56‑1 (MQ=255)                                                                                                                                                                                     ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAAttt                                                                                 <  1:1109901/201‑1 (MQ=255)                                                                                                                                                                                      tctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATtgt                                                                                                                    <  1:1055559/166‑1 (MQ=255)                                                                                                                                                                                           gCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCctct                                                                                                                                                                                 <  1:266510/101‑1 (MQ=255)                                                                                                                                                                                            cccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTa                                                                                                                                                                                                             <  1:19907/72‑1 (MQ=255)                                                                                                                                                                                            cccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTAc                                                                               <  1:930080/196‑1 (MQ=255)                                                                                                                                                                                             ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTggg                                                                                                                 <  1:893514/162‑1 (MQ=255)                                                                                                                                                                                             ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCaaaa                                                                                    <  1:349904/190‑1 (MQ=255)                                                                                                                                                                                             ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGT‑GGT‑GGGGGGGGGGGA‑GTATTCGCGATGGGTTACTCCACAACACCCCGCGGGGAATTCTCCGGGTCAAAAACTCCCCCCTTACCAAAAGCTGGCAAGTGATGCACTGGCAGGGGGAAATTTGGGGTTGTTCT‑CTGGCGGGGGGCGACCCAAAATTACTCTA‑TATAGCACAAACGCCGCCAGGCTATAAACAAGCGTTggga                                   >  2:1500328/1‑237 (MQ=255)                                                                                                                                                                                               tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGT‑GGGGGGGGGGAA‑tat                                                                                                                                                                                                                           >  2:761916/1‑54 (MQ=255)                                                                                                                                                                                               tGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGG‑CGGGGGGGGGGA‑TTATCAGCGGTGGGGTTCTGTCCAACACCCTTCGGGAATTTCTCTGGGTCAAAAAATCAACCCTTACGAACAACTGGCGAACGATGCACCGGCA‑GGAGAAAATGTGGGGTCTTTTGCGGGGGGGGGGG‑CGCCCAAAATTActcacttat                                                                       >  2:1109901/1‑193 (MQ=255)                                                                                                                                                                                                   catATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACAGCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTgg                                                                                                                                  <  2:294836/138‑1 (MQ=255)                                                                                                                                                                                                     tATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTg                                                                                                                                           <  1:700497/129‑1 (MQ=255)                                                                                                                                                                                                     tATTTAACGGTGAATATACTGGCGGTTGGGT‑GGG‑CGGGGGGGGGGA‑ATATCAGCGGTGGATTTATGCCCAACAACCTGCGGGAAATCCTCCTGATCaaaaa                                                                                                                                                                       >  2:266510/1‑101 (MQ=255)                                                                                                                                                                                                      aTTTAACGGTGAATATACTGGCGGTTTGTT‑GGT‑GGGGGGGGGGGA‑TTTTTAGCGATGGGTTTCTGCacaaca                                                                                                                                                                                                   >  2:19907/1‑72 (MQ=255)                                                                                                                                                                                                       tttAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGGACAGCTGGCAAGTg                                                                                                                                           <  1:412157/127‑1 (MQ=255)                                                                                                                                                                                                       tttAACGGTGAATATACTGGCGGTTGGGT‑GGT‑CGGGGGGGGGGA‑TTTTACCCGACGGGTTACCGCCCAACACCCCGCGGGAATTACTCTGGGGCCAAAACTCACCCCCTTCCCAAAGCCGGCAAATGATGCACTGGCA‑GGGGGAAATTTGGGGTCTTCTgcg                                                                                                       >  2:893514/1‑162 (MQ=255)                                                                                                                                                                                                            cGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑Ggag                                                                                                                           <  1:1491652/137‑1 (MQ=255)                                                                                                                                                                                                             ggTGAATATACTGGCGGTTGGGT‑GGG‑CCGGGGGGGGGA‑TTTTTATGGGTGGGTTACTTCTCAAAAACCTCCGGAAACTCctct                                                                                                                                                                                 >  1:294836/1‑83 (MQ=255)                                                                                                                                                                                                               tGAATATACTGGCGGTTGGTT‑GGT‑CGGGGGGGGGGA‑TTATTAGCGAGGGGTTACTGGACAAAACCCCGCGGGAATTTCTCTTGGTCAAAAAATCACCCATTTCGAAAAACCGGCAAGTGATGCACTGGc                                                                                                                                 >  2:700497/1‑129 (MQ=255)                                                                                                                                                                                                               tGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAc                                                                                                                                                                  <  1:344717/96‑1 (MQ=255)                                                                                                                                                                                                                   tataCTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGaaa                                                                                                                        <  1:346124/133‑1 (MQ=255)                                                                                                                                                                                                                   tataCTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAAtt                              <  1:473470/221‑1 (MQ=255)                                                                                                                                                                                                                   tataCTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTAc                      <  1:224401/229‑1 (MQ=255)                                                                                                                                                                                                                        tGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCgg       <  1:204316/239‑1 (MQ=255)                                                                                                                                                                                                                           cGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACcgcg        <  1:906553/235‑1 (MQ=255)                                                                                                                                                                                                                              ttggtt‑ggT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGc  <  2:1601507/238‑1 (MQ=255)                                                                                                                                                                                                                              ttggtt‑ggT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGc  <  1:138525/238‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    TGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTT‑GGT‑CGGGGAGGGGAA‑TTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCA‑GGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGC‑AGCCAAAATTTACCTACTAT‑GCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGC  >  NC_000913/4622015‑4622473

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