Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	668,970	Δ1 bp	coding (67/318 nt)	rsfS ←	ribosomal silencing factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	668,970	0	T	.	100.0%	333.9 / NA	69	coding (67/318 nt)	rsfS	ribosomal silencing factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (26/43);  total (26/43)

CCTTCTACGCCGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAATTGT  >  NC_000913/668825‑669115                                                                                                                                                  |                                                                                                                                                  ccTTCTACGCCGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTc                                                                                                                                             <  2:396362/151‑1 (MQ=255)    tCTACGCCGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCt                                                                                                                                           <  2:522280/150‑1 (MQ=255)       aCGCCGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGa                                                                                                                                         <  1:185792/149‑1 (MQ=255)         gCCGGGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACACCGTGGTCAGCAATGGACATAACATGTCGGCTGGACGTACCCGTAGAGATGATCATGCAGTCTGTTATGCTGGATTTGCCCTGAACGTCTAAAGCGA‑GATGTACTGGCCt                                                                                                                                      <  2:1022422/150‑1 (MQ=255)          ccGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCttt                                                                                                                                    <  2:338393/151‑1 (MQ=255)            gAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTg                                                                                                                                   >  1:390711/1‑150 (MQ=255)            gAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGa                                                                                                                                  <  1:835436/151‑1 (MQ=255)             aGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGa                                                                                                                                  <  1:954042/150‑1 (MQ=255)               cGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGt                                                                                                                               <  1:207508/151‑1 (MQ=255)                   aaCAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTTGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGtcat                                                                                                                            <  1:902818/150‑1 (MQ=255)                   aaCAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGtcat                                                                                                                            <  1:32420/150‑1 (MQ=255)                     cAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCaa                                                                                                                         <  1:941992/151‑1 (MQ=255)                            gCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTc                                                                                                                  >  1:591197/1‑151 (MQ=255)                               gcgcGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCg                                                                                                                 <  2:62076/149‑1 (MQ=255)                                   gagaCTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACa                                                                                                           <  1:347149/151‑1 (MQ=255)                                         ccTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATcc                                                                                                     >  1:565133/1‑151 (MQ=255)                                             cacaACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGa                                                                                                 <  2:919563/151‑1 (MQ=255)                                              acaACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAg                                                                                                <  2:66530/151‑1 (MQ=255)                                               caACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGt                                                                                               <  1:378793/151‑1 (MQ=255)                                                aaCGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTg                                                                                              <  1:986432/151‑1 (MQ=255)                                                  cGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCt                                                                                            >  2:354394/1‑151 (MQ=255)                                                  cGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCt                                                                                            >  1:788065/1‑151 (MQ=255)                                                     ggTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCttt                                                                                          <  2:250824/150‑1 (MQ=255)                                                        cAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTAcc                                                                                       <  2:637631/150‑1 (MQ=255)                                                         aGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCt                                                                                     <  2:349285/151‑1 (MQ=255)                                                         aGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCt                                                                                     >  1:896163/1‑151 (MQ=255)                                                          gCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTg                                                                                    <  1:442889/151‑1 (MQ=255)                                                               ggACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAg                                                                                <  2:715387/150‑1 (MQ=255)                                                               ggACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGt                                                                               <  1:727455/151‑1 (MQ=255)                                                                     aaCATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTccc                                                                         >  1:892822/1‑151 (MQ=255)                                                                     aaCATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTcc                                                                          <  2:313078/150‑1 (MQ=255)                                                                      aCATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTcccc                                                                        <  2:986025/151‑1 (MQ=255)                                                                          gACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTgg                                                                    <  1:1027612/151‑1 (MQ=255)                                                                             ggCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTg                                                                 <  2:1020360/151‑1 (MQ=255)                                                                                tGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAAt                                                              >  1:991679/1‑151 (MQ=255)                                                                                tGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAAt                                                              >  2:156899/1‑151 (MQ=255)                                                                                 ggACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATc                                                             >  1:953653/1‑151 (MQ=255)                                                                                 ggACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATc                                                             <  1:602508/151‑1 (MQ=255)                                                                                   aCGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaa                                                           <  2:739983/151‑1 (MQ=255)                                                                                   aCGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaa                                                           <  1:378606/151‑1 (MQ=255)                                                                                    cGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCaaa                                                          <  2:448309/151‑1 (MQ=255)                                                                                     gTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAAt                                                         <  2:351839/151‑1 (MQ=255)                                                                                       aCCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAAt                                                         >  2:755688/1‑149 (MQ=255)                                                                                       aCCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAAt                                                         >  1:695754/1‑149 (MQ=255)                                                                                        cccGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAAt                                                         >  1:241809/1‑148 (MQ=255)                                                                                          cGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGAt                                                     >  2:180371/1‑150 (MQ=255)                                                                                          cGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATa                                                    <  2:875564/151‑1 (MQ=255)                                                                                           gTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATa                                                    >  2:1002903/1‑150 (MQ=255)                                                                                            tACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATaa                                                   >  2:250148/1‑150 (MQ=255)                                                                                                     tCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg                                         >  2:499677/1‑151 (MQ=255)                                                                                                      cATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg                                         <  1:336248/150‑1 (MQ=255)                                                                                                      cATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg                                         <  1:416635/150‑1 (MQ=255)                                                                                                      cATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTg                                         <  1:738589/150‑1 (MQ=255)                                                                                                         gCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCa                                      <  2:309720/150‑1 (MQ=255)                                                                                                              cGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATAc                                <  2:855292/151‑1 (MQ=255)                                                                                                                    tGCGGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGTTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGaa                           <  2:199691/150‑1 (MQ=255)                                                                                                                       tGGATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAAcc                         >  2:852198/1‑149 (MQ=255)                                                                                                                         gATTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCagag                     <  2:565125/151‑1 (MQ=255)                                                                                                                          aTTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAgg                    <  1:625325/151‑1 (MQ=255)                                                                                                                          aTTTGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAgg                    >  1:256755/1‑151 (MQ=255)                                                                                                                             tGCCCTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCg                  >  1:172003/1‑150 (MQ=255)                                                                                                                               cccTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGAtt               >  2:750384/1‑151 (MQ=255)                                                                                                                                ccTGAACGTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGAttt              >  2:697271/1‑151 (MQ=255)                                                                                                                                       gTCTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACa       <  1:156900/151‑1 (MQ=255)                                                                                                                                         cTAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAAt     >  1:1026602/1‑151 (MQ=255)                                                                                                                                          tAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAAtt    >  1:780806/1‑151 (MQ=255)                                                                                                                                          tAAGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAAtt    >  1:99838/1‑151 (MQ=255)                                                                                                                                            aGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAATTGt  <  2:317896/151‑1 (MQ=255)                                                                                                                                            aGGCGA‑GATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAATTGt  >  1:662990/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCTTCTACGCCGAGCGGTAACAGGCCCGCTGCGCGAGACTCCTGCACAACGTGGTCAGCAATGGACATAACATGACGGCTGGACGTACCCGTACAGATGATCATGCAGTCGGTGATGCTGGATTTGCCCTGAACGTCTAAGGCGATGATGTCCTGACCTTTGAGGTCATCAATTTTGTCGATAACAAAATCCTGGAGTGCTTTACCCTGCAAGTTTTCCCCCTGGGTGAATCAAATAGATAAAAATGGTCTGTCAGTATACCTGAACCAGAGGCGATTTCGGGACAATTGT  >  NC_000913/668825‑669115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑