Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,903,783	G→C	6.1%	T382S (ACT→AGT)	PP_2556 ←	chromate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,903,783	0	G	C	6.1%	209.4 / 4.2	66	T382S (ACT→AGT)	PP_2556	chromate transporter
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (2/2);  total (35/31)
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.88e-01

GCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGG  >  NC_002947/2903640‑2903929                                                                                                                                                |                                                                                                                                                    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)                                                                                                                                                |                                                                                                                                                    GCTTGTAGCGGAACAGCGCCACAGCCGCGGCGAGCGCTATCACCGCCGAGGACCAGTGGAACGCGCCATCGAACCCGTCCGGCCACAGCACGTGGTAGCCGAAGAACAGCGCAAGGTTGAGGATTACCCCCACCACCGCCGCAGTGATCCCGGTCAGTGGGGCGGTGAACTTCAGT‑TCGTTGTGGGTCGACTCCACCAGCGGGCCTCCGGCGAGGATGAACAGGAACGAGGGCAGGAAGGTGAACCAGGTCACCAGGCCCGCTGCCACTGCGCCCGCGAGGAACGCCTGG  >  NC_002947/2903640‑2903929

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