Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,903,787	T→C	6.0%	I381V (ATC→GTC)	PP_2556 ←	chromate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,903,787	0	T	C	6.0%	196.8 / 2.6	67	I381V (ATC→GTC)	PP_2556	chromate transporter
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (2/2);  total (35/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 = 7.61e-01

GCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGG  >  NC_002947/2903640‑2903932                                                                                                                                                    |                                                                                                                                                   gcTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATc                                                                                                                                                   <  5:245999/149‑1 (MQ=255)     gTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCgg                                                                                                                                               >  4:40900/1‑149 (MQ=255)            aaCAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTgg                                                                                                                                        >  6:298521/1‑149 (MQ=255)                   ccACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTg                                                                                                                                 <  8:145328/149‑1 (MQ=255)                     acaGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGaa                                                                                                                               <  4:324740/149‑1 (MQ=255)                     acaGCCGCGGCGAGAGCTAACACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGCCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCGCGGTCAGTCGGGCGGTGaa                                                                                                                               >  2:318156/1‑149 (MQ=255)                      ccgccgccGCGAGCGCTATCCCCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACACCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAAc                                                                                                                              <  6:218433/147‑1 (MQ=255)                          cgcgGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCa                                                                                                                          <  2:2686/149‑1 (MQ=255)                          cgcgGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCa                                                                                                                          <  2:8502/149‑1 (MQ=255)                                 aGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGtt                                                                                                                  <  3:112872/149‑1 (MQ=255)                                   cgcTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTtgt                                                                                                                >  8:254814/1‑149 (MQ=255)                                      tATCACCGCCGAGGACCAGTGGAACGCGCCTTCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCTCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTggg                                                                                                             >  6:283189/1‑149 (MQ=255)                                         cACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTc                                                                                                           >  7:187398/1‑148 (MQ=255)                                           ccgccgAGGACCAGTGGAACGCGCCATAGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCTCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGAc                                                                                                        <  5:283189/149‑1 (MQ=255)                                             gccgAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTc                                                                                                      <  5:298521/149‑1 (MQ=255)                                                  ggACCAGTGGAACGCGCCATCGAACCCGTCCGGCCCCAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTccac                                                                                                   >  1:77567/1‑147 (MQ=255)                                                  ggACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCa                                                                                                                                            <  7:220016/107‑1 (MQ=255)                                                  ggACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCa                                                                                                                                            >  8:220016/1‑107 (MQ=255)                                                   gACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTggtgg                                                                                                                             >  5:92529/1‑117 (MQ=37)                                                   gACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCCCACTGACCGGGATCACTGCGGCGGTggtgg                                                                                                                             <  6:92529/121‑5 (MQ=37)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTc                                                                                                           <  5:231545/135‑1 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTc                                                                                                           >  6:231545/1‑135 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTcca                                                                                                    >  1:220502/1‑142 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTcca                                                                                                    <  2:220502/142‑1 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCgg                                                                                             >  1:300453/1‑149 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCgg                                                                                             >  4:150483/1‑149 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGAATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTgg                                                                                                              >  2:217641/1‑132 (MQ=255)                                                      cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGAATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTgg                                                                                                              <  1:217641/132‑1 (MQ=255)                                                        gTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTc                                                                                                      <  4:100101/138‑1 (MQ=255)                                                        gTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTc                                                                                                      >  3:100101/1‑138 (MQ=255)                                                          ggAACGCGCCATCGAACCCGTCCGGCCACGGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTTCCCCCACCACCGCCGCAGTGATCCCGGTGAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCt                                                                                         <  1:318156/149‑1 (MQ=255)                                                            aaCGCGCCATCGAACCCGTCCGGCCACAGCTCGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTcc                                                                                       >  1:28824/1‑149 (MQ=255)                                                            aaCGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTcc                                                                                       >  8:282677/1‑149 (MQ=255)                                                              cgcgCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCgg                                                                                     <  7:254814/149‑1 (MQ=255)                                                                cgcCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGc                                                                                    >  4:286787/1‑148 (MQ=255)                                                                    aTCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGGGCAAGGTTGAGGATTACCCCCTCCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGa                                                                               <  3:150483/149‑1 (MQ=255)                                                                    aTCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGa                                                                               >  7:151035/1‑149 (MQ=255)                                                                       gAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGa                                                                            <  2:300453/149‑1 (MQ=255)                                                                          cccGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACa                                                                         <  3:40900/149‑1 (MQ=255)                                                                                 ggCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAAc                                                                    >  7:213192/1‑147 (MQ=255)                                                                                  gCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAg                                                                 <  6:67087/149‑1 (MQ=255)                                                                                   ccACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAgg                                                                >  3:240000/1‑149 (MQ=255)                                                                                        gCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAg                                                            >  7:303156/1‑148 (MQ=255)                                                                                                  gCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAcc                                                 <  7:282677/149‑1 (MQ=255)                                                                                                  gCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAc                                                  <  2:28824/148‑1 (MQ=255)                                                                                                  gCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAc                                                  <  8:187398/148‑1 (MQ=255)                                                                                                    cGAAGAACAGCGCAAGGTTGAGGATTAACCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGCTCGACTCAACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGCTGAACCAg                                               >  7:200962/1‑149 (MQ=255)                                                                                                     gaagaaCAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAgg                                              >  4:12531/1‑149 (MQ=255)                                                                                                        gaaCAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCa                                           >  3:108978/1‑149 (MQ=255)                                                                                                          aCAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCAcc                                         <  4:240000/149‑1 (MQ=255)                                                                                                                  aGGTGGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGc                                 <  8:303156/149‑1 (MQ=255)                                                                                                                        aGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGACGGTGAACCAGGTCACCAGGCCCGCTGCCAc                           >  5:297425/1‑149 (MQ=255)                                                                                                                        aGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAg                                                        >  4:45340/1‑120 (MQ=255)                                                                                                                        aGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAg                                                        <  3:45340/120‑1 (MQ=255)                                                                                                                         ggATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACt                          <  8:213192/149‑1 (MQ=255)                                                                                                                          gATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCATGCCCGATGCCGCTg                         >  4:275600/1‑149 (MQ=255)                                                                                                                               cccccACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGTGTCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGcc                     >  1:268662/1‑149 (MQ=255)                                                                                                                               cccccACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGccc                    >  5:139047/1‑149 (MQ=255)                                                                                                                                  ccaccaCCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGc                              >  6:339389/1‑136 (MQ=255)                                                                                                                                  ccaccaCCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGc                              <  5:339389/136‑1 (MQ=255)                                                                                                                                  ccaccaCCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCcgc                  >  2:216060/1‑148 (MQ=255)                                                                                                                                      caccGCCCCACTGACCGGGATCACTGCGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGa                                                                               >  7:181052/1‑83 (MQ=21)                                                                                                                                      caccGCCCCACTGACCGGGATCACTGCGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGa                                                                               <  8:181052/83‑1 (MQ=21)                                                                                                                                         cgccgcAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGa                                                    <  4:205668/107‑1 (MQ=255)                                                                                                                                         cgccgcAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGa                                                    >  3:205668/1‑107 (MQ=255)                                                                                                                                              cAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCtgg     >  6:52430/1‑149 (MQ=255)                                                                                                                                                  gATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCtggtgg  <  8:151035/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGGTGG  >  NC_002947/2903640‑2903932

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