Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,450,204	T→C	5.7%	L28L (TTG→CTG)	ybeY →	endoribonuclease YbeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,450,204	0	T	C	5.7%	214.7 / 2.1	70	L28L (TTG→CTG)	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTC  >  NC_002947/5450061‑5450347                                                                                                                                                |                                                                                                                                                 tATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgcgc                                                                                                                                             <  4:195212/149‑1 (MQ=255)        gAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCAGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCcagcg                                                                                                                                        >  7:149760/1‑147 (MQ=255)            ccTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCagcgcagcgc                                                                                                                                  <  8:125389/149‑1 (MQ=255)            ccTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCagcgcagcgc                                                                                                                                  >  3:267422/1‑149 (MQ=255)                  aCCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTc                                                                                                                            <  8:197414/149‑1 (MQ=255)                   ccGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGCTGGCGCAACGcc                                                                                                                                 >  7:196367/1‑143 (MQ=255)                   ccGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGCTGGCGCAACGcc                                                                                                                                 <  8:196367/143‑1 (MQ=255)                    cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCcagcgcagc                                                                                                                                    >  3:98718/1‑139 (MQ=255)                    cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCcagcgcagc                                                                                                                                    <  4:98718/139‑1 (MQ=255)                    cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGcc                                                                                                                                 >  8:155316/1‑142 (MQ=255)                    cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGcc                                                                                                                                 <  7:155316/142‑1 (MQ=255)                    cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCg                                                                                                                           >  7:202357/1‑148 (MQ=255)                      cTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGaa                                                                                                                        <  8:238434/149‑1 (MQ=255)                         cGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGaaa                                                                                                                       <  1:194920/147‑1 (MQ=255)                          gAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGa                                                                                                                    <  2:74485/149‑1 (MQ=255)                           aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGGGTTgc                                                                                                                                               <  7:194219/121‑1 (MQ=255)                           aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTgc                                                                                                                                               >  8:194219/1‑121 (MQ=255)                           aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGAc                                                                                                                   <  8:149760/149‑1 (MQ=255)                           aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGAc                                                                                                                   >  7:266825/1‑149 (MQ=255)                            aGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGAcc                                                                                                                  <  3:180020/149‑1 (MQ=255)                                    cttgataGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCt                                                                                                          <  2:154043/143‑1 (MQ=255)                                      cccAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTgg                                                                                                        >  5:222629/1‑149 (MQ=255)                                       ccAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGt                                                                                                       <  4:267422/149‑1 (MQ=255)                                        cAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCgg                                                                                                                          >  2:10022/1‑129 (MQ=255)                                        cAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCgg                                                                                                                          <  1:10022/129‑1 (MQ=255)                                        cAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTc                                                                                                      >  3:180373/1‑149 (MQ=255)                                        cAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACAGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTc                                                                                                      >  6:34217/1‑149 (MQ=255)                                           gCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTcgac                                                                                                   <  6:222629/149‑1 (MQ=255)                                               gAGGCGTCCGGCAAAGATGCTTGAACTTCATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGCAAATGACCATTCGCCTGGTCGACGAAg                                                                                               >  8:276905/1‑149 (MQ=255)                                                aGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGc                                                                                              >  1:148536/1‑149 (MQ=255)                                                                  cTTGAACTTGATATCCAACCGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGaa                                                                            <  5:34217/149‑1 (MQ=255)                                                                     gAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAc                                                                           <  7:108756/147‑1 (MQ=255)                                                                           gATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑a                                                                  >  2:60713/1‑149 (MQ=255)                                                                           gATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑a                                                                  >  4:239977/1‑149 (MQ=255)                                                                              atCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCacac                                                                      <  8:120416/143‑1 (MQ=255)                                                                              atCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCacac                                                                      >  7:120416/1‑143 (MQ=255)                                                                               tCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCTAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCg                                                               <  8:266825/148‑1 (MQ=255)                                                                                 cAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGca                                                            >  7:61468/1‑149 (MQ=255)                                                                                  aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACGCCT‑ACCGGca                                                            >  5:187865/1‑148 (MQ=255)                                                                                  aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGca                                                            >  3:148062/1‑148 (MQ=255)                                                                                  aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGca                                                            >  1:141478/1‑148 (MQ=255)                                                                                    cGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACaa                                                         <  3:250037/149‑1 (MQ=255)                                                                                     gggCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCTCCGCTGGTGCGAATTGGCGTTGCGCCATCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAg                                                        <  8:202357/149‑1 (MQ=255)                                                                                     gggCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAg                                                        <  4:180373/149‑1 (MQ=255)                                                                                       gCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGAGTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGa                                                      >  1:101445/1‑149 (MQ=255)                                                                                          aCGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAcc                                                                          >  4:3183/1‑127 (MQ=255)                                                                                          aCGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAcc                                                                          <  3:3183/127‑1 (MQ=255)                                                                                            ggATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACGCCT‑ACCGGCACAAGGACTAcg                                                 <  6:187865/149‑1 (MQ=255)                                                                                             gATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTAcgc                                                <  2:101445/149‑1 (MQ=255)                                                                                                 ctgcCCCTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGAc                                             >  4:292870/1‑148 (MQ=255)                                                                                                    ccACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAAt                                         <  4:148062/149‑1 (MQ=255)                                                                                                        tCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATgtgt                                     >  5:117754/1‑149 (MQ=255)                                                                                                          cAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATgtgt                                     <  2:141478/147‑1 (MQ=255)                                                                                                          cAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACAAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTAc                                                  <  5:272912/134‑1 (MQ=255)                                                                                                          cAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACAAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTAc                                                  >  6:272912/1‑134 (MQ=255)                                                                                                           agaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAc                                                                           >  7:38962/1‑109 (MQ=255)                                                                                                           agaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAc                                                                           <  8:38962/109‑1 (MQ=255)                                                                                                            gaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTc                                 <  7:276905/149‑1 (MQ=255)                                                                                                            gaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTc                                 <  3:292870/149‑1 (MQ=255)                                                                                                                  gccgccTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAAtgtg                                      >  4:16335/1‑138 (MQ=255)                                                                                                                  gccgccTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAAtgtg                                      <  3:16335/138‑1 (MQ=255)                                                                                                                   ccgccTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCTAACCGGCACAAGGACTACGCGACCAATGTGTTGTCGtt                              <  1:52581/146‑1 (MQ=255)                                                                                                                   ccgccTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGtt                              >  2:52581/1‑145 (MQ=255)                                                                                                                     gccTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGc                        >  5:128428/1‑149 (MQ=255)                                                                                                                          gcgcCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGAt                    <  8:61468/148‑1 (MQ=255)                                                                                                                             ccgcTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGc                >  1:181327/1‑149 (MQ=255)                                                                                                                                    tGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGa         >  7:110705/1‑149 (MQ=255)                                                                                                                                     ccgagtcggCGCTGCGCTGGCGCAACGCCAACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGAc        >  5:261455/8‑149 (MQ=255)                                                                                                                                      cgagtcggCGCTGCGCTGGCGCAACGCCAACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGAc        <  6:261455/142‑1 (MQ=255)                                                                                                                                            ggCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTc  >  5:218751/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                 TATCGTCGAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTC  >  NC_002947/5450061‑5450347

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