Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,903,780	G→C	6.3%	A383G (GCG→GGG)	PP_2556 ←	chromate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,903,780	0	G	C	6.3%	202.1 / 4.3	63	A383G (GCG→GGG)	PP_2556	chromate transporter
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

CCCCGCGCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGC  >  NC_002947/2903634‑2903916                                                                                                                                                   |                                                                                                                                          ccccGCGCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCAccgccg                                                                                                                                           >  4:112872/1‑147 (MQ=255)       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                                                                             >  4:150483/1‑149 (MQ=255)                                                            cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCgg                                                                             >  1:300453/1‑149 (MQ=255)                                                            cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGAATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTgg                                                                                              <  1:217641/132‑1 (MQ=255)                                                            cAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGAATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTgg                                                                                              >  2:217641/1‑132 (MQ=255)                                                              gTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTc                                                                                      >  3:100101/1‑138 (MQ=255)                                                              gTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTc                                                                                      <  4:100101/138‑1 (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                                  <  8:187398/148‑1 (MQ=255)                                                                                                        gCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAAc                                  <  2:28824/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)                                                                                                                                                   |                                                                                                                                          CCCCGCGCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGC  >  NC_002947/2903634‑2903916

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