Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,346	T→G	5.4%	T13P (ACG→CCG)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,346	0	T	G	5.4%	260.1 / 4.3	88	T13P (ACG→CCG)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base G (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

CGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGG  >  NC_002947/6199‑6490                                                                                                                                                    |                                                                                                                                                      cGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGt                                                                                                                                                       <  5:234273/148‑1 (MQ=255)    ccAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                    >  4:15233/1‑148 (MQ=255)     cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCg                                                                                                                                                   <  4:302180/148‑1 (MQ=255)     cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGa                                                                                                                                                  >  8:313341/1‑149 (MQ=255)     cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGa                                                                                                                                                  >  6:291004/1‑149 (MQ=255)       gCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGTAATGCCGCGGGGTCAGGGTTCGGCCGGTGATCATCAGGCCGGCCTTGGCAGCCAGAGGGCCTGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCGGTg                                                                                                                                                <  2:241726/149‑1 (MQ=255)        cACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATgg                                                                                                                                               <  5:75706/149‑1 (MQ=255)         aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggc                                                                                                                                              >  7:183430/1‑149 (MQ=255)           cAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcgg                                                                                                                                            >  3:251510/1‑149 (MQ=255)             gCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcggc                                                                                                                                           >  1:204801/1‑148 (MQ=255)             gCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcggc                                                                                                                                           >  5:6681/1‑148 (MQ=255)             gCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcggc                                                                                                                                           >  3:211974/1‑148 (MQ=255)             gCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGgcggcg                                                                                                                                          <  8:126530/149‑1 (MQ=255)             gCCCTTCATCGTCGCGGAATGGAACATAATGGGAATGCCGCGGGGTAAGGGTGCGACCGGAAAACAAAAGGCCGGCATTGACAGCCAGAGGGCCGAAAACCATAAAAATACCCACACACCCCCGgccgtgggccgtgg                                                                                                                                                     >  7:235141/1‑138 (MQ=255)             ccccTTCATCGTCGCGGAATGTACCGTAATGGGCATGCCGCGGGTTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTgg                                                                                                                                                     <  8:235141/137‑1 (MQ=255)              cccTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGa                                                                                                                                         <  4:194550/149‑1 (MQ=255)               ccTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGAt                                                                                                                                        <  3:178331/149‑1 (MQ=255)                 ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGAACAACAGGCCGGCCATGGCAGCCAGAGGGCCGGACACCCTGACAATTCCCCCACAG‑CCCGGCCCTGTGCCGTGGCGCTGGCGGCGATCGt                                                                                                                                     >  8:141021/1‑149 (MQ=255)                       gtcgCGGAATGGACTCTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                <  1:317206/149‑1 (MQ=255)                        tcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACg                                                                                                                               >  7:156320/1‑149 (MQ=255)                            ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                           <  2:97063/149‑1 (MQ=255)                            ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                           <  7:157032/149‑1 (MQ=255)                               aTGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACTGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgtg                                                                                                                        <  1:70096/149‑1 (MQ=255)                                  gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGca                                                                                                                             >  6:18091/1‑141 (MQ=255)                                  gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGca                                                                                                                             <  5:18091/141‑1 (MQ=255)                                    ccGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGTCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCGGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCAt                                                                                                                   >  8:280350/1‑149 (MQ=255)                                       tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCt                                                                                                                >  1:310167/1‑149 (MQ=255)                                         aTGGGCATGCCGCGGGGTCAGGGGGCGGCCGGTGATCACCAGGCCGGCCTTGGCAGGCTGAGGGCCGGACAACCTGACAATTCCCACACCGCGCCGGCCCTGGGCCGTGGCGCTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                               <  7:336160/148‑1 (MQ=255)                                              cATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                         <  5:291004/149‑1 (MQ=255)                                                tGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCct                                                                                                       <  6:256493/149‑1 (MQ=255)                                                   cggggggTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGTCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                                    <  1:179074/146‑1 (MQ=255)                                                      ggggTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑a                                                                                              <  4:251510/149‑1 (MQ=255)                                                            aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACaaaaa                                                                                        >  8:29799/1‑149 (MQ=255)                                                               gTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGa                                                                                     >  3:91204/1‑149 (MQ=255)                                                                tGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGAc                                                                                    <  3:235232/149‑1 (MQ=255)                                                                tGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGAc                                                                                    <  4:91204/149‑1 (MQ=255)                                                                   ggCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAACACGACAg                                                                                  >  6:193751/1‑148 (MQ=255)                                                                   ggCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGa                                                                                     >  5:283917/1‑145 (MQ=255)                                                                      cGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAg                                                                              >  6:187129/1‑149 (MQ=255)                                                                          gATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCaa                                                                           >  2:213798/1‑148 (MQ=255)                                                                              aaCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACg                                                                      >  7:145904/1‑149 (MQ=255)                                                                               aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGc                                                                     <  7:313341/149‑1 (MQ=255)                                                                                 aGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCTGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGccc                                                                   <  7:29799/149‑1 (MQ=255)                                                                                   gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  5:209929/1‑149 (MQ=255)                                                                                   gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCCGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  7:164425/1‑149 (MQ=255)                                                                                     cggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCAc                                                                >  2:206008/1‑148 (MQ=255)                                                                                         cTTGGCAGCCAGAGGGCCGGACAACCTGAC‑ATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑ggg                                                         >  8:266890/1‑149 (MQ=255)                                                                                            ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGc                                                       >  5:178713/1‑149 (MQ=255)                                                                                            ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAAG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACGAg‑ggggg                                                       >  8:348659/1‑148 (MQ=255)                                                                                               aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑tt                                                   >  8:226157/1‑149 (MQ=255)                                                                                                 ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑tttt                                                 <  7:209605/149‑1 (MQ=255)                                                                                                 ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCG‑tttt                                                 <  6:316892/149‑1 (MQ=255)                                                                                                  cAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑ttttt                                                <  5:125251/149‑1 (MQ=255)                                                                                                    gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTa                                              >  7:109553/1‑149 (MQ=255)                                                                                                    gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTa                                              <  7:46721/149‑1 (MQ=255)                                                                                                    gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTa                                              <  6:84587/149‑1 (MQ=255)                                                                                                    gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTa                                              <  6:46067/149‑1 (MQ=255)                                                                                                          cGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTcaca                                        <  7:280350/149‑1 (MQ=255)                                                                                                            gACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATAGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAgg                                      <  6:6681/149‑1 (MQ=255)                                                                                                             aCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCAACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GCGGCG‑TTTTTTATTCACAGGc                                     >  3:119584/1‑149 (MQ=255)                                                                                                              cAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑ttttt                                                <  1:325452/137‑1 (MQ=255)                                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAACACGACAGATAGCACAAACGCCCCACTAG‑GGGGCG‑TGTTTTATTACCAGGCTAGc                                 >  2:290469/1‑149 (MQ=255)                                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGc                                 >  3:47450/1‑149 (MQ=255)                                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGc                                 <  3:15233/149‑1 (MQ=255)                                                                                                                    gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTa                              >  3:248352/1‑149 (MQ=255)                                                                                                                     aCAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGAt                                                                                <  2:351359/100‑1 (MQ=255)                                                                                                                      cAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTTATTCACAGGCTAGCGTa                              <  8:129526/148‑1 (MQ=255)                                                                                                                      cAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTAt                             <  2:204801/148‑1 (MQ=255)                                                                                                                       aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTAt                             >  3:274904/1‑147 (MQ=255)                                                                                                                       aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTa                           <  6:209929/149‑1 (MQ=255)                                                                                                                        aTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAg                          <  7:266890/149‑1 (MQ=255)                                                                                                                         ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGc                         >  3:147566/1‑149 (MQ=255)                                                                                                                           cccACACCGCCCCGGCCCTTGGCCGTGGCGATTGCGGCGCTGGTTTCACGCCCAGTTTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAACACTACAGATAGCAAAACCGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGccc                       <  2:171336/149‑1 (MQ=255)                                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGccc                       <  8:109553/148‑1 (MQ=255)                                                                                                                             cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGt                     >  8:82255/1‑149 (MQ=255)                                                                                                                              ccatcGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGt                     >  3:268713/5‑148 (MQ=255)                                                                                                                              acacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                    <  2:310167/149‑1 (MQ=255)                                                                                                                               cacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGCTTACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCg                      <  1:115959/149‑1 (MQ=255)                                                                                                                                acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGc                                                                     <  3:272389/100‑1 (MQ=255)                                                                                                                                acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGc                                                                     >  4:272389/1‑100 (MQ=255)                                                                                                                                 ccGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                 >  8:333938/1‑149 (MQ=255)                                                                                                                                   gCCCCGGCCCTGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                 >  5:74243/1‑147 (MQ=255)                                                                                                                                   gCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                >  2:260522/1‑148 (MQ=255)                                                                                                                                   gCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAAG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGTGGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                <  1:260522/149‑1 (MQ=255)                                                                                                                                       cGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag           >  6:16589/1‑149 (MQ=255)                                                                                                                                            cTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttt      <  4:147566/149‑1 (MQ=255)                                                                                                                                             cgggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt     <  8:164425/148‑1 (MQ=255)                                                                                                                                               gggggTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGCCAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg   >  7:46073/5‑149 (MQ=255)                                                                                                                                               ggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg   <  8:145904/149‑1 (MQ=255)                                                                                                                                                 ccGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg  >  4:203128/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCG‑TTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGG  >  NC_002947/6199‑6490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑