Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,550,896	C→G	6.7%	R89G (CGT→GGT)	pydA →	NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,550,896	0	C	G	6.7%	226.8 / 3.8	71	R89G (CGT→GGT)	pydA	NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base G (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCC  >  NC_002947/4550750‑4551043                                                                                                                                                   |                                                                                                                                                      gcggcgTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGt                                                                                                                                                     <  5:336850/149‑1 (MQ=255) gcggcgTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGt                                                                                                                                                     <  4:192900/149‑1 (MQ=255)     cgTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGaaa                                                                                                                                                 <  6:90021/149‑1 (MQ=255)              gAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAgg                                                                                                                                        >  2:227614/1‑149 (MQ=255)              gAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGCTACAGGGGATCAACAATATCGAGTTGATCACCGCCCGCTCTATGGAAATCAACCTGCGTGAAATCACCCAgg                                                                                                                                        >  2:190310/1‑149 (MQ=255)                  aCCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTg‑aa                                                                                                                                   <  2:76794/149‑1 (MQ=255)                    ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCAGGAAATCAACCTGCGTGAAATCACCCAGGTg‑aag                                                                                                                                  >  4:280944/1‑148 (MQ=255)                        ggCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGa                                                                                                                             >  5:141894/1‑149 (MQ=255)                            aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTgg                                                                                                                         <  7:130414/149‑1 (MQ=255)                            aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCGCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTgg                                                                                                                         >  2:271674/1‑149 (MQ=255)                             aGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGc                                                                                                                        <  6:286890/149‑1 (MQ=255)                              gACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGcc                                                                                                                       >  8:297670/1‑149 (MQ=255)                                 ccTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGa                                                                                                                    >  2:227187/1‑149 (MQ=255)                                      ggcggTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACcgcg                                                                                                               <  4:38201/149‑1 (MQ=255)                                        cggTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCa                                                                                                             <  1:227614/149‑1 (MQ=255)                                             aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACACGGGAACAACAATATCGAGTTGATCACCGACCGCTCACTGGAAATCAACCTGCCAGAAATCACCCAGGTG‑AAGAAGGACAGGCCCGACCGCGCACTCAt                                                                                                        >  2:299993/1‑149 (MQ=255)                                                gtgtCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCgt                                                                                                     <  1:227187/149‑1 (MQ=255)                                                        gcgcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                             <  1:202813/149‑1 (MQ=255)                                                        gcgcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                             <  7:85321/149‑1 (MQ=255)                                                          gcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATg                                                                                           >  7:282892/1‑149 (MQ=255)                                                              cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGt                                                                                         <  6:262783/147‑1 (MQ=255)                                                              cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGc                                                                                       <  6:141894/149‑1 (MQ=255)                                                                    gcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGCAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGc                                                                                  <  1:258605/149‑1 (MQ=255)                                                                    gcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCa                                                                                 >  1:64519/1‑149 (MQ=255)                                                                      aCTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATa                                                                               <  2:279726/149‑1 (MQ=255)                                                                          tGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATaga                                                                             <  2:64519/147‑1 (MQ=255)                                                                             ccccAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAt                                                                        >  2:208248/1‑149 (MQ=255)                                                                               ccAACCGTCTGGTACAGGGGATCAACAATATCGAGTTTATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑Cg                                                                     <  1:190310/149‑1 (MQ=255)                                                                                   ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGa                                                                 <  7:297670/149‑1 (MQ=255)                                                                                     gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAg                                                               <  6:168573/149‑1 (MQ=255)                                                                                     gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAg                                                               >  5:168573/1‑149 (MQ=255)                                                                                     gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAg                                                               <  1:208248/149‑1 (MQ=255)                                                                                      tCAGGTACAGGGGATCAACAATGTCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGt                                                              >  5:316237/1‑149 (MQ=255)                                                                                             cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCt                                                       >  3:311640/1‑149 (MQ=255)                                                                                             cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCt                                                       >  3:106185/1‑149 (MQ=255)                                                                                                gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCAc                                                                                                            >  5:6728/1‑94 (MQ=255)                                                                                                gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCAc                                                                                                            <  6:6728/94‑1 (MQ=255)                                                                                                       acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                            <  4:139506/119‑1 (MQ=255)                                                                                                       acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                            >  3:139506/1‑119 (MQ=255)                                                                                                       acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCtggtg                                             >  6:141806/1‑149 (MQ=255)                                                                                                       acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCtggtg                                             >  1:116303/1‑149 (MQ=255)                                                                                                             tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  4:201525/1‑148 (MQ=255)                                                                                                             tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  7:205333/1‑148 (MQ=255)                                                                                                             tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  5:162645/1‑148 (MQ=255)                                                                                                              cGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                      >  6:179171/1‑149 (MQ=255)                                                                                                                      tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                              >  1:253704/1‑149 (MQ=255)                                                                                                                      tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                              >  8:176212/1‑149 (MQ=255)                                                                                                                          cgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AATAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                          <  4:32209/149‑1 (MQ=255)                                                                                                                              cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCAt                      <  8:282892/149‑1 (MQ=255)                                                                                                                               gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggg                                  >  5:324319/1‑136 (MQ=255)                                                                                                                               gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                     >  1:317078/1‑149 (MQ=255)                                                                                                                               gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGGGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggg                                  <  6:324319/136‑1 (MQ=255)                                                                                                                                 tctctGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGAGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                     >  4:123969/1‑147 (MQ=255)                                                                                                                                   tctGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGa                   <  1:299993/147‑1 (MQ=255)                                                                                                                                   tctGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 >  2:102814/1‑149 (MQ=255)                                                                                                                                   tctGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 >  5:181722/1‑149 (MQ=255)                                                                                                                                     tGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 <  7:176212/147‑1 (MQ=255)                                                                                                                                     tGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                >  6:274632/1‑148 (MQ=255)                                                                                                                                     tGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                >  7:125412/1‑148 (MQ=255)                                                                                                                                        ttcttcaCCTGGGTGATTTCACGCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAAc            >  4:236150/7‑149 (MQ=255)                                                                                                                                          aTCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCCCTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACt           <  1:271674/149‑1 (MQ=255)                                                                                                                                             caCCTGGGTGATTTCACGCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCgg       >  5:37094/2‑149 (MQ=255)                                                                                                                                             caCCTGGGTGATTTCACGCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCgg       >  4:91447/2‑149 (MQ=255)                                                                                                                                              aCCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGt      <  6:181722/149‑1 (MQ=255)                                                                                                                                               ccTGGGTGATTTCACGCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt     <  6:37094/149‑1 (MQ=255)                                                                                                                                               ccTGGGTGATTTCACGCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt     <  3:236150/149‑1 (MQ=255)                                                                                                                                                cTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCtggt                                              >  6:205838/1‑107 (MQ=255)                                                                                                                                                cTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCtggt                                              <  5:205838/107‑1 (MQ=255)                                                                                                                                                 tGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTgg                                                                  <  2:132257/86‑1 (MQ=255)                                                                                                                                                 tGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTgg                                                                  >  1:132257/1‑86 (MQ=255)                                                                                                                                                 tGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTAATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGc   >  7:315885/1‑149 (MQ=255)                                                                                                                                                   cGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc  <  2:317078/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCGGCGTGGTCTGGAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTG‑AAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAT‑CGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCC  >  NC_002947/4550750‑4551043

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