Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,877,162	A→T	8.6%	Q84L (CAG→CTG)	PP_3425 →	RND family transporter MFP subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,877,162	0	A	T	8.6%	197.4 / 6.5	70	Q84L (CAG→CTG)	PP_3425	RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base T (2/4);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGTGACA  >  NC_002947/3877016‑3877308                                                                                                                                                    |                                                                                                                                                    gcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCT‑TCGGTGAAGGCGACCTgg                                                                                                                                                   <  1:319000/149‑1 (MQ=255)  ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCT‑TCGGTGAAGGCGACCTGGt                                                                                                                                                  <  3:310324/149‑1 (MQ=255)  ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGt                                                                                                                                                  >  8:93948/1‑149 (MQ=255)  ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACAGTCGAAGTGCGCCCGCGCGTCTCCGGCCAGAACGACACGGTCGACT‑TCAGCGAAGGCGCACAGGt                                                                                                                                                  >  6:346645/1‑149 (MQ=255)     aaGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCaa                                                                                                                                               <  7:73304/149‑1 (MQ=255)        gTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACAGGCCGCCTCGAAGCCCCGGAGACAGTCGAAGTGCGCCCGCGGTTCTCCGGCCAGATCGAACAGGTCGCCT‑TCACCGAAGTCGCCCAGGTCAAGa                                                                                                                                             >  3:225677/1‑148 (MQ=255)         tAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGaa‑a                                                                                                                                          <  7:179629/149‑1 (MQ=255)           aGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑Ag                                                                                                                                         <  5:200636/148‑1 (MQ=255)           aGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑Ag                                                                                                                                         <  6:114722/148‑1 (MQ=255)               tgGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGAc                                                                                                                                    <  2:137704/149‑1 (MQ=255)               tgGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGAc                                                                                                                                    <  7:93948/149‑1 (MQ=255)                 gCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCT‑TCGGTGAAGGCGACCTGGTc                                                                                                                                                 <  5:138087/135‑1 (MQ=38)                 gCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCT‑TCGGTGAAGGCGACCTGGTc                                                                                                                                                 >  6:138087/1‑135 (MQ=38)                     aGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACctgctg                                                                                                                              <  2:21544/149‑1 (MQ=255)                      ggTGATCGAACAACCGATCACCAAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACctgctg                                                                                                                              >  7:55323/1‑148 (MQ=255)                          aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCa                                                                                                                         <  8:45958/149‑1 (MQ=255)                             ggttgttcGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGa                                                                                                                       <  3:257266/141‑1 (MQ=255)                             gAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGTCACCGACGAAGTGCGCCCCCGCGTCTCCGGCCCGATCGACCACGTCGCAT‑TCGCCGAAGGCGCGCAGGTCATCGA‑AGGCTCCGTGCTGGTCCAgg                                                                                                                       >  8:335625/1‑147 (MQ=255)                             gAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCTT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGa                                                                                                                       <  7:142788/148‑1 (MQ=255)                               acaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTGTCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATc                                                                                                                     <  1:79818/149‑1 (MQ=255)                                     gATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCcgc                                                                                                              >  3:63615/1‑149 (MQ=255)                                          ccGACTGGGACGAGTGCACCGGCCGCCTCGAATCCCCGGAAACCTTCGAAGTGCGCCCGCGTGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTc                                                                                                                                                 <  6:167849/110‑1 (MQ=255)                                          ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGACACCCACGAAGTGCGCCCGCGGGTCGGCGGCGCGATCGACCCGCTCCCAG‑TCACCGACGTCGCCGACGTCAAGAA‑ACGCGACCGCCTGTTCCAGATAGACCCGCGACCg                                                                                                         >  6:213420/1‑149 (MQ=255)                                          ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCg                                                                                                         >  3:347220/1‑149 (MQ=255)                                          ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGACAGAGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTc                                                                                                                                                 >  5:167849/1‑110 (MQ=255)                                             aaTGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTc                                                                                                      <  7:293854/149‑1 (MQ=255)                                                 ggACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAgg                                                                                                  <  8:55323/149‑1 (MQ=255)                                                      aGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGaa                                                                                             <  2:145336/149‑1 (MQ=255)                                                          cACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTcc                                                                                         >  4:315445/1‑149 (MQ=255)                                                           aCCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                               >  3:183860/1‑142 (MQ=255)                                                           aCCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                               <  4:183860/142‑1 (MQ=255)                                                            ccGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCC‑TCACCGAAGGCGCCCAGGTCAAGAA‑ACGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTccgc                                                                                       <  6:163061/149‑1 (MQ=255)                                                              ggCCGCCTCGAAGCCCCGGAACCCCTCGAAGTGCGCCCGCGGCTCGCCGGCCAGATCGACCAGGTCGGGG‑TCCCCGAAGGCGCCCAGGGCAAGCA‑ACGCGCCCTGCAGTTCCAGCACGACCCGCGCCCTGTCCAGGCTGAAGAcctccg                                                                                     >  2:221526/1‑149 (MQ=255)                                                                ccgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCgccgcc                                                                                   <  3:315445/149‑1 (MQ=255)                                                                 cgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGCCCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCt                                                                                  >  5:222134/1‑149 (MQ=255)                                                                       gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGTTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCAGTTCCAGGATGAATTCCGCCGCCTCGaa                                                                              <  2:146346/147‑1 (MQ=255)                                                                       gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                              <  4:88718/147‑1 (MQ=255)                                                                       gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGc                                                                            >  2:115552/1‑149 (MQ=255)                                                                             ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCg                                                                                                                                      <  6:341983/85‑1 (MQ=255)                                                                             ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCg                                                                                                                                      >  5:341983/1‑85 (MQ=255)                                                                             ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                              >  3:3252/1‑141 (MQ=255)                                                                             ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                              <  4:3252/141‑1 (MQ=255)                                                                                   cgggcataaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACa                                                                <  3:53087/142‑1 (MQ=255)                                                                                          aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaag                                                         >  8:325937/1‑149 (MQ=255)                                                                                            gTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaaggc                                                       >  8:251411/1‑149 (MQ=255)                                                                                                gcCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTTAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAcc                                                   <  1:44023/149‑1 (MQ=255)                                                                                                     cgGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGCCCCGCGCCCGTTCCAGGCGGAAGTCCGCCGCCTCGAAGCCCAACTGCAACag                                                               >  1:91652/1‑132 (MQ=255)                                                                                                      gggTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGAGTCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                             <  4:63615/149‑1 (MQ=255)                                                                                                         tctcCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCg                                           >  3:319169/1‑148 (MQ=255)                                                                                                            ccGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                       >  2:77787/1‑149 (MQ=255)                                                                                                            ccGGCCAGATCGACCAGGTCACCT‑TCACCGAAGGCGCCCAGGTCAAGAAGAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAg                                        >  1:85214/1‑149 (MQ=255)                                                                                                               gCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                             >  4:189511/1‑140 (MQ=255)                                                                                                               gCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                             <  3:189511/140‑1 (MQ=255)                                                                                                                 cAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGTTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCCCGCCCGTTCCAGGCTTAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCaa                                  <  5:346645/149‑1 (MQ=255)                                                                                                                        aCCTGGGCGCCT‑TCGGTGAAGGCGACCTGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                       >  4:112416/1‑137 (MQ=255)                                                                                                                        aCCTGGGCGCCT‑TCGGTGAAGGCGACCTGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                       <  3:112416/137‑1 (MQ=255)                                                                                                                        aCCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGcc                           >  3:39553/1‑149 (MQ=255)                                                                                                                                gCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTgg                   <  1:77787/149‑1 (MQ=255)                                                                                                                                  cT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                          >  1:178180/1‑108 (MQ=255)                                                                                                                                  cT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                          <  2:178180/108‑1 (MQ=255)                                                                                                                                  cT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc                  >  2:211915/1‑148 (MQ=255)                                                                                                                                   t‑tCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc                  >  8:58799/1‑147 (MQ=255)                                                                                                                                      cACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAg               <  1:211915/148‑1 (MQ=255)                                                                                                                                       aCCGAAGGCGCCCATGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGAGCGGCGCGATTTGCAGGCTGAATTCGGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgc              <  2:87953/148‑1 (MQ=255)                                                                                                                                       aCCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgcg             <  7:251411/149‑1 (MQ=255)                                                                                                                                           aaGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTg         >  4:114031/1‑149 (MQ=255)                                                                                                                                            aGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGc        <  7:58799/149‑1 (MQ=255)                                                                                                                                            aGGCGACCAGGTCAAGAA‑AGTCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGc        <  6:222134/149‑1 (MQ=255)                                                                                                                                              gcgcCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGt      >  2:192789/1‑149 (MQ=255)                                                                                                                                                  ccAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGTGACa  <  4:39553/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCT‑TCACCGAAGGCGCCCAGGTCAAGAA‑AGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGTGACA  >  NC_002947/3877016‑3877308

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