Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,658	T→C	D46D (GAT→GAC)	creA →	PF05981 family protein CreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,658	0	T	C	100.0%	182.6 / NA	57	D46D (GAT→GAC)	creA	PF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/22);  total (35/22)

TAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCG  >  NC_000913/4635418‑4635879                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              tAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGATCGTGTAATAGATAAAAATGGAAACAATATGAAATACAAGCATTTGAACCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAACATTGTTGTGGAAGCCTTTGACGATCCCGAtg                                                                                                                                                                                                                     >  2:119765/1‑251 (MQ=255)  aGGGCTTGCCGGAAATAAAAGTATTGAGCTTTTGTTCTTCATCAATATGTTAATTACCGTGACGAACTTATTGCTCGTGTAATAGATACAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATACTGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAACATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGCTCCCGAtgt                                                                                                                                                                                                                    >  2:372489/1‑251 (MQ=255)   gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATGTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGAtgtg                                                                                                                                                                                                                   <  2:548836/251‑1 (MQ=255)   gcgCTTGACGGAAATAAGAGTATTGAGATTTTTTTCTTAATCAATATGTTCTGTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACATGCCTTTGATCCTGTCTTTAAGACTGATAATGCTGGGGCCATTGGCTAATGAAGAAGAGATTGGTTAGGTCGACACCGTATTTCAAATGATCGGCCCGGATAACAAAATGTTTGTGAAAGCCTTTGACGATCCCGAtggg                                                                                                                                                                                                                   >  2:7405/3‑249 (MQ=255)           cGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGt                                                                                                                                                                                                           <  1:228971/251‑1 (MQ=255)                    aaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGtt                                                                                                                                                                                                   >  1:325479/1‑250 (MQ=255)                      gTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATAGGCCCGGATCACAAAATTGTTGTGGAAGCTTTTGACGATACCGATGGGCAAAATGTCACCTGTTAtg                                                                                                                                                                                                >  2:24098/1‑251 (MQ=255)                                      cTTAATAAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACACAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAAc                                                                                                                                                                                >  2:503419/1‑251 (MQ=255)                                        tAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGaaaa                                                                                                                                                                                 <  1:63016/248‑1 (MQ=255)                                              atatGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAACATGTCACCTGTTATGTGAGCCGGGCGAAAACAggtggt                                                                                                                                                                         >  1:579038/1‑250 (MQ=255)                                              atatGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTa                                                                                                                                                                        <  2:438670/251‑1 (MQ=255)                                                    tATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGTTAAGAATATTAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTTGGTCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGTCCCGGATCACAAAATGGTTGTGGAAGCCTTTGACGATCCCGATATGAAGAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTaaa                                                                                                                                                                   <  1:284485/250‑1 (MQ=255)                                                     aTTTACCGTGCCGAACTAATTGCTCGTGTAATAGATAACAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGCCGATCACTATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAgg                                                                                                                                                                 >  2:4952/1‑251 (MQ=255)                                                         aCCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAggggg                                                                                                                                                              <  1:338572/250‑1 (MQ=255)                                                                      aaTTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCTTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCg                                                                                                                                                 <  1:403288/250‑1 (MQ=255)                                                                        ttGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGACTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGGGGAATTAAAGGGGGATTGGATCTGGCGGa                                                                                                                                               >  2:363465/1‑250 (MQ=255)                                                                              gtgtAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCGCCTGTGATGTGAGGCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATAcc                                                                                                                                        <  2:419771/251‑1 (MQ=255)                                                                                     aGATAAAAAAGGTAACAATATGACATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGGTACCTCCGATg                                                                                                                                 >  2:329426/1‑251 (MQ=255)                                                                                        tAAAAATGGTAACAATATGAAATACAAGCATTTGCTCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGGCTGGCGGAAGATACCTCCGATGCgg                                                                                                                              >  2:257354/1‑251 (MQ=255)                                                                                         aaaaaTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGc                                                                                                                             >  1:622964/1‑251 (MQ=255)                                                                                           aaaTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGct                                                                                                                            >  1:540417/1‑249 (MQ=255)                                                                                            aaTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGAATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGctat                                                                                                                          >  2:133230/1‑248 (MQ=255)                                                                                                     aaTATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                  >  2:310764/1‑250 (MQ=255)                                                                                                       tatGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTAGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGc                                                                                                               >  1:522916/1‑251 (MQ=255)                                                                                                        atGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGc                                                                                                               >  1:197345/1‑250 (MQ=255)                                                                                                            aaTACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTc                                                                                                          >  1:419632/1‑251 (MQ=255)                                                                                                               aCAAGCATTTGATCCTGTCTTTAAACCTGATAATGCTGGGGCCATTGGATCATGCAGAAGAGATTGGTTCGGACGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTCTTTCTTGTCAGCAAGTCggg                                                                                                       >  1:151184/1‑251 (MQ=255)                                                                                                                cAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGc                                                                                                      <  2:622964/251‑1 (MQ=255)                                                                                                                           tCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACATGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGCAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTg                                                                                           >  2:226245/1‑251 (MQ=255)                                                                                                                                   ttAAGCCTGATAATGCTGTGGCCATTGGCTCATGCAGAATAGATTGGTGCGATCGTCACCGTATTGAAAATGATCGGCGCGGATCACAAAATTGTTGTGGAAGACTGTGACGATCCCGATGTGAAAAATGACACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAGTAAAGGGGGAGTGGGTCTGGCGGAAGATACCTCCGATGCGGCGATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATc                                                                                    <  2:230497/250‑1 (MQ=255)                                                                                                                                   ttAAGCCTGATAATGCTGGGGCCAGTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAAGTGTTGTGGAAGCCGTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGAt                                                                                     <  2:325479/249‑1 (MQ=255)                                                                                                                                    tAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGTTCGt                                                                                  >  2:198970/1‑251 (MQ=255)                                                                                                                                        ccTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTAtt                                                                               >  2:322586/1‑250 (MQ=255)                                                                                                                                         cTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTa                                                                              >  1:499041/1‑250 (MQ=255)                                                                                                                                                     gggCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGc                                                                 >  1:240213/1‑251 (MQ=255)                                                                                                                                                       gCCATAGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGGATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTACAAACGGCAAAGCTc                                                               >  2:592637/1‑251 (MQ=255)                                                                                                                                                          aTTGGCTCATGCAGAAGAGCTTGGTTAGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAACATTGTTGAGGAAGACTTTGACGATACCGATGTGAAAACTGTCAACTGTTATGTGAGCCGGGCGTACACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTACGATGCGGCTATTTCTTGTCAGCCAGTCGGACCGATTGAACTGTAGGATCGCATTCAAAACGGCAAAGCTCAgg                                                            >  1:277831/1‑251 (MQ=255)                                                                                                                                                              gCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCg                                                         >  1:301744/1‑250 (MQ=255)                                                                                                                                                                 cATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGgt                                                     <  1:463326/251‑1 (MQ=255)                                                                                                                                                                    gCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACTGCAAAGCTCAGGGCGAGgtag                                                   <  1:527042/250‑1 (MQ=255)                                                                                                                                                                        aagaGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCGGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTc                                              <  2:103694/251‑1 (MQ=255)                                                                                                                                                                            gaTTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGGCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCaaaa                                          <  1:615399/251‑1 (MQ=255)                                                                                                                                                                                 gTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGc                                     >  1:566822/1‑251 (MQ=255)                                                                                                                                                                                 gTTCGGTCAACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGc                                     <  2:449529/251‑1 (MQ=255)                                                                                                                                                                                         gACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCt                             >  1:546515/1‑251 (MQ=255)                                                                                                                                                                                           cacCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTgg                           <  2:499041/251‑1 (MQ=255)                                                                                                                                                                                               gTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCtt                       <  1:550220/251‑1 (MQ=255)                                                                                                                                                                                                      aaaTGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCg                >  1:136470/1‑251 (MQ=255)                                                                                                                                                                                                          gATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTa             <  2:197345/250‑1 (MQ=255)                                                                                                                                                                                                             cGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAgg         >  1:547477/1‑251 (MQ=255)                                                                                                                                                                                                               gCCCGGATCACGAAATTGTTGTGCGAGCATTTGACGATCCCGATGAGAAAAATGTAAACTGTTATGTGCGCCGGGCGAAAACAGGCGGTATTAAAGAGGGCTTAGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCCAGTCGGGCCGATTGACCTGTCGGATCGTCTTCAAACCGGCAAAGCTCAGGGCGCGGTAGAATTCAAATAACGCACGACCCTGGTCTTTAAGTCGTTACAGGTc       >  2:627172/1‑251 (MQ=255)                                                                                                                                                                                                                  cGGATCACAAAATTGTTGTGGAAGCCTTTGGCGAATCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAGTAAAGGGGGAGTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGTATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTg    <  2:544093/251‑1 (MQ=255)                                                                                                                                                                                                                  cGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTg    <  2:457123/251‑1 (MQ=255)                                                                                                                                                                                                                   ggATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGCAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTgc   >  2:185055/1‑251 (MQ=255)                                                                                                                                                                                                                    gATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTgcg  >  1:352171/1‑251 (MQ=255)                                                                                                                                                                                                                    gATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTAGTTACAGGTCGTgcg  >  1:528296/1‑251 (MQ=255)                                                                                                                                                                                                                    gATCACAAAATTGTTGGGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGGATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCGGGTCTTTAAGTCGTTACAGGTCGTgcg  <  2:252056/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              TAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCG  >  NC_000913/4635418‑4635879

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