Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,450,210	C→T	6.0%	Q30* (CAG→TAG)	ybeY →	endoribonuclease YbeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,450,210	0	C	T	6.0%	199.1 / 2.4	66	Q30* (CAG→TAG)	ybeY	endoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGC  >  NC_002947/5450068‑5450357                                                                                                                                               |                                                                                                                                                     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                                                                                                                                           <  8:196367/143‑1 (MQ=255)            ccGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGCTGCGCTGGCGCAACGcc                                                                                                                                           >  7:196367/1‑143 (MQ=255)             cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCcagcgcagc                                                                                                                                              <  4:98718/139‑1 (MQ=255)             cGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCcagcgcagc                                                                                                                                              >  3:98718/1‑139 (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)                    aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGAc                                                                                                                             >  7:266825/1‑149 (MQ=255)                    aaGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGAc                                                                                                                             <  8:149760/149‑1 (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                                                                                                                                    <  1:10022/129‑1 (MQ=255)                                 cAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCgg                                                                                                                                    >  2:10022/1‑129 (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                                                                            >  4:239977/1‑149 (MQ=255)                                                                    gATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑a                                                                            >  2:60713/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                                                                      >  1:141478/1‑148 (MQ=255)                                                                           aaCGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGca                                                                      >  3:148062/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                                                                                    <  3:3183/127‑1 (MQ=255)                                                                                   aCGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAcc                                                                                    >  4:3183/1‑127 (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                                                            >  6:272912/1‑134 (MQ=255)                                                                                                   cAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACAAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTAc                                                            <  5:272912/134‑1 (MQ=255)                                                                                                    agaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAc                                                                                     <  8:38962/109‑1 (MQ=255)                                                                                                    agaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAAc                                                                                     >  7:38962/1‑109 (MQ=255)                                                                                                     gaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTc                                           <  3:292870/149‑1 (MQ=255)                                                                                                     gaagaCGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTc                                           <  7:276905/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)                                                                                                                                            cgcCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATcc     <  4:290165/148‑1 (MQ=255)                                                                                                                                             gcCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCg   >  5:94382/1‑149 (MQ=255)                                                                                                                                              ccAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCgc  <  8:110705/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     GAAGCCTACGACCGCTTCGAAGCCCGCCAGCCCAAGCCCGAGGCGTCCGGCAAAGATGCTTGAACTTGATATCCAACGGGCCACGGATGCTGCCACTCCAGAAGACGCCGCCTTGCGCCGCTGGTGCGAATTGGCGTTGCGCCAGCGCAGCGCCGACTCGGAAATGACCATTCGCCTGGTCGACGAAGCCGAAGGCCGCGAGCTGAACCACACCT‑ACCGGCACAAGGACTACGCGACCAATGTGTTGTCGTTCCCGGCCGATGTGCCCGATGACCTGCTCGATATCCCGC  >  NC_002947/5450068‑5450357

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