Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,677,945	A→C	9.3%	A20A (GCT→GCG)	PP_4983 ←	amine oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,677,945	0	A	C	9.3%	158.0 / 2.2	58	A20A (GCT→GCG)	PP_4983	amine oxidase
Reads supporting (aligned to +/- strand):  ref base A (18/34);  new base C (4/2);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGTCAGGTCCGAG  >  NC_002947/5677800‑5678091                                                                                                                                                  |                                                                                                                                                   ccAGGACCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCAGCAAGCCTTAGTTCTGCACCTCGTCGCCTGCCGTCAGGCCTGCACGCCCCCCCCCCCCCAAAATCCCCCTGCGCGGCGGGACTTTTGCcaacccagcg                                                                                                                                                 >  4:177902/1‑149 (MQ=255)   aGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCGTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGGCTGCCCGCCCCCCCCCCACGCCAAACGCCGTGGGCGGCTTGGCTTTTGCCACCGCCGcccc                                                                                                                                               >  6:350065/1‑149 (MQ=255)      gccGAACGCCCGCCCATGCCCGCCTTGGGCTCCAGCATGGTGACCTTCCAGCCTTTTGTGTTCAGCTAGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGGGTCTTGTCTTTTGCCAGCGCAGCGCCGCt                                                                                                                                            <  3:177902/148‑1 (MQ=255)         gAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGGCGCTGCGCTggc                                                                                                                                         <  4:360646/149‑2 (MQ=255)           aCGCCCGCCCTTGCCCGCCTTGGCCTCCTGCATTGTCACCTGCCAGCCTTTGTTCTGCTGCTCGTAGGGTGCCGTGAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGaca                                                                                                                                       <  5:42689/149‑1 (MQ=255)           aCGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGGCGCTGCGCTggcaa                                                                                                                                       <  8:332350/149‑4 (MQ=255)               ccgcccATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGacacac                                                                                                                                    >  6:211973/1‑148 (MQ=255)                  cccATGCCCGCCTTAGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACcaat                                                                                                                                <  3:89304/149‑1 (MQ=255)                        cccGCCTTGGCCTCCAGCCGGGTCACCTGCCAGCCTTTGTGCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATa                                                                                                                           <  5:350065/148‑1 (MQ=255)                                ggCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCCGGCCTGCCAGCCCCCCCCCCACGACAAACGCCGTGGGCGTGATGTCTTTTGCCAGCGCAGGGCCGATGGACACACCAATCAAAACCAAcccc                                                                                                                   >  6:388827/1‑145 (MQ=255)                                     ccAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTg                                                                                                                                           >  1:400211/1‑119 (MQ=255)                                     ccAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTg                                                                                                                                           <  2:400211/119‑1 (MQ=255)                                        gcagGGGCCCCTGCCACCCTCTGTCCTTCCGCTCGTAGGCGGCCGCCCCGCCTCCCCGCCCCCCCCCCCCGACAATCGCCGTGGGCGTCTTGTCTTTTGCCATCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGgcgc                                                                                                          <  2:387536/149‑1 (MQ=255)                                          agGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCAc                                                                                                        <  7:380827/149‑1 (MQ=255)                                          agGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGGCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCAc                                                                                                        <  2:358203/149‑1 (MQ=255)                                              tCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCaa                                                                                                    >  6:339303/1‑149 (MQ=255)                                                  cTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAagcag                                                                                                >  5:296148/1‑149 (MQ=255)                                                  cTGCCAGCCGTTGTTCTGCCCCTCGTAGGGTGCCGTCAGGCCTGCCACCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAagcag                                                                                                <  8:127008/149‑1 (MQ=255)                                                     ccAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAagcagca                                                                                              <  7:83002/148‑1 (MQ=255)                                                      cAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAagcagcagc                                                                                            <  5:211973/149‑1 (MQ=255)                                                       aGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCCGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAagcagcagc                                                                                            >  4:121100/1‑148 (MQ=255)                                                           tttGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGt                                                                                       >  5:341462/1‑149 (MQ=255)                                                            ttGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTc                                                                                      >  3:145550/1‑149 (MQ=255)                                                                 cTGCAGCTCGTAGGCTTCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCg                                                                                  <  3:34234/148‑1 (MQ=255)                                                                         cGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAg                                                                         >  5:324686/1‑149 (MQ=255)                                                                          gTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGt                                                                        <  6:324686/149‑1 (MQ=255)                                                                          gTAGGCGGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGt                                                                        <  6:314259/149‑1 (MQ=255)                                                                             ggCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCag                                                                     <  2:298328/149‑1 (MQ=255)                                                                              gCTGCCGTCAGGCCTGCCAG‑CCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCagag                                                                   >  1:182823/1‑149 (MQ=255)                                                                                tGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCCCCGTGGGCGTGTTGTCTTTTGCCAGCGCAGCGCCGGTGGACCCACCAATAAATACCAACGCCCACAGGCCCACCCAAGCAGCACCCATGACGGCCAACTCCGAGTCAGAgg                                                                  >  6:142770/1‑149 (MQ=255)                                                                                  ccGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGt                                                                 >  2:372251/1‑148 (MQ=255)                                                                                    gTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTAc                                                              <  6:341462/149‑1 (MQ=255)                                                                                    gTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTAc                                                              <  5:142770/149‑1 (MQ=255)                                                                                      cAGGCCTGCCAGCCCCCCCCCCCCGACAAACGCCGTGGGCGGGTTGGGTTTTGCCAGCGCCGCGCCGCTGGACCCCCCAAACAATACCAACGCGCACAGCCGCACCCAAGCAGCAGCCATGACGCCGAACGCCCAGGCAGAGGTTa                                                               >  2:282086/1‑146 (MQ=255)                                                                                          ccTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCCGCGCCGCTGGACACCCCAATCAATACCAACCCGCACAGGCGCCCCCAAGCAGCCCCCATGTCGGCGCACGCCGAGTCAGCGGTGACAGAAAc                                                        >  6:314736/1‑149 (MQ=255)                                                                                              ccAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACggg                                                     <  1:372251/148‑1 (MQ=255)                                                                                               cAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGga                                                   <  2:63358/149‑1 (MQ=255)                                                                                                  ccccccccccACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGagag                                                <  5:388827/149‑1 (MQ=255)                                                                                                    ccccccccACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGCCACACCACTCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACTGAAACGGGGAGAGGc                                              >  7:371990/1‑149 (MQ=255)                                                                                                      ccccccACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAg                                            <  3:121100/149‑1 (MQ=255)                                                                                                          ccACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACAAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAAt                                        >  3:170178/1‑149 (MQ=255)                                                                                                           cACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATg                                       >  6:80019/1‑149 (MQ=255)                                                                                                                 gattgCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAg                                 >  6:249087/5‑149 (MQ=255)                                                                                                                 aaTCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAg                                 >  2:140299/1‑149 (MQ=255)                                                                                                                   tCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAgg                                >  4:6551/1‑148 (MQ=255)                                                                                                                        gTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATg                          <  4:145550/149‑1 (MQ=255)                                                                                                                         tGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGc                         <  3:6551/149‑1 (MQ=255)                                                                                                                           ggCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCg                        <  6:296148/148‑1 (MQ=255)                                                                                                                           ggCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGt                       <  1:282086/149‑1 (MQ=255)                                                                                                                           ggCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGt                       <  5:339303/149‑1 (MQ=255)                                                                                                                             cGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTa                     <  3:158752/149‑1 (MQ=255)                                                                                                                              gtcttgtcttTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTaa                    <  5:314736/149‑1 (MQ=255)                                                                                                                              gtcttgtcttTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTaa                    <  8:169623/149‑1 (MQ=255)                                                                                                                                   gtcttTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAAt                <  2:230530/148‑1 (MQ=255)                                                                                                                                   gtcttTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGAAATGCCCCAGGAAGATGCGTTAAGAAt                <  6:292024/148‑1 (MQ=255)                                                                                                                                   gtcttTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCACCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAAt                <  1:140299/148‑1 (MQ=255)                                                                                                                                   gtcttTTGCCAGCGCAGCGCCGCTGGACACACAAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGa                  <  4:170178/146‑1 (MQ=255)                                                                                                                                      ttttGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGNGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGt            >  1:140262/1‑149 (MQ=255)                                                                                                                                         tGCCAGCGCAGCGCCGCGGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGTCa          <  6:98872/148‑1 (MQ=255)                                                                                                                                              gcgcagcgcCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGTCAGGTcc     >  2:160808/1‑148 (MQ=255)                                                                                                                                                gcagcgcCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGTCAGGTCCGAg  >  7:213865/1‑149 (MQ=255)                                                                                                                                                 cagcgcCGATGGACACACCAATCAATACAAACGCGCACAGGCGCACCCAAGCAGAAGCCATGTCGGCGAACTACGAGTCAGAGGTTACAGAAACGGTGAGAGGCAGGAATGACCCAGTAAGATGCGTTAAGAATACGTCAGGTCCGat  >  4:37008/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCAGGCCCGAACGCCCGCCCATGCCCGCCTTGGCCTCCAGCAGGGTCACCTGCCAGCCTTTGTTCTGCAGCTCGTAGGCTGCCGTCAGGCCTGCCAGCCCCCCCCCCACGACAATCGCCGTGGGCGTCTTGTCTTTTGCCAGCGCAGCGCCGCTGGACACACCAATCAATACCAACGCGCACAGGCGCACCCAAGCAGCAGCCATGTCGGCGAACTCCGAGTCAGAGGTTACAGAAACGGGGAGAGGCAGGAATGCCCCAGGAAGATGCGTTAAGAATACGTCAGGTCCGAG  >  NC_002947/5677800‑5678091

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: ‑