Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	850,995	C→G	8.0%	intergenic (‑575/‑230)	pcs ← / → hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	850,995	0	C	G	8.0%	184.9 / 6.6	63	intergenic (‑575/‑230)	pcs/hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCA  >  NC_002947/850854‑851144                                                                                                                                                |                                                                                                                                                      gCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                                               <  2:566073/150‑1 (MQ=255)    tttGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                            <  1:1494030/150‑2 (MQ=255)      tGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATGGATGCCCGTCaagc                                                                                                                                          <  1:1620322/150‑4 (MQ=255)      tGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCt                                                                                                                                          <  1:1475414/150‑1 (MQ=255)        ctctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATGGGGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTcc                                                                                                                                        <  1:1185607/150‑1 (MQ=255)         tctTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTccc                                                                                                                                       >  1:555324/1‑150 (MQ=255)           ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTg                                                                                                                                     >  2:1913240/1‑150 (MQ=255)                   gCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                             >  1:1918308/1‑150 (MQ=255)                    cTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGTA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAgcc                                                                                                                             <  1:1364782/150‑3 (MQ=11)                     ttCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCt                                                                                                                           <  2:1701067/150‑1 (MQ=255)                        tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TGGGG‑CCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtcc                                                                                                                                    <  2:754008/138‑3 (MQ=14)                        tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TGGGG‑CCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtcc                                                                                                                                    >  1:754008/1‑136 (MQ=14)                           gcgcgcGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                     >  1:1683841/1‑150 (MQ=255)                            cgcgcgGTTGCTGACGCGCTGTATATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCg                                                                                                                    <  2:1049449/150‑1 (MQ=21)                            cgcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGTTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCg                                                                                                                    <  2:1610688/150‑1 (MQ=255)                            cgcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCg                                                                                                                    <  2:313936/150‑1 (MQ=255)                                   tggATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCTGG‑CACAGCTCTTGACGTGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctcct                                                                                                             <  2:565688/148‑1 (MQ=255)                                     gATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGa                                                                                                           >  2:171067/1‑150 (MQ=255)                                      aTCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGaa                                                                                                          <  1:201056/150‑1 (MQ=255)                                       tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                         >  2:156323/1‑150 (MQ=255)                                       tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCAAAGTACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                         >  1:1453583/1‑150 (MQ=25)                                         aGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCACGGCCTGCGTTCGGCCTCCTGAAGTc                                                                                                       >  1:726663/1‑150 (MQ=255)                                               tGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCACGGCCTGAGTTGTTCCTCCTGAAGTCGCGaag                                                                                                 >  1:575229/1‑150 (MQ=21)                                                 cAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagat                                                                                               >  1:1861680/1‑150 (MQ=255)                                                  aGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                  >  2:1802683/1‑146 (MQ=255)                                                  aGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                  <  1:1802683/146‑1 (MQ=255)                                                      tGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaaga                                                                                          <  2:1013701/150‑1 (MQ=255)                                                       gCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagat                                                                                         <  2:332826/150‑1 (MQ=255)                                                            ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGa                                                                                    <  2:512930/150‑1 (MQ=255)                                                                ggaggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCaa                                                                                >  1:245993/1‑150 (MQ=255)                                                                  aggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCaaaa                                                                              <  2:1135674/150‑1 (MQ=255)                                                                    gCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGc                                                                            <  2:1918308/150‑1 (MQ=255)                                                                    gCCGAGCGTAAGGCTTGCTCCCGGGAGGACCGCGATGCAGAACTATTCA‑TAGCG‑CACATTGCTGGACGTGCATTAGTCTTCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCTTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGc                                                                            <  1:889893/150‑1 (MQ=12)                                                                                 cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                               >  2:564251/1‑150 (MQ=255)                                                                                    gCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGAATGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGc                                                            >  1:1294581/1‑150 (MQ=255)                                                                                        ccGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                          >  2:1037745/1‑148 (MQ=255)                                                                                        ccGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                          <  1:1037745/148‑1 (MQ=255)                                                                                        ccGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                        <  1:1892326/150‑1 (MQ=255)                                                                                        ccGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTAACCTCCCTAATCAAGGCGTGGGTGCTGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                        >  1:704303/1‑150 (MQ=21)                                                                                          gggAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAcc                                                      <  2:1683841/150‑1 (MQ=255)                                                                                               gaCCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaa                                                 >  1:477227/1‑150 (MQ=255)                                                                                                  cgcgCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAaca                                              >  1:1902958/1‑150 (MQ=255)                                                                                                   gcgcTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                             <  1:1492792/150‑1 (MQ=255)                                                                                                             ccccATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                        <  1:517529/97‑1 (MQ=255)                                                                                                             ccccATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                        >  2:517529/1‑97 (MQ=255)                                                                                                              cccATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAg                                  >  1:1014278/1‑150 (MQ=255)                                                                                                                    c‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcg                                                                                                      >  1:1661916/1‑76 (MQ=255)                                                                                                                    c‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcg                                                                                                      <  2:1661916/76‑1 (MQ=255)                                                                                                                      tCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGcc                                                                           >  2:1761103/1‑102 (MQ=255)                                                                                                                      tCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGcc                                                                           <  1:1761103/102‑1 (MQ=255)                                                                                                                            cacaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                      <  1:171067/150‑1 (MQ=255)                                                                                                                             acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGAt                                                                                   >  1:225658/1‑88 (MQ=255)                                                                                                                             acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGAt                                                                                   <  2:225658/88‑1 (MQ=255)                                                                                                                               aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                   <  2:939787/150‑1 (MQ=255)                                                                                                                                 cgcTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATa                 <  2:575229/150‑1 (MQ=255)                                                                                                                                    gtGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGTCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAACCACCGCTAGGGCAGGCTACGAAAAGGCGATAGcc              >  2:1618928/2‑150 (MQ=255)                                                                                                                                     tGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg             >  1:1794158/1‑150 (MQ=255)                                                                                                                                     tGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg             >  1:1043358/1‑150 (MQ=255)                                                                                                                                     tGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg             >  1:1441113/1‑150 (MQ=255)                                                                                                                                      gACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGc            >  1:69237/1‑150 (MQ=255)                                                                                                                                           gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGATACGAAAAGGCGATAGCCGCTgcgc       >  1:392162/1‑150 (MQ=255)                                                                                                                                            cATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCt      >  1:1836487/1‑150 (MQ=255)                                                                                                                                              tCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAt    <  2:1294581/150‑1 (MQ=255)                                                                                                                                                cATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCa  <  2:1861680/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      GCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCC‑TCGGG‑CACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCA  >  NC_002947/850854‑851144

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