Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	132,939	A→T	6.1%	Q121L (CAG→CTG)	dsbA →	thiol:disulfide interchange protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	132,939	0	A	T	6.1%	202.6 / 2.9	66	Q121L (CAG→CTG)	dsbA	thiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (2/2);  total (34/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.89e-01

AAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCGAAGAAGTATGAGAT  >  NC_002947/132792‑133086                                                                                                                                                    |                                                                                                                                                    aaaaGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCa                                                                                                                                                     <  3:226273/148‑1 (MQ=255) aaaaGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCa                                                                                                                                                     <  5:387757/148‑1 (MQ=255) aaaaGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCa                                                                                                                                                     <  5:381004/148‑1 (MQ=255)  aaaGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAg                                                                                                                                                    >  8:298286/1‑148 (MQ=255)  aaaGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGc                                                                                                                                                   <  5:64241/149‑1 (MQ=255)               gACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGAc                                                                                                                                      <  6:189116/149‑1 (MQ=255)                 cGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGAcc                                                                                                                                     >  4:172134/1‑148 (MQ=255)                   tCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGAt                                                                                                                                  <  3:172134/149‑1 (MQ=255)                   tCAACTTCAAGCGCGTACCGGACATGTTCGGTGGCCGGTGGGTTGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGAt                                                                                                                                  <  2:356305/149‑1 (MQ=255)                    cAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGa                                                                                                                                       >  5:32645/1‑143 (MQ=255)                    cAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGa                                                                                                                                       <  6:32645/143‑1 (MQ=255)                          cAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAAc                                                                                                                            <  2:30775/148‑1 (MQ=255)                           aaGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAAc                                                                                                                            >  2:65531/1‑147 (MQ=255)                            aGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACgag                                                                                                                         <  1:285886/149‑1 (MQ=255)                               gcgTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGCCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATg                                                                                                                      >  7:119509/1‑149 (MQ=255)                                            gTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCg                                                                                                         >  3:286803/1‑149 (MQ=255)                                                    gCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGCCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCaa                                                                                                 <  8:119509/149‑1 (MQ=255)                                                     cccGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAg                                                                                                <  4:286803/149‑1 (MQ=255)                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTc                                                                                                   >  5:146171/1‑140 (MQ=255)                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTc                                                                                                   <  6:146171/140‑1 (MQ=255)                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGtgtg                                                                                           >  2:308327/1‑148 (MQ=255)                                                           ggATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGtgtg                                                                                           >  2:292121/1‑148 (MQ=255)                                                            gATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTgga                                                                                         <  5:307369/149‑1 (MQ=255)                                                                 acacGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaag                                                                                     >  8:333353/1‑148 (MQ=255)                                                                   acGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaaggac                                                                                  <  7:333353/149‑1 (MQ=255)                                                                     ggCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaaggacaa                                                                                <  7:217523/149‑1 (MQ=255)                                                                         aGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGtt                                                                             >  4:292264/1‑148 (MQ=255)                                                                                 cTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTTCTCGCCACTCAAGGTGTGGACAAGGACAAGTTcc                                                                           >  3:375314/1‑142 (MQ=255)                                                                                 cTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTTCTCGCCACTCAAGGTGTGGACAAGGACAAGTTcc                                                                           <  4:375314/142‑1 (MQ=255)                                                                                 cTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCAc                                                                    <  7:298286/149‑1 (MQ=255)                                                                                  tGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCAcc                                                                   >  3:110175/1‑149 (MQ=255)                                                                                     cccTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTc                                                                <  1:65531/149‑1 (MQ=255)                                                                                                aTGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGc                                                     >  3:234992/1‑149 (MQ=255)                                                                                                aTGCGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGc                                                     <  5:300993/149‑1 (MQ=255)                                                                                                  gggCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCa                                                   >  3:365920/1‑149 (MQ=255)                                                                                                    gCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATc                                                 <  3:292264/149‑1 (MQ=255)                                                                                                            cacaAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGc                                           >  5:270053/1‑147 (MQ=255)                                                                                                            cacaAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCa                                         <  1:292121/149‑1 (MQ=255)                                                                                                             acaAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAg                                        >  1:373657/1‑149 (MQ=255)                                                                                                                  gTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCaa                                   >  2:333526/1‑149 (MQ=255)                                                                                                                   atcttgctgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCaaa                                  <  7:53801/145‑1 (MQ=255)                                                                                                                    tCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAACCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGAc                                                             <  2:114945/121‑1 (MQ=255)                                                                                                                    tCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAACGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGAc                                                             >  1:114945/1‑121 (MQ=255)                                                                                                                      aTGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAg                               <  2:373657/149‑1 (MQ=255)                                                                                                                       tgctgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAgg                              >  1:280576/1‑149 (MQ=255)                                                                                                                       tgctgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAgg                              >  6:349563/1‑149 (MQ=255)                                                                                                                        gatgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGACATCAAGGGCCAGGTCaaa                                  >  8:53801/3‑144 (MQ=255)                                                                                                                          tgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTc                                     >  6:256537/1‑139 (MQ=255)                                                                                                                          tgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTc                                     <  5:256537/139‑1 (MQ=255)                                                                                                                          tgcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGcc                            >  5:325345/1‑148 (MQ=255)                                                                                                                           gcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCt                                                                                                           >  1:42323/1‑68 (MQ=255)                                                                                                                           gcGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCt                                                                                                           <  2:42323/68‑1 (MQ=255)                                                                                                                            cGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAg                         <  6:270053/149‑1 (MQ=255)                                                                                                                            cGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAg                         <  2:280576/149‑1 (MQ=255)                                                                                                                              gtgtTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAgg                        >  8:92143/1‑148 (MQ=255)                                                                                                                                      tGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCg                >  2:108555/1‑148 (MQ=255)                                                                                                                                         cATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGTCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCgaag             >  6:163669/1‑148 (MQ=255)                                                                                                                                          aTCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCgaagaa           <  4:365920/149‑1 (MQ=255)                                                                                                                                          aTCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCgaagaa           <  5:349563/149‑1 (MQ=255)                                                                                                                                          aTCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCgaagaa           <  7:92143/149‑1 (MQ=255)                                                                                                                                            ccAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCGaagaag          >  5:320595/1‑148 (MQ=255)                                                                                                                                             cAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaaggac                                                                                  <  6:383897/75‑1 (MQ=255)                                                                                                                                             cAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTggacaaggac                                                                                  >  5:383897/1‑75 (MQ=255)                                                                                                                                               gAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCTGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCGAAGAAGTAt       >  5:86729/1‑148 (MQ=255)                                                                                                                                                aaCCTGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGgacaaggacaag                                                                               >  7:113213/1‑75 (MQ=255)                                                                                                                                                aaCCTGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGgacaaggacaag                                                                               >  7:107902/1‑75 (MQ=255)                                                                                                                                                aaCCTGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGgacaaggacaag                                                                               <  8:107902/75‑1 (MQ=255)                                                                                                                                                aaCCTGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGgacaaggacaag                                                                               <  8:113213/75‑1 (MQ=255)                                                                                                                                                 aCCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAACTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCGAAGAAGTAt       >  5:48848/1‑146 (MQ=255)                                                                                                                                                    aGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAACCAGGCCAAGGAATTGGCGAAGAAGTATGAGAt  >  8:207932/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAAAGCTGCCGAAGGACGTCAACTTCAAGCGCGTACCGGCCATGTTCGGTGGCCCGTGGGATGCACACGGCCAGATGTTCCTGACCCTCGAAGCCATGGGCGTCGAGCACAAGGTTCATGCTGCGGTGTTCGATGCCATCCAGAACCAGCACAAGCGTCTGACCGATAAGAACGAGATGGCCGACTTCCTCGCCACTCAAGGTGTGGACAAGGACAAGTTCCTCGCCACCTTCGACTCGTTCGCCATCAAGGGCCAGGTCAAACAGGCCAAGGAATTGGCGAAGAAGTATGAGAT  >  NC_002947/132792‑133086

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