Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,214,325	T→C	F427L (TTC→CTC)	rpoD →	RNA polymerase sigma factor RpoD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,214,325	0	T	C	100.0%	230.9 / NA	72	F427L (TTC→CTC)	rpoD	RNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/34);  total (38/34)

CAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATG  >  NC_000913/3214079‑3214567                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cagcagATTGAAGAAGAAGCCGGCCTGACCATCGAGCAGGTTAAAGATAGCACCCGTCGTATGTCCATCGGTGAAGCGAAACCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGAGGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCtctc                                                                                                                                                                                                                                                <  1:22773/251‑1 (MQ=255)         tGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATACATTCGAATCCCGCCGTGGTTACAAGCTCTCCCCCTACg                                                                                                                                                                                                                                        >  2:146988/1‑251 (MQ=255)         tGAAGAAAAAACCGGCCTGACCACCGAGCAGGTAAAAGATATCACCCGGCGTATGCCCATCGTTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAAAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTAc                                                                                                                                                                                                                                         <  2:370977/250‑1 (MQ=255)           aagaagaaACCGGCATGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGc                                                                                                                                                                                                                                       >  2:395523/1‑250 (MQ=255)                 aaaCCGGCCTGACCATCGAGCAGGTTCAAGATATCAACCGTCGTATGTCCATCGGTGGAGCGAAAGACCGCCGTGCGCAGAAAGAGATGGTTGCAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGACATACACCCACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGCAGGCAAACTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGGGGTTCCAAGCTCTCCACCTCAGCAACCtgg                                                                                                                                                                                                                                >  2:378721/1‑251 (MQ=255)                  aaCCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCtggt                                                                                                                                                                                                                               <  1:178936/251‑1 (MQ=255)                   accggcctcccgaTCGACCAGGTTAAAGCTATACACAGTCTTATGTCACTCGGTGAAGCGAAAGTCCGCCGTGCGAAGAACGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTCTGAAATAAACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCGGAAGAAAGCGGTTGACAACTTCCACTATCGCCGTGGTTACAAGCTCTCCACCTACGCAACCtggtg                                                                                                                                                                                                                              >  2:253288/13‑251 (MQ=255)                    ccGGCCTGACCATCGAGCCGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGCACTTACGGCTGGTTATTTCTATCGCTAAGCAATACACCAACCGAGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATCCTGCCGGGGTTCCAAGCTCTCCCCCTACGCAACCtggtgg                                                                                                                                                                                                                             >  2:74640/1‑251 (MQ=255)                      ggCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGAt                                                                                                                                                                                                                           <  1:104032/251‑1 (MQ=255)                         cTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCg                                                                                                                                                                                                                        >  2:109047/1‑251 (MQ=255)                                   aGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCTTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATc                                                                                                                                                                                                              <  1:207775/251‑1 (MQ=255)                                           aaaGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGctc                                                                                                                                                                                                      >  1:290759/1‑251 (MQ=255)                                                   cAACCTTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATCCACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCgg                                                                                                                                                                                              <  1:395523/251‑1 (MQ=255)                                                   cAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGTAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCTCTCTATCGCgg                                                                                                                                                                                              >  2:236110/1‑251 (MQ=255)                                                             aTGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGgcgcg                                                                                                                                                                                    >  1:392877/1‑251 (MQ=255)                                                               gTCCTTCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTAAAGCGAACTGACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGgcgcgc                                                                                                                                                                                   <  2:100905/250‑1 (MQ=255)                                                                 ccaggggggAAGCAAAACCCCGCCGTGAGAAGAAAGTGATGGTTGAAGCGACCTTACTTCTGGTTTTTTCTAACGCTAAGAAATCCACCCACCGTGG‑TTTGGATTTCCTTGCCTTTCATCAGGAAGGCGATATGGTTCGTATAAAAGCGCTTGATAAAGTCGAAAACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCTCAcc                                                                                                                                                                                <  2:333368/243‑1 (MQ=255)                                                                 ccATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATGCCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCAc                                                                                                                                                                                 <  1:238061/250‑1 (MQ=255)                                                                       gTGAAGCAAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCTCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGAGGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGt                                                                                                                                                                          <  1:389609/251‑1 (MQ=255)                                                                           agcgaaagaaagccgTGCGAAGAAAGAGATGGTTGACGAGAACTTAAGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTc                                                                                                                                                                      >  1:345583/12‑251 (MQ=255)                                                                             cGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCg                                                                                                                                                                    <  1:284988/251‑1 (MQ=255)                                                                               aaaGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTAAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGt                                                                                                                                                                  >  1:416788/1‑251 (MQ=255)                                                                                     cgccgTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCAAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCatat                                                                                                                                                            <  2:131503/251‑1 (MQ=255)                                                                                     cgcagtgaaaccaacgcgcaggTTGACGCGAACTCATGTCTGCTTATTTCAATAGCCAAGAAATACACCACCCGTGG‑CTTCCAGTTCATTGACCTGTTTCAGGCAGGCAACATCGGTCTGATGGAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCAACTACGCAAACAGGTGGACTCGCCAGGCGATAACCAGCTCTATCGCGGAACGGGCGCGAGCCATCCGCTTTCCGGTGCatat                                                                                                                                                            >  2:47021/21‑251 (MQ=255)                                                                                          tGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTg                                                                                                                                                       >  1:345858/1‑251 (MQ=255)                                                                                             gaagaaAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACACGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTg                                                                                                                                                       >  2:1691/1‑248 (MQ=255)                                                                                               agaaagaGAGGGTTGAAGCGAACTGACTTCGGGTTATTGTTACCGCTAAGAAATACCCCACCCGTGG‑CTGGCATTTCCTTGCCCTGATTCAGGAAGGCACCTACGGTCTGATAAAGGCGTTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCTCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGAcc                                                                                                                                                  <  2:157652/251‑1 (MQ=255)                                                                                                  aagaGATGGTTGAAGCGAACTTACGTCGTGTGATTTCTATCGCTAAGAAATACACCAACC‑TGGTCTTGCTGTTCCTGGACCTGATTCAGGAGGGCAATACCGGTCTGATGAAAGCGGTTGACAAATTCGAATACCGCCGTGGTTACAAGCTCTCCTCCTACGCAACTTGGTGGATCCGTCAGGCGATCTCCCGCTCTATCGCGTATGAGGCTCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATc                                                                                                                                               <  1:419264/251‑1 (MQ=255)                                                                                                    gagaTGGTTGAAGCGACCTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATca                                                                                                                                              >  2:403400/1‑250 (MQ=255)                                                                                                       aTGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTAGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATcaaca                                                                                                                                           >  1:394557/1‑250 (MQ=255)                                                                                                              aaGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAAc                                                                                                                                   <  1:131881/251‑1 (MQ=255)                                                                                                               aGCGAACTTACGTCTGGTTATTTATATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCCCCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAAcc                                                                                                                                  >  1:109038/1‑251 (MQ=255)                                                                                                                 cGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACAACAACCGTGG‑CTTGCAGTTCCTTGACCTGAATCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGAGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCAACCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGt                                                                                                                                >  1:386913/1‑251 (MQ=255)                                                                                                                  gAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAACGCGGTTGATCCATTCGACTACCGCCGTGGTTCCAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGt                                                                                                                                >  1:55264/1‑250 (MQ=255)                                                                                                                  gAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGt                                                                                                                                >  1:91988/1‑250 (MQ=255)                                                                                                                       tACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑TTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTAttt                                                                                                                            <  1:403400/249‑1 (MQ=255)                                                                                                                                  aTTTCTATCGCTAAGAAATACTCCATCCGTGG‑CTTGCAGTTCTTTGACCGGATTCAGGAAGGCAACAGCGGTCTGATGAAAGCGGTTTATAAATTCGAACACCTCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGAGCCGGCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCTTATTCCGTTGCATACGATTTAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAtgct                                                                                                               <  1:431401/251‑1 (MQ=255)                                                                                                                                      cTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGCTCCGTCAGGCGATCACCCCCTCTATCGCGGACCAGGCGCGCACCATCCGTATTCAGGTGCATATGATTGAGACCGTCAACAAGCTCAACCGTATTTCTCGCCAGATgcagca                                                                                                            >  2:320074/1‑250 (MQ=255)                                                                                                                                      cTATCGCTAAGAAATACACAAACCGAGA‑CTTGCAGTTCATTGACCTGCTTCAGGAAGGCAACATCGGTCTGCTGAAAGCGGCTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCAACCTACGCAACCTGGTGGCTCCGTCAGGCGCTCACCCGCCCTATCGCGGATCGGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAAAAAGCTCAACCGTCTTTCTCGCCAGCTGCTGCa                                                                                                            >  2:289494/1‑250 (MQ=255)                                                                                                                                        aTCGCTAAGAAATCCACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAag                                                                                                          <  1:389646/250‑1 (MQ=255)                                                                                                                                        aTCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGGATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAAACGTATTTGTCGCCAGAtgcggcc                                                                                                            >  2:149161/1‑244 (MQ=255)                                                                                                                                        aTCGCTAAGAAATACAACAAACGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTAGCCAGATGCTGCa                                                                                                            >  1:361670/1‑248 (MQ=255)                                                                                                                                          cGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAag                                                                                                          >  2:2264/1‑248 (MQ=255)                                                                                                                                              aaGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATggg                                                                                                   <  1:2264/251‑1 (MQ=255)                                                                                                                                                   aTACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTg                                                                                              <  1:49343/251‑1 (MQ=255)                                                                                                                                                            ccGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGc                                                                                     <  1:1691/251‑1 (MQ=255)                                                                                                                                                            ccGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGc                                                                                     >  2:161003/1‑251 (MQ=255)                                                                                                                                                                     tGCAGTTCCTTGTCCTGATTCAGGAAGGCAACATCGGTCTGATGTACGCGGTTGATAAATTCGAATACCGCCGTGGGTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACACGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCATCAAGATCAACCGTATTTCTCGCCAGATGCTGAAAGAGATGGGCCGTGAACCGACGCCGgaaga                                                                              >  2:7998/1‑250 (MQ=255)                                                                                                                                                                         gTTCCTGGACCTGATTCAGGAAGGCAACATCGGGCTGAGGAAAGCGGTTGATAAATTGAAATCCCGCCGTGGTGGCAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCTGGGTTCCGGGGCATATGAGTGAGACGATCATCAAGCTCGACCGTAGTTCGCGCCAGAGGGTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTgg                                                                         <  1:253288/251‑1 (MQ=255)                                                                                                                                                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGACAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGCATTTCTCGCCAGATGCTGCAAGAGATGGGACGTGAACCGACGCCGGAAGAACTGGCTg                                                                      >  2:187510/1‑251 (MQ=255)                                                                                                                                                                                   cTGATTCAGGAAGGCAACATCGGTCTGATGCATGCGGTTGATGAATTCGACTACCGCCGTGGTTACAAGCTCTCAACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGGATTCCGGTGCATATGATTGAGACCATCAACAAGCCCAACCGTATTTCTCGCCAGCTGCTGCACTAGATGGGCCGCGAACCGACGCCGGCAGAACGGGGTGAACGTAt                                                               >  2:191217/1‑251 (MQ=255)                                                                                                                                                                                   cTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTAt                                                               >  1:394532/1‑251 (MQ=255)                                                                                                                                                                                      aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCt                                                            >  1:7996/1‑251 (MQ=255)                                                                                                                                                                                             aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCg                                                     <  1:207547/251‑1 (MQ=255)                                                                                                                                                                                                  aaCATCGGTCGGATGAAAGCGGTTGAAAAATTGGAATGCCGCCGTGGTTACAAGCTCTCCACCTGCGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGGGCAGAGGATTGAGACCATCATCAAGCTCAACCGTATTTCTCGCCAGATGCTTCGAGAGATGGGCGTTTAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGa                                                <  1:402093/251‑1 (MQ=255)                                                                                                                                                                                                         gTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACa                                              >  1:295509/1‑246 (MQ=255)                                                                                                                                                                                                          ccTGATGAAAGCGGTGAATAAATCCGAATACCGCCGTGGTTACAAGCTCTCCCCCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATcc                                        <  2:233976/250‑1 (MQ=255)                                                                                                                                                                                                              atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCaa                                    <  1:304171/251‑1 (MQ=255)                                                                                                                                                                                                                gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCGGAACGGATGCTGCTGCCGGAAGAAAAGATCCGCaa                                    >  2:239419/1‑249 (MQ=255)                                                                                                                                                                                                                aaaaGCGGTTGATAAATTCTAATCCCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCCCCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCAGATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAg                                  <  1:85470/250‑1 (MQ=255)                                                                                                                                                                                                                   cgcgttggATAAATTCGAATCCCTCCGGGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCGGATGCCGGAAGACAAGTTCCGCAAAGTGc                               <  2:416788/244‑1 (MQ=255)                                                                                                                                                                                                                        ttGATAAATTCGAATACCGCAGCGGTTACAAGCTCTCCACCTAAGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTACTGa                            >  1:155589/1‑249 (MQ=255)                                                                                                                                                                                                                           aTAAATTCGAATACCGCCGTGGTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGa                            >  1:175099/1‑246 (MQ=255)                                                                                                                                                                                                                               aTTCGAACAGCGCCGTGGTTACAAGCTCTCCTCCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGAGCAGGCGCGCGCCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCATCCGTAGTTCTCGCCAGAGGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTTATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCa                   <  2:353404/251‑1 (MQ=255)                                                                                                                                                                                                                                        ccgccgTGGTTACAAGCTCCCCACCTACGCAACCTGGTGGATCCGTCAGGCGAGCACCCGCTCCATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAATCTCAACCGTATGTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTTATGCCGGAAGACAAGATCCGCAAAGTGCTTAAGAGCGCCAAAGAGCCaa          <  2:379798/251‑1 (MQ=255)                                                                                                                                                                                                                                         cgccggagTTACAAGCTCTCCACCTACGCAACCTGTTGGCTCCGTCAGGCGATCCCCCGCTCGATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATGGCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGAGCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAAt         <  2:290759/244‑1 (MQ=255)                                                                                                                                                                                                                                          gccgTGGTTACAAGCTCTCCACCTACGCAACCTGGGGGATCCGTCAGGCGATCACCCGCTCTATCGCAGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCCTCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAAtc        >  2:352400/1‑251 (MQ=255)                                                                                                                                                                                                                                          gccgTGGTTACAAGCTCTCCACCTACGCAACATGGTGGACCCGTCCGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTACGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAAtc        >  2:214825/1‑251 (MQ=255)                                                                                                                                                                                                                                                gTTACAAGCTCTCCACCTACGCAACCTGTTGGATCCGTCAGGCGATCACCCGCTCTATCGCTGATCAGGCGCGCTCCAGCCGTAGTCCGGTGCATAGGATGGAGACCATCAACAAGCTCAACCGTAGTGCTCGCCAGATGCTGCAAGAGATGGGGCGTGAACCGACGCCGGAAGAACTGGCGGAACGTAGGCGGATGGCGGAAGATAAGATCCGCAAAGTGCGGAAGATTGCTAAAGAGcccttcgccatg  <  2:396785/251‑11 (MQ=255)                                                                                                                                                                                                                                                gTTACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATg  <  1:381156/251‑1 (MQ=255)                                                                                                                                                                                                                                                  tACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCgg                                                    <  2:125467/199‑1 (MQ=255)                                                                                                                                                                                                                                                  tACAAGCTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCgg                                                    >  1:125467/1‑199 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGG‑CTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATG  >  NC_000913/3214079‑3214567

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