Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,920,317	C→T	6.8%	R498Q (CGG→CAG)	acs ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,920,317	0	C	T	6.8%	251.0 / 4.4	85	R498Q (CGG→CAG)	acs	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base T (4/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTTGCT  >  NC_002947/3920169‑3920462                                                                                                                                                     |                                                                                                                                                   gcgcATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATTCTTGCCCTCGCCCACTCACGCAGCACTACCTCACTGACAGCATGCCCCGGGTGCAGCCg                                                                                                                                                   >  1:989393/1‑150 (MQ=255) gcgcATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCg                                                                                                                                                   >  1:559451/1‑150 (MQ=255)             tttCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCACACCCAGGCCAc                                                                                                                                       <  2:1615291/150‑1 (MQ=255)                 tGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGAtc                                                                                                                                   >  2:1163388/1‑150 (MQ=255)                     ccTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCt‑t                                                                                                                              >  2:188393/1‑150 (MQ=255)                       ttGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCa                                                                                                                            <  2:34249/150‑1 (MQ=255)                       ttGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCa                                                                                                                            >  1:1353384/1‑150 (MQ=255)                                  cgtcaTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTtgct                                                                                                                 <  2:814158/150‑1 (MQ=255)                                      aTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAg                                                                                                             <  1:1163388/150‑1 (MQ=255)                                      aTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAg                                                                                                             >  2:574428/1‑150 (MQ=255)                                       tCGGGAACTCGTCGACAAAACGGAAATAACGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAgg                                                                                                            <  1:1208581/150‑1 (MQ=255)                                       tCGGAAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGAGCGATCCGCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAgg                                                                                                            >  2:228675/1‑150 (MQ=255)                                          ggAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCa                                                                                                         >  2:1030123/1‑150 (MQ=255)                                           gAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGCGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACTACATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCAc                                                                                                        <  2:740112/150‑1 (MQ=255)                                             aCTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGc                                                                                                      <  1:139746/150‑1 (MQ=255)                                                 gtcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCa                                                                                                              <  2:1036732/138‑1 (MQ=255)                                                 gtcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCa                                                                                                              >  1:1036732/1‑138 (MQ=255)                                                  tcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTCGCTGCAg                                                                                                             <  2:395937/138‑1 (MQ=255)                                                  tcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTCGCTGCAg                                                                                                             >  1:395937/1‑138 (MQ=255)                                                     aCAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCAcc                                                                                              >  2:1264188/1‑150 (MQ=255)                                                     aCAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATCCCCCAAGCATAGCTGATTGCAGCCCTC‑ATCCCTGTCACCATACTTGCTGCAGGGCACGCCAATCAcc                                                                                              >  2:766297/1‑150 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTg                                                                                            >  2:1175288/1‑150 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCATTGCCCTCGCCCACTAACGCGTTTATACCTCACTGACTGCATGCTCGGAGTGCAGCCGCATACAGGACAAGATTTCT‑TACCCGTATTGAGTTCAGGGCACGCCAATCACTTg                                                                                            >  2:283623/1‑150 (MQ=255)                                                         aaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGca                                                                                          >  2:1328746/1‑150 (MQ=255)                                                            ggAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACAt                                                                                       >  1:171971/1‑150 (MQ=255)                                                             gAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATc                                                                                      <  1:754516/150‑1 (MQ=255)                                                              aaaTACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCg                                                                                     <  2:665959/150‑1 (MQ=255)                                                               aagaCCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCgg                                                                                    <  2:397300/147‑1 (MQ=255)                                                                aTACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGc                                                                                   <  2:1353384/150‑1 (MQ=255)                                                                  aCCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCa                                                                                 >  1:1555458/1‑150 (MQ=255)                                                                      gggCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCAccgc                                                                             >  1:1452379/1‑150 (MQ=255)                                                                      gggCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCAccgc                                                                             >  2:427558/1‑150 (MQ=255)                                                                       ggCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGAtctc                                                                                                                                 >  2:77140/1‑98 (MQ=255)                                                                       ggCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGAtctc                                                                                                                                 <  1:77140/98‑1 (MQ=255)                                                                         cACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgg                                                                          <  2:559451/150‑1 (MQ=255)                                                                         cACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgg                                                                          <  2:989393/150‑1 (MQ=255)                                                                         cACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgg                                                                          >  1:186224/1‑150 (MQ=255)                                                                                  gTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaa                                                                 <  1:1030123/150‑1 (MQ=255)                                                                                   tGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaag                                                                >  2:1510622/1‑150 (MQ=255)                                                                                     cgcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaa                                                              <  1:188393/150‑1 (MQ=255)                                                                                        aaTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaagaa                                                           >  2:739247/1‑150 (MQ=255)                                                                                        aaTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaagaa                                                           <  1:574428/150‑1 (MQ=255)                                                                                         aTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAAc                                                          >  1:465976/1‑150 (MQ=255)                                                                                            cTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCaa                                                                                                   <  2:1047811/106‑1 (MQ=255)                                                                                            cTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCaa                                                                                                   >  1:1047811/1‑106 (MQ=255)                                                                                            cTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActcc                                                       <  2:1555458/150‑1 (MQ=255)                                                                                              tGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActcctc                                                     <  1:739247/150‑1 (MQ=255)                                                                                               gCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCa                                                                                                                            <  1:1648377/78‑1 (MQ=255)                                                                                               gCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCa                                                                                                                            >  2:1648377/1‑78 (MQ=255)                                                                                               gCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACtcctcc                                                    <  1:766297/150‑1 (MQ=255)                                                                                               gCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑GCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACtcctcc                                                    <  1:763802/150‑1 (MQ=255)                                                                                                cccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCa                                                   >  2:1268977/1‑150 (MQ=255)                                                                                                cccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCGCCa                                                   <  1:228675/150‑1 (MQ=255)                                                                                                 ccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAg                                                  >  1:678515/1‑150 (MQ=255)                                                                                                   tCGCCCACTCACGCAACTCTACATCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCt                                                >  2:324483/1‑150 (MQ=255)                                                                                                     gCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgg                                                                          >  2:1347554/1‑122 (MQ=255)                                                                                                     gCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgg                                                                          <  1:1347554/122‑1 (MQ=255)                                                                                                        cACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGc                                                                                   >  1:1092041/1‑110 (MQ=255)                                                                                                        cACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGc                                                                                   <  2:1092041/110‑1 (MQ=255)                                                                                                        cACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTcgcgc                                           >  1:642974/1‑150 (MQ=255)                                                                                                        cACTCACGCAACTCTACCCCACTGACAGCATGCGCCGGGTGCAGCCGCACCCCGGGCCCGGCACCT‑TGCTGGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTcgcgc                                           <  1:1264188/150‑1 (MQ=25)                                                                                                           tCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCggg                                        >  1:560115/1‑150 (MQ=255)                                                                                                               gCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                             >  2:1598198/1‑93 (MQ=255)                                                                                                               gCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                             <  1:1598198/93‑1 (MQ=255)                                                                                                               gCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTaaa                                    >  2:1238695/1‑150 (MQ=255)                                                                                                                cAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAAt                                   <  2:642974/150‑1 (MQ=255)                                                                                                                cAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAAt                                   <  1:427558/150‑1 (MQ=255)                                                                                                                cAACTCTACCGCTCTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCATATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAAt                                   <  1:711946/150‑1 (MQ=255)                                                                                                                    tctACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtc                               >  1:1272683/1‑150 (MQ=255)                                                                                                                     ctACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtct                              <  1:1238695/150‑1 (MQ=255)                                                                                                                                aCAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTCGGTGAAGAAGGACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGa                   >  2:1119566/1‑150 (MQ=255)                                                                                                                                 cAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGaa                  >  1:1342429/1‑150 (MQ=255)                                                                                                                                    cATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                         >  2:979344/1‑124 (MQ=255)                                                                                                                                    cATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                         <  1:979344/124‑1 (MQ=255)                                                                                                                                      tGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATc             <  2:465976/150‑1 (MQ=255)                                                                                                                                        ccccGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCAt           <  1:1175288/150‑1 (MQ=255)                                                                                                                                         cccGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgc                           >  2:353765/1‑133 (MQ=255)                                                                                                                                         cccGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgc                           <  1:353765/133‑1 (MQ=255)                                                                                                                                          ctgggTGCGGCTGCACCCGGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggtgg                                                                      >  2:99427/3‑89 (MQ=38)                                                                                                                                          ctgggTGCGGCTGCACCCGGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggtgg                                                                      <  1:99427/87‑1 (MQ=38)                                                                                                                                          ctgggTGCGGCTGCACCCAGGCCACGATCTCT‑TCACCATACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGa               >  2:682179/3‑144 (MQ=255)                                                                                                                                          ctgggTGCGGCTGCACCCAGGCCACGATCTCT‑TCACCATACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGa               <  1:682179/142‑1 (MQ=255)                                                                                                                                          ccGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGt         <  1:1472167/150‑1 (MQ=255)                                                                                                                                           tgggTGCGGCTGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGGCATCCAATCACCTGCACATCGGCCACCGCCGGGTGGGGGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTc        >  1:1350752/2‑150 (MQ=18)                                                                                                                                                 cAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTtgct  <  2:560115/150‑1 (MQ=255)                                                                                                                                                    ccGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTtg    <  1:828322/145‑1 (MQ=255)                                                                                                                                                    ccGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTtg    >  2:828322/1‑145 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCT‑TCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTTGCT  >  NC_002947/3920169‑3920462

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