Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	77,027	A→G	D91D (GAT→GAC)	sgrR ←	DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	77,027	0	A	G	97.8%	147.7 / ‑4.3	46	D91D (GAT→GAC)	sgrR	DNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/19);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGA  >  NC_000913/76784‑77271                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ccAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATAc                                                                                                                                                                                                                                               <  1:183770/251‑1 (MQ=255)    gggttgggcgagttcttattCCAGATCCCCCTTGTCTTCAGTTATGCGCGTTAGCGAACTGAATATTTGCCGGGACATATGGGTTTCGGAACGGCGCTATGCGTTGCCAGGTAGCATATTACGCAACGGACGATATTATAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGTTGCGGTTCAGATGAGAAACCAGCATTTGCCGCACAGTAGTTTTGTCGCTTACCAACTGCTCTAGTTGTTCGATAccat                                                                                                                                                                                                                                            <  2:286609/233‑4 (MQ=255)           gCGATGTCTGCTTCCAGTTTCCCATGTTCCTCATGTATGCGCGTTAGCGAACTGAAGATTTGCCGTGCAATATGGTTTTCGGAACGGCGCATTGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGTAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTc                                                                                                                                                                                                                                     <  2:394285/251‑1 (MQ=255)                    gCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCCCTAGTTGGTCGATACGATCCTTCTCCCGCAGGTc                                                                                                                                                                                                                            >  2:385699/1‑251 (MQ=255)                                 cATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGt                                                                                                                                                                                                               <  1:173953/251‑1 (MQ=255)                                                 gcgcGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGcc                                                                                                                                                                                               >  1:358411/1‑251 (MQ=255)                                                    cgTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCgg                                                                                                                                                                                            >  1:304362/1‑251 (MQ=255)                                                     gTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGCACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAAAAACCAGCATTTGCCGCACCGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCCGCTCCAGCAGGTCTTCCGCACGCTGTTGCTGAAGCGCCAGCCCGGt                                                                                                                                                                                           >  2:271149/1‑251 (MQ=255)                                                                   gATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCGGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCa                                                                                                                                                                             <  2:437594/251‑1 (MQ=255)                                                                      ttGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGa                                                                                                                                                                           <  2:75692/250‑1 (MQ=255)                                                                                   tGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGCAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTCTAGAGGACTGTCAGAAGCACGCGTTTAcc                                                                                                                                                             >  2:16275/1‑251 (MQ=255)                                                                                        ttCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGcc                                                                                                                                                        >  2:59370/1‑251 (MQ=255)                                                                                                     aaTGCGCTGCCAGGTAGCAGATTACGCCACGGACGATAGTATAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTc                                                                                                                                           >  1:337410/1‑251 (MQ=255)                                                                                                          gcTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACg                                                                                                                                      >  2:196715/1‑251 (MQ=255)                                                                                                            tGCCAGGTAGCAGATTCCGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGGCAGCCGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTc                                                                                                                                    >  2:341847/1‑251 (MQ=255)                                                                                                            tGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCACAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTc                                                                                                                                    >  2:196622/1‑251 (MQ=255)                                                                                                               cAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTcagc                                                                                                                                 >  1:293557/1‑251 (MQ=255)                                                                                                                     gCAGATTACGCAACGGAAGATAGGAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTGGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCcg                                                                                                                           <  2:94464/251‑1 (MQ=255)                                                                                                                            aCGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCt                                                                                                                    <  1:319276/251‑1 (MQ=255)                                                                                                                              gCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGc                                                                                                                  >  1:291/1‑251 (MQ=255)                                                                                                                                      cGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGtt                                                                                                          <  1:64232/251‑1 (MQ=255)                                                                                                                                                     cgcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGTCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCat                                                                                           <  1:68235/251‑1 (MQ=255)                                                                                                                                                      gcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTGTAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCata                                                                                          <  1:96128/251‑1 (MQ=255)                                                                                                                                                      gcAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCat                                                                                           >  1:364899/1‑250 (MQ=255)                                                                                                                                                        aGGATGTGCCGCCCCTGTCGGAAGCTGCGTCCCAGATGAGCAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCAATAGTTGGTCGAGACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGGATAGAGAAATGTCAGACGCGACCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCTGGCATGGGGTTGAGCAGGGTggcgatatg                                                                                        >  2:140402/1‑243 (MQ=255)                                                                                                                                                        aGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATg                                                                                        >  2:100199/1‑251 (MQ=255)                                                                                                                                                             gtgCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATgacgac                                                                                   <  1:53573/251‑1 (MQ=255)                                                                                                                                                                 cgccCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGAcgcgg                                                                               >  1:386939/1‑250 (MQ=255)                                                                                                                                                                     ccTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGagc                                                                             >  1:2480/1‑249 (MQ=255)                                                                                                                                                                           ggAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATCCGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAAt                                                                     <  1:1994/251‑1 (MQ=255)                                                                                                                                                                              aGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGCAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCTGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGATCAATacc                                                                  >  2:458973/1‑249 (MQ=255)                                                                                                                                                                                       ccAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGCATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCCCGTCAGCCAGCCGCGATCCTGCATGGCGTTGCGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATACCGCTTCCAGt                                                         >  2:12092/1‑251 (MQ=255)                                                                                                                                                                                       ccAGATGAGAAACCAGCATGTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGt                                                         <  1:331804/251‑1 (MQ=255)                                                                                                                                                                                         aGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTc                                                       <  1:100199/251‑1 (MQ=255)                                                                                                                                                                                         aGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTc                                                       >  1:103125/1‑251 (MQ=255)                                                                                                                                                                                               gaAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTg                                                   >  1:428851/1‑249 (MQ=255)                                                                                                                                                                                                      gCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTtgtg                                          <  1:341847/251‑1 (MQ=255)                                                                                                                                                                                                        aTTTGCCGCCCAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCGTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTc                                        <  2:428851/251‑1 (MQ=255)                                                                                                                                                                                                             ccGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCg                                   <  1:361037/251‑1 (MQ=255)                                                                                                                                                                                                             ccGCACAGTCGCTTTGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCg                                   >  1:373509/1‑251 (MQ=255)                                                                                                                                                                                                                        cttttccGCCAACCAACTGCTCTAGTTGGTCGTTCCGTTCCTGCCTCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTAccc                           <  2:386939/242‑1 (MQ=255)                                                                                                                                                                                                                          ttGTCGCCAACCAACTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCgc                       >  1:145219/1‑250 (MQ=255)                                                                                                                                                                                                                                    ccaacTGCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCcac            <  1:59370/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCACTAGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCg       <  1:391418/250‑1 (MQ=255)                                                                                                                                                                                                                                            actaGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGAt    >  1:133971/4‑251 (MQ=255)                                                                                                                                                                                                                                               aGTTGGTCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCACGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGa  >  1:407462/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGAGAAACCAGCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCCGACTTCCGCTTCCCACGTCAGCCAGCCGCGATCCTGCATGGTGTTGAGCAGGGTGCGCATATGACGACGCGAGCAGCTCAATAACGCTGCCAGTTCGTTGAGCGTTGTGTCCTGCGATTTACCCTCGCAGCATTGCCACAGGCGGATGA  >  NC_000913/76784‑77271

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