Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,340,953	C→T	G226R (GGG→AGG)	melR ←	DNA‑binding transcriptional dual regulator MelR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,340,953	0	C	T	96.2%	156.1 / ‑2.0	53	G226R (GGG→AGG)	melR	DNA‑binding transcriptional dual regulator MelR
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (28/22);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

ATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTT  >  NC_000913/4340711‑4341193                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ttatCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGGTGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACACGCTGGAATATCCTCATTGCAt                                                                                                                                                                                                                                          >  2:165715/2‑251 (MQ=255)    tCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATAGAGAATACTTTTATCGGTATCGCTCCGTAACGCGCGAACGTGGTTGAGGCGCATCGCGGTAATGTACTGTTTCATCGTAAATTGCATGACCCGCTGGATTATCCTCATTGCATAGt                                                                                                                                                                                                                                       >  2:136418/1‑251 (MQ=255)                      cGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGt                                                                                                                                                                                                                     <  1:536517/251‑1 (MQ=255)                                         cTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGAt                                                                                                                                                                                                  >  1:266794/1‑251 (MQ=255)                                           aaGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATg                                                                                                                                                                                                 <  1:165715/250‑1 (MQ=255)                                            aGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGt                                                                                                                                                                                               >  1:234130/1‑251 (MQ=255)                                             gttgtcGGTATTTTTGCTGCAACATGCCGACAGTTTTTCCGAAAGGTCTGTAAAATCGACTACTTGAACGAAAGCCTGCCGTCAGTGCAATATCGAGAATACTTTTATCGGGATCGCTCAGGAACGCGCGAACGTGGTTGATGCGCAGCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTc                                                                                                                                                                                              <  2:659012/246‑1 (MQ=255)                                                 gCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACTCGCTGGAATATCCCCCTTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAgcg                                                                                                                                                                                          >  1:330271/1‑251 (MQ=255)                                                   ggTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGcc                                                                                                                                                                                        <  2:652422/251‑1 (MQ=255)                                                            gcggcgACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATa                                                                                                                                                                                <  1:458070/250‑1 (MQ=255)                                                                 gACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGTCGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGtttt                                                                                                                                                                           >  2:154558/1‑250 (MQ=255)                                                                 gACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGtttt                                                                                                                                                                           <  1:154558/250‑1 (MQ=255)                                                                      gcctaCATATTTGCCGAACGTTCTGTAAAAACGACTTCGTGAACGAAAGCCTGCCGTCAGGGCAAGATCGAGAATCCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCTTCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACAGCGTTGATGGTCAGCGCCTGAGCATAGTTGTCGGCaa                                                                                                                                                                     <  2:582452/247‑1 (MQ=255)                                                                            atatTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGc                                                                                                                                                               >  1:84466/1‑251 (MQ=255)                                                                                 tGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAg                                                                                                                                                           <  1:437411/250‑1 (MQ=255)                                                                                    cGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATc                                                                                                                                                       <  2:266794/251‑1 (MQ=255)                                                                                         gTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCt                                                                                                                                                  >  2:160951/1‑251 (MQ=255)                                                                                         gTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCt                                                                                                                                                  >  1:107830/1‑251 (MQ=255)                                                                                          tGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCACTAAAGCCCAGCATCTGGCt                                                                                                                                                  >  2:666596/1‑250 (MQ=255)                                                                                                aaaaactccttctTGATCGATAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTACCGGTATCGCTCAGTAACGCGCGAACGTTTTTGTTGCGCACCGCGGTAATGTACTGTTTCATCGTCAATTGCATGTCCCGCTGGAATATCCTCATTGCTTAGTTTGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGTTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATaa                                                                                                                                           <  2:21493/240‑1 (MQ=255)                                                                                                     cGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTg                                                                                                                                      >  1:338515/1‑251 (MQ=255)                                                                                                          actTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCAt                                                                                                                                 >  1:385819/1‑251 (MQ=255)                                                                                                                  gAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGcgcgcg                                                                                                                         >  1:670549/1‑249 (MQ=255)                                                                                                                  gAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGcgcg                                                                                                                           >  1:470152/1‑249 (MQ=255)                                                                                                                  gAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCgag                                                                                                                         <  1:234060/251‑1 (MQ=255)                                                                                                                   aaaGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCgaga                                                                                                                        <  1:664082/251‑1 (MQ=255)                                                                                                                       ccTGCCGTCAGGGCAATATCGAGAATACTTTTATCTGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGGTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCt                                                                                                                    <  2:347210/251‑1 (MQ=255)                                                                                                                            cGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGtttt                                                                                                               >  1:407796/1‑251 (MQ=255)                                                                                                                             gTCAGGGCAATATCGAGAATACTTTTATCGGAATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATGTCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGGTCATAGTTTTCGGCAATAAGGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGttttt                                                                                                              >  1:478642/1‑251 (MQ=255)                                                                                                                                    cAATATCGAGAATACTTTTATCGGTCTCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGGTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTtgtgtgtg                                                                                                       >  1:473209/1‑251 (MQ=255)                                                                                                                                     aaTATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATGAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTgc                                                                                                      >  1:207736/1‑251 (MQ=255)                                                                                                                                        atCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCa                                                                                                                                                                                                                <  2:360803/142‑1 (MQ=255)                                                                                                                                        atCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCa                                                                                                                                                                                                                >  1:360803/1‑142 (MQ=255)                                                                                                                                            agaATACTTTTATCGGGATCGCTCAGTAACGCGCGACCGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCACAGTTGGCGTTAAGTTTGCCGTGCTCAGCCACATCGTTGCTGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGAAGCTGTTTTTGTGTGTGCGCGAggg                                                                                               >  2:594193/1‑250 (MQ=255)                                                                                                                                                             aTCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCga                                                                                                                          <  1:241688/207‑1 (MQ=255)                                                                                                                                                                gCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGGAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTTGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGACGAGAGACGCTGGTTTTGTGTGTGCGCGAGGGTTTATTGTCCAGAATCGGtt                                                                           >  2:59075/1‑251 (MQ=255)                                                                                                                                                                    tctAACGCGCGAACGTGGATGATGCTCTACGAGGTACTGTACTGTTTCATCGTCAACTGCATGACCCACTGGAATATCCTACTTGCATAGTTAGCGTTACTGTTGACGTGCTCAGACACATCGTAGATGGTCAGCGCCTGCTCATAGTTTTCGGCAATAAAGCCCAGCATCTGCCTAACATAAAATTGCTCATGGCGCGCGACGCTGTTTTTTAGTGTGCGAGAGGTTTTATTGTTCAGAATCGGTTccc                                                                        >  2:35749/3‑250 (MQ=255)                                                                                                                                                                             cgAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCagag                                                                >  1:533333/1‑249 (MQ=255)                                                                                                                                                                              gAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAgg                                                             >  2:55436/1‑251 (MQ=255)                                                                                                                                                                                 cGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTa                                                          <  1:400779/251‑1 (MQ=255)                                                                                                                                                                                  gTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGATAACATAAAATTGCGCATGGCGCGATACGCTGTTTTTGGGTGTGAGCGAGGTTTTATTGACCAGAATCGGTTCCCAGACAGAGAGGCTa                                                          >  2:203409/1‑250 (MQ=255)                                                                                                                                                                                           gcgcATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTGTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTg                                                 <  2:640313/250‑1 (MQ=255)                                                                                                                                                                                                  gcgGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAg                                         <  1:55436/251‑1 (MQ=255)                                                                                                                                                                                                    ggTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATc                                                                               >  1:393421/1‑211 (MQ=255)                                                                                                                                                                                                    ggTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAAGTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATc                                                                               <  2:393421/211‑1 (MQ=255)                                                                                                                                                                                                        aTGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAAt                                   >  1:6273/1‑251 (MQ=255)                                                                                                                                                                                                         tGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAAtt                                  >  2:7769/1‑251 (MQ=255)                                                                                                                                                                                                               gTTTCATCGTCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTtcatc                             <  1:7769/250‑1 (MQ=255)                                                                                                                                                                                                               gTTTCATCGTCAAGTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGGTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGTCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGATGCTAAATCGCTTGAGCATCAGGCCAATTtcatc                             <  1:398393/250‑1 (MQ=255)                                                                                                                                                                                                                        tCAATTGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCATCAAGAGAGGCTAAATCGCTTGAGCATACGGCCAATTTCATCAATGGCGAGc                   >  2:370410/1‑251 (MQ=255)                                                                                                                                                                                                                             tGCATGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCg              >  1:418168/1‑251 (MQ=255)                                                                                                                                                                                                                                   gcccGCTGGAATAGCCGCATGGCATAGGTGGGGTGAAGTTGGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTGTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCAGCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTg        <  2:224126/250‑1 (MQ=255)                                                                                                                                                                                                                                          ggAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGtt  <  1:666596/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTT  >  NC_000913/4340711‑4341193

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