Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,884,949	G→C	5.4%	I143M (ATC→ATG)	yedI ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,884,949	0	G	C	5.4%	240.1 / 4.4	74	I143M (ATC→ATG)	yedI	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base C (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

ACCAACCCATACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGC  >  NC_002947/4884803‑4885094                                                                                                                                                    |                                                                                                                                                  aCCAACCCATACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGAt                                                                                                                                                 >  4:116693/1‑149 (MQ=255)   cAACCCATACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCt                                                                                                                                               >  8:214352/1‑149 (MQ=255)     aCCCATACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTg                                                                                                                                             >  3:113537/1‑149 (MQ=255)          tACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCttt                                                                                                                                        <  3:99082/149‑1 (MQ=255)            cacacCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCttttt                                                                                                                                      <  3:116693/149‑1 (MQ=255)             acacCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCttttt                                                                                                                                      >  5:7534/1‑148 (MQ=255)               acCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCa                                                                                                                                    >  2:209185/1‑148 (MQ=255)                ccAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaa                                                                                                                                  <  2:370456/149‑1 (MQ=255)                 cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                 <  4:113537/149‑1 (MQ=255)                 cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                 <  2:132472/149‑1 (MQ=255)                 cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                 <  7:221530/149‑1 (MQ=255)                 cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                 <  4:76114/149‑1 (MQ=255)                       tcatcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGGCGAGTCAGTCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGccac                                                                                                                            >  1:86322/1‑149 (MQ=255)                         atcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAg                                                                                                                         >  1:317418/1‑149 (MQ=255)                         atcaCCACCGCAATGCCGGAAAGCACGATGAGCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGAGCTTGGCTTTTTCAAAAGccacca                                                                                                                          <  4:183935/148‑1 (MQ=255)                          tcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAgg                                                                                                                        >  7:254702/1‑149 (MQ=255)                            accaccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGt                                                                                                                       >  3:126070/1‑148 (MQ=255)                                ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGACGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAaca                                                                                                                  >  7:312686/1‑149 (MQ=255)                                ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAg                                                                                                                         <  4:245308/142‑1 (MQ=255)                                ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAg                                                                                                                         >  3:245308/1‑142 (MQ=255)                                 cGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGcca                                                                                                                             >  2:196171/1‑137 (MQ=255)                                 cGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGcca                                                                                                                             <  1:196171/137‑1 (MQ=255)                                  gCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACt                                                                                                                <  6:239094/149‑1 (MQ=255)                                         cGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATc                                                                                                          >  2:211852/1‑148 (MQ=255)                                         cGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTACAAAAGCCACCAGGTCAACACTGGCATc                                                                                                          <  1:211852/148‑1 (MQ=255)                                                            ctgGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGGTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCtt                                                                                           <  1:274176/144‑1 (MQ=255)                                                                   aGC‑GGCGAGTCAGCCCCGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTg                                                                               >  7:268692/1‑149 (MQ=255)                                                                    gc‑ggcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGt                                                                                >  4:96780/1‑147 (MQ=255)                                                                    gc‑ggcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGt                                                                                <  3:96780/147‑1 (MQ=255)                                                                     c‑ggcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCA‑CAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtg                                                                            >  3:164435/1‑149 (MQ=255)                                                                        gcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                           <  6:7534/149‑1 (MQ=255)                                                                        gcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                           <  8:117582/149‑1 (MQ=255)                                                                          gAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                         <  7:214352/149‑1 (MQ=255)                                                                           aGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGcc                                                                        >  2:140574/1‑149 (MQ=255)                                                                           aGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGcc                                                                        <  2:317418/149‑1 (MQ=255)                                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTc                                                                      <  5:104546/146‑1 (MQ=255)                                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTc                                                                      >  6:104546/1‑146 (MQ=255)                                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                   >  2:77038/1‑149 (MQ=255)                                                                                   aCGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                           <  8:118782/138‑1 (MQ=255)                                                                                   aCGATACCCAAGGTAATGGCGATGATTTCTGCCGAAAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                           >  7:118782/1‑138 (MQ=255)                                                                                     gATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCtt                                                              >  4:77071/1‑149 (MQ=255)                                                                                        aCCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCt                                                           >  1:70097/1‑149 (MQ=255)                                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                         >  2:159124/1‑149 (MQ=255)                                                                                           cAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCA‑CAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCt                                                       <  4:164435/149‑1 (MQ=255)                                                                                            aaGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCt                                                       >  7:239711/1‑149 (MQ=255)                                                                                              ggTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTCATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCacc                                                                   >  3:241706/1‑133 (MQ=255)                                                                                              ggTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTCATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCacc                                                                   >  3:236886/1‑133 (MQ=255)                                                                                              ggTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTCATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCacc                                                                   <  4:236886/135‑3 (MQ=255)                                                                                              ggTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTCATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCacc                                                                   <  4:241706/135‑3 (MQ=255)                                                                                                      cGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtg                                                                            <  6:312356/118‑1 (MQ=255)                                                                                                      cGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtg                                                                            >  5:312356/1‑118 (MQ=255)                                                                                                      cGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                   >  5:119605/1‑127 (MQ=255)                                                                                                      cGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                   <  6:119605/127‑1 (MQ=255)                                                                                                      cGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                             >  3:334067/1‑149 (MQ=255)                                                                                                           atTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTtgt                                         <  2:70097/148‑1 (MQ=255)                                                                                                           atTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTtgt                                         <  1:140574/148‑1 (MQ=255)                                                                                                               cTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgcgc                                     <  3:77071/148‑1 (MQ=255)                                                                                                                  ccGACAGGTTGAAGTCGGGACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                 <  8:312686/149‑1 (MQ=255)                                                                                                                  ccGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                         <  2:170902/109‑1 (MQ=255)                                                                                                                  ccGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                         >  1:170902/1‑109 (MQ=255)                                                                                                                  ccGACAGGAGGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAg                                  <  8:211023/148‑1 (MQ=255)                                                                                                                       aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                            <  4:126070/149‑1 (MQ=255)                                                                                                                       aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                            >  2:213404/1‑149 (MQ=255)                                                                                                                       aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                            <  1:159124/149‑1 (MQ=255)                                                                                                                               gTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgt                                                     >  7:224598/1‑116 (MQ=255)                                                                                                                               gTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgt                                                     <  8:224598/116‑1 (MQ=255)                                                                                                                                   gTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGaa                <  8:254702/149‑1 (MQ=255)                                                                                                                                        cACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGca             >  2:53180/1‑147 (MQ=255)                                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAg          <  6:162854/149‑1 (MQ=255)                                                                                                                                            gCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTTCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGaa       <  1:77038/149‑1 (MQ=255)                                                                                                                                               cccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                         >  2:228948/1‑80 (MQ=255)                                                                                                                                               cccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                         <  1:228948/80‑1 (MQ=255)                                                                                                                                                ccTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAAcgcg   >  2:308600/1‑149 (MQ=255)                                                                                                                                                  ttGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACgcgc  >  5:158011/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ACCAACCCATACACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGC‑GGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGC  >  NC_002947/4884803‑4885094

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