Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,578,340	T→A	7.3%	L161Q (CTG→CAG)	rhdA →	thiosulfate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,578,340	0	T	A	7.3%	200.7 / 5.4	68	L161Q (CTG→CAG)	rhdA	thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base A (2/3);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCA  >  NC_002947/5578194‑5578486                                                                                                                                                   |                                                                                                                                                     ccTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTgg                                                                                                                                                    >  2:170501/1‑149 (MQ=255)    cAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCa                                                                                                                                                 >  2:152068/1‑149 (MQ=255)     aaTGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCAt                                                                                                                                                >  3:109422/1‑149 (MQ=255)       tGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCt                                                                                                                                              <  5:107514/149‑1 (MQ=255)        ggCGGTATCCAGGCCTTGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTg                                                                                                                                             <  5:224752/149‑1 (MQ=255)        ggCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCAACCGAAGTGCCTGCCATTGGCGACGCCCCGGTGCCCCTGCACATAGACAGCGCACCCACTGCCCCACGCGAGTACCTGCAAAGCCCAATCGGCGCAGATGACCTGGTCATCTg                                                                                                                                             >  6:298208/1‑149 (MQ=255)         gCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCAGGTCATcgg                                                                                                                                             <  5:99059/148‑3 (MQ=255)                      cTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGcc                                                                                                                               <  1:152068/149‑1 (MQ=255)                         gCCGGCAGACAAACTGTCCACCGACGTGCCTGCCATTTGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGccc                                                                                                                              <  6:147869/147‑1 (MQ=255)                                  cAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑c                                                                                                                           >  6:228926/1‑140 (MQ=255)                                  cAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑c                                                                                                                           <  5:228926/140‑1 (MQ=255)                                     aCTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAgcg                                                                                                               >  5:44966/1‑149 (MQ=255)                                     aCTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAgcg                                                                                                               <  6:42254/149‑1 (MQ=255)                                      ctgtagaCCGAAGTGCCTGCGATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCAGCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACGTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAgcg‑g                                                                                                             <  8:128430/143‑1 (MQ=255)                                           cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGa                                                                                                                                          >  4:234062/1‑117 (MQ=255)                                           cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGa                                                                                                                                          <  3:234062/117‑1 (MQ=255)                                           cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGaaa                                                                                                        <  6:44966/149‑1 (MQ=255)                                           cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGaa                                                                                                         >  5:3713/1‑148 (MQ=255)                                           cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGaa                                                                                                         >  7:11973/1‑148 (MQ=255)                                            aCCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGCCAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCCCGCGGCCCG‑CAGGAATTCAGCG‑GCGaaa                                                                                                        >  5:240589/1‑148 (MQ=255)                                             ccGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGaaa                                                                                                        >  5:269104/1‑147 (MQ=255)                                                     cctgccATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCCAGGAATTCAGCG‑GCGAAAAAGTGCTg                                                                                                <  1:170501/148‑1 (MQ=255)                                                          cATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGcc                                                                                          >  5:89391/1‑148 (MQ=255)                                                           gttGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCa                                                                                         >  2:173050/2‑148 (MQ=255)                                                           aTTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCa                                                                                         >  2:223619/1‑148 (MQ=255)                                                           aTTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCa                                                                                         >  2:190503/1‑148 (MQ=255)                                                               gCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAggc                                                                                    >  4:89171/1‑149 (MQ=255)                                                                cAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAggcg                                                                                   >  6:200034/1‑149 (MQ=255)                                                                cAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAggcg                                                                                   <  3:89171/149‑1 (MQ=255)                                                                cAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAggcg                                                                                   <  1:223619/149‑1 (MQ=255)                                                                 aaCGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAggcgg                                                                                  >  2:67601/1‑149 (MQ=255)                                                                      cccGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCcaca                                                                             >  4:254368/1‑149 (MQ=255)                                                                         ggTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTTGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATcc                                                                          <  5:74601/149‑1 (MQ=255)                                                                                   cacaTAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGAt                                                                >  3:83966/1‑149 (MQ=255)                                                                                    acaTAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCGCGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGAt                                                                <  1:67601/149‑1 (MQ=255)                                                                                     caTAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCa                                                              >  3:216582/1‑149 (MQ=255)                                                                                        aGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                           >  4:174694/1‑149 (MQ=255)                                                                                          aCAGCGCACCCACTGCCACCCGCGAGTACCTGCAACGCCGCCTGGGCGCCTATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGTAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAACGGCGGCCACATCCCCGGCGCTATCAACTTc                                                         <  2:293016/149‑1 (MQ=255)                                                                                          aCAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTc                                                         <  5:200034/149‑1 (MQ=255)                                                                                          aCAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTc                                                         >  5:41806/1‑149 (MQ=255)                                                                                          aCAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTc                                                         <  6:3713/149‑1 (MQ=255)                                                                                             gcgcACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAg                                                      <  3:254368/149‑1 (MQ=255)                                                                                               gcACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCg                                                        <  5:107240/145‑1 (MQ=255)                                                                                               gcACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCg                                                        >  6:107240/1‑145 (MQ=255)                                                                                                 aCCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGa                                                  >  8:151112/1‑149 (MQ=255)                                                                                                     aCTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                           >  3:146077/1‑136 (MQ=255)                                                                                                     aCTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                           <  4:146077/136‑1 (MQ=255)                                                                                                      cTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGAccgcc                                             >  7:184429/1‑149 (MQ=255)                                                                                                      cTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGAccgcc                                             <  6:240589/149‑1 (MQ=255)                                                                                                                     aCCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCAG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAAc                              >  5:215537/1‑149 (MQ=255)                                                                                                                      ccTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAAc                              >  8:154326/1‑148 (MQ=255)                                                                                                                        tGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCAt                                        <  1:100301/136‑1 (MQ=255)                                                                                                                        tGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCAt                                        >  2:100301/1‑136 (MQ=255)                                                                                                                         gCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGcc                                                                                >  4:177758/1‑95 (MQ=255)                                                                                                                         gCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGcc                                                                                <  3:177758/95‑1 (MQ=255)                                                                                                                                     gggCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAACTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGAGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCAtac              >  5:214700/1‑147 (MQ=255)                                                                                                                                      gtcccagATGACCAGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCg             >  4:43634/7‑149 (MQ=255)                                                                                                                                      gtcccagATGACCAGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCg             >  8:192850/7‑149 (MQ=255)                                                                                                                                         gcCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCa          <  4:216582/149‑1 (MQ=255)                                                                                                                                          ccGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAg         <  1:17476/149‑1 (MQ=255)                                                                                                                                           cGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAg         >  6:176027/1‑148 (MQ=255)                                                                                                                                           cGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAg         >  4:287499/1‑148 (MQ=255)                                                                                                                                                aCCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATc   <  3:174694/149‑1 (MQ=255)                                                                                                                                                aCCAGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATc   <  3:43634/149‑1 (MQ=255)                                                                                                                                                aCCAGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATc   <  7:192850/149‑1 (MQ=255)                                                                                                                                                 gcTGGTCATCGGGGACGCTCGAGGCCCG‑CAGGTATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCa  <  6:215537/148‑1 (MQ=255)                                                                                                                                                 ccTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCa  <  8:184429/149‑1 (MQ=255)                                                                                                                                                  cTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCa  >  5:204265/1‑148 (MQ=255)                                                                                                                                                   tGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATc   <  3:80808/146‑1 (MQ=255)                                                                                                                                                   tGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATc   >  4:80808/1‑146 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCG‑CAGGAATTCAGCG‑GCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCA  >  NC_002947/5578194‑5578486

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