Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,474,160	G→A	7.4%	intergenic (‑65/+115)	PP_2167 ← / ← tal	PhoD‑like phosphatase/transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,474,160	0	G	A	7.4%	185.5 / 3.0	67	intergenic (‑65/+115)	PP_2167/tal	PhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (2/3);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474014‑2474308                                                                                                                                                   |                                                                                                                                                       gccagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCGTAGATGGATCAGGGTAggc                                                                                                                                                      >  4:162633/1‑148 (MQ=255)  ccagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAg                                                                                                                                                     <  2:168588/149‑1 (MQ=255)    agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                                   <  1:181889/149‑1 (MQ=37)     gacggccatacgCTGGGGTTCCAGGCGGGGGAGTACCGGGCCGGCGAGTACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGGAAGCCCGCTCCTACACTGATCGATGTAAGGGCgg                                                                                                                                                  <  3:215329/149‑1 (MQ=25)       cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGTGCCTGCGAGAACGAGGGGCAACGGGTTAGGTGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCggg                                                                                                                                                 <  7:171554/148‑1 (MQ=255)       cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGTTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGc                                                                                                                                                <  1:281669/149‑1 (MQ=255)       cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAgggcgggc                                                                                                                                                <  7:103282/149‑1 (MQ=37)         gccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCtt                                                                                                                                              <  2:214096/149‑1 (MQ=255)          ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑g                                                                                                                                            <  1:104736/149‑1 (MQ=255)          ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCCGGGCAGCCCCGCCAGACATCGCTGAGTGTTGGGGCgggggtg                                                                                                                                             >  8:283590/1‑145 (MQ=255)           cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑Gc                                                                                                                                           >  4:46993/1‑149 (MQ=255)             gacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑Gccc                                                                                                                                         <  2:52787/149‑1 (MQ=255)               cgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCcgc                                                                                                                                       <  3:46993/149‑1 (MQ=255)                gCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGCCAAACCCGCCCGTTCATAGCTAAGTTTAGGAGCCGGCGT‑Gccccc                                                                                                                                       >  1:300949/1‑146 (MQ=255)                      gTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAgg                                                                                                                                <  6:139399/149‑1 (MQ=255)                       ttCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCGTACCTGGGGGAGGGTTGGAGCCGGCTT‑GCCCGCGAACCGGc                                                                                                                               >  2:1398/1‑149 (MQ=255)                               ggcgCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTACGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTGCAGAGAGCTGTGTTGGGGCGGGGTT‑GGCCGCGGACCGGCCC‑GAGAAt                                                                                                                      >  7:317891/1‑149 (MQ=255)                                gcgcAGTACCGGGCCGGCGAGAACGAGGGCCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAAtt                                                                                                                     <  3:162633/149‑1 (MQ=255)                                  gcAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCGTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTc                                                                                                                    >  1:216490/1‑148 (MQ=255)                                    aGTACCGGGCCTGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCACTTACATCTATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAAt                                                                                                                 <  7:264070/149‑1 (MQ=255)                                     gTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGCGGGCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGAGCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATa                                                                                                                <  3:228235/149‑1 (MQ=255)                                     gTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATa                                                                                                                <  4:298843/149‑1 (MQ=255)                                        ccGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGGCAAGCCCGCCCGTTAAAGGCTGACTGTTGGCTCGGGCTT‑TGCAGCGAACACGCCC‑GAGAATTaaaaaac                                                                                                              >  4:229554/1‑144 (MQ=255)                                           ggCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAACCCCCCCCATACAAAAAAAAGTGTTAGATCGGGGTT‑TACAGCCGACACGCCC‑GAGAATTCAATaaaaaa                                                                                                           >  4:144768/1‑143 (MQ=255)                                            gCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCGT‑TCCCGCGCACAGGCCC‑GAGAATTCAATACACAAc                                                                                                          >  6:138795/1‑148 (MQ=255)                                                 cGAGAACGAGGGGCAACGGGTTAGGGTAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGc                                                                                                    <  3:272158/149‑1 (MQ=255)                                                    gaACGTGGGGCAACTGGTTAGGGGAGGTCATCGACAGCTCGGGGGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTg                                                                                                 <  5:16810/149‑1 (MQ=255)                                                    gaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTg                                                                                                 >  1:85510/1‑149 (MQ=255)                                                      aCGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCCTGAGAATTCAATACACAACTAAGGCGTgg                                                                                                <  1:205600/149‑1 (MQ=255)                                                         aGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                                            <  7:182391/149‑1 (MQ=255)                                                               aaCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAAc                                                                                      >  1:217274/1‑149 (MQ=255)                                                                    gTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCCCCCATACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GACAAATAAAAAAACAACTAAGGCGTGGCAAAAAAAACAGGCg                                                                                 >  4:301112/1‑149 (MQ=38)                                                                     ttAGGGGAGGTCATCGACAGCTCGGGTGGCTTTGGCGACCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCaa                                                                                                                  >  1:37368/1‑115 (MQ=255)                                                                     ttAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGACCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCaa                                                                                                                  <  2:37368/115‑1 (MQ=255)                                                                        ggggAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGaaa                                                                             <  3:301112/149‑1 (MQ=255)                                                                            aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGc                                                                                  <  7:175165/140‑1 (MQ=255)                                                                            aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGGCAAGCCCCCCCGAACATCGATCAGTGTAGGAGCCGGCTT‑GCCCGCGGACAGGCCC‑GAGGAATCAAAACACAAAAAACGCGGGGGCAACAAAACCGGCGGTAAATaa                                                                         >  6:270422/1‑148 (MQ=255)                                                                                 aTCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGaa                                                                     <  3:144768/148‑1 (MQ=255)                                                                                         cTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCATAAAAAGATAAGTGTAGGAGCGGGCTT‑GCCCCCGAACAGGGCC‑CAGAATTAAAAAAACAACAAATGCGTGGCAAACAAAACGGGCGGAAAATACTGAAAACTTcacg                                                            >  4:283180/1‑149 (MQ=255)                                                                                                ggCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCa                                                     <  1:1398/149‑1 (MQ=255)                                                                                                ggCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCa                                                     <  3:36658/149‑1 (MQ=255)                                                                                                            cTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                                         <  7:283590/149‑1 (MQ=255)                                                                                                            cTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCAGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                                         >  1:282456/1‑149 (MQ=255)                                                                                                               ttCGCGGGCAAGCCCGCGCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAGGGCGTGGCCAGCAAAACGGGCGGAAGATACTGAACACTTGCCGCCCGTCAGGCCAGCACGAGCTg                                      >  8:208332/1‑149 (MQ=38)                                                                                                               ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGc                                                                                               >  2:80532/1‑92 (MQ=255)                                                                                                               ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGc                                                                                               <  1:80532/92‑1 (MQ=255)                                                                                                                 cgcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGt                                     <  1:106544/148‑1 (MQ=255)                                                                                                                 cgcgGGCAACCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGGAAACACGAACTGt                                     >  2:32632/1‑148 (MQ=255)                                                                                                                  gcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTa                                    <  2:217274/148‑1 (MQ=255)                                                                                                                    gggCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGGCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTAt                                   >  1:68874/1‑148 (MQ=255)                                                                                                                      gCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAgg                               >  1:97518/1‑149 (MQ=255)                                                                                                                           cccgctcCTACACTGATCGATGTAAGGGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt                          >  5:118268/1‑149 (MQ=37)                                                                                                                             cgccaTTACTTCGTTCAGTGTTGGTCGGGGCTT‑GCCCGCTAAATGGCCC‑TTTATTTAATTACACAATTAAGGGGTGGCAAAAAAACCGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTc                        <  2:97518/149‑1 (MQ=255)                                                                                                                               cccTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAg                      <  2:306018/149‑1 (MQ=255)                                                                                                                                 cTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGGAAACACGAACTGTATCAGGCCTTTTCAGCt                    <  1:32632/149‑1 (MQ=255)                                                                                                                                    aCATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt                  >  8:314805/1‑148 (MQ=255)                                                                                                                                      aTCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACggg                                                                                   <  2:255476/81‑1 (MQ=255)                                                                                                                                      aTCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAAATCAATACACAACTAAGGCGTGGCAAACAAAACggg                                                                                   >  1:255476/1‑81 (MQ=255)                                                                                                                                       tCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc              >  7:217755/1‑149 (MQ=255)                                                                                                                                        cGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACt             >  8:325483/1‑149 (MQ=255)                                                                                                                                        cGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACt             >  1:144303/1‑149 (MQ=255)                                                                                                                                          aTCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCaaa                                               >  7:216416/1‑113 (MQ=255)                                                                                                                                          aTCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCaaa                                               <  8:216416/113‑1 (MQ=255)                                                                                                                                          aTCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc              <  1:99755/146‑1 (MQ=255)                                                                                                                                          aTCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc              >  2:99755/1‑146 (MQ=255)                                                                                                                                              gtgtAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCaa       <  7:314805/149‑1 (MQ=255)                                                                                                                                               tgtAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAg      <  3:229554/149‑1 (MQ=255)                                                                                                                                                 tAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCt    >  6:364254/1‑149 (MQ=255)                                                                                                                                                   ggAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCttt  <  2:85510/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTT‑GCCCGCGAACAGGCCC‑GAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474014‑2474308

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