Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,735,476	C→A	6.0%	G350C (GGC→TGC)	csdA ←	cysteine desulfurase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,735,476	0	C	A	6.0%	211.1 / 3.5	67	G350C (GGC→TGC)	csdA	cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

AGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAGGCTGGCTTCAT  >  NC_002947/1735328‑1735622                                                                                                                                                     |                                                                                                                                                    aGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGcc                                                                                                                                                     <  8:257040/149‑1 (MQ=255)      cAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGt                                                                                                                                                 >  7:98753/1‑148 (MQ=255)        ggCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtca                                                                                                                                              <  8:98753/149‑1 (MQ=255)        ggCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtca                                                                                                                                              <  2:237626/149‑1 (MQ=255)         gCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcag                                                                                                                                             >  8:252925/1‑149 (MQ=255)          ccTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagt                                                                                                                                            <  7:43966/149‑1 (MQ=255)           cTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagt                                                                                                                                            >  2:111615/1‑148 (MQ=255)              gggccAGCGCGTCGAAGAAACGCTGGACGTCATCGCTGGCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTCGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagtca                                                                                                                                          <  4:236128/144‑1 (MQ=255)                  cAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAAt                                                                                                                                    <  6:148450/149‑1 (MQ=255)                   aGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATg                                                                                                                                   >  4:282294/1‑149 (MQ=255)                    gcgcGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGa                                                                                                                                  <  4:39787/149‑1 (MQ=255)                           gaagaaACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCCCGGCAATGCCCTGTTCAGTCAGCAAATGAGCAatat                                                                                                                           >  3:246726/1‑149 (MQ=255)                            aagaaACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCACTGCCCTGATCAGTCAGCAAATGAGCAAtctc                                                                                                                          >  8:255524/1‑149 (MQ=255)                            aagaaACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACACACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATc                                                                                                                          >  1:65200/1‑149 (MQ=255)                                 aCGCTGCACGTCATCGCTGTCGTTGTACAACCCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCATTACGCAAATGAGAAATATCGGCGt                                                                                                                     >  7:259639/1‑149 (MQ=255)                                  cGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGAATCTCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGtt                                                                                                                    <  5:72452/149‑1 (MQ=255)                                   gCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGtt                                                                                                                    <  4:100524/148‑1 (MQ=255)                                   gCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGtt                                                                                                                    >  3:100524/1‑148 (MQ=255)                                    cTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTtgt                                                                                                                  >  7:51711/1‑149 (MQ=255)                                    cTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTtgt                                                                                                                  <  7:252925/149‑1 (MQ=255)                                       cACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCa                                                                                                               <  2:225551/149‑1 (MQ=255)                                       cACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCa                                                                                                               >  8:268862/1‑149 (MQ=255)                                         cGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACACACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACg                                                                                                             <  2:65200/149‑1 (MQ=255)                                                    tCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGa                                                                                                   <  3:291099/148‑1 (MQ=255)                                                      gTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGa                                                                                                <  8:91487/149‑1 (MQ=255)                                                       ttGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGa                                                                                                >  2:212629/1‑148 (MQ=255)                                                          tACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGCCCGGCCCGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACt                                                                                            >  2:269483/1‑149 (MQ=255)                                                                cccAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAg                                                                                      <  1:111615/149‑1 (MQ=255)                                                                      gATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGt                                                                                >  6:291715/1‑149 (MQ=255)                                                                           ccGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCg                                                                           >  7:214235/1‑149 (MQ=255)                                                                           ccGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCACTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCg                                                                           >  6:301009/1‑149 (MQ=255)                                                                           ccGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCACTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCg                                                                           >  6:296155/1‑149 (MQ=255)                                                                                  gcgcCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGTGCAGTGATGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGcc                                                                    >  6:307949/1‑149 (MQ=255)                                                                                  gcgcCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTgcgc                                                                     <  8:174335/148‑1 (MQ=255)                                                                                  gcgcCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTgcgc                                                                     <  8:179062/148‑1 (MQ=255)                                                                                    gcCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGccc                                                                   >  3:223220/1‑148 (MQ=255)                                                                                    gcCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCa                                                                  <  7:268862/149‑1 (MQ=255)                                                                                      ccTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCaa                                                                 >  6:112823/1‑148 (MQ=255)                                                                                      ccTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCaa                                                                 >  6:117938/1‑148 (MQ=255)                                                                                          cagccccagcccCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACAccc                                                            >  5:107101/1‑149 (MQ=255)                                                                                             ccccagccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCa                                                         >  6:29271/1‑149 (MQ=255)                                                                                               ccagccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGt                                                                                >  2:27604/1‑124 (MQ=255)                                                                                               ccagccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGt                                                                                <  1:27604/124‑1 (MQ=255)                                                                                               ccagccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACg                                                       >  1:16190/1‑149 (MQ=255)                                                                                               ccagccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACg                                                       >  5:211508/1‑149 (MQ=255)                                                                                                  gccccTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGccc                                                    <  3:282294/149‑1 (MQ=255)                                                                                                      cTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑cgc                                               <  7:255524/149‑1 (MQ=255)                                                                                                      cTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑cgc                                               >  6:97254/1‑149 (MQ=255)                                                                                                       ttCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑cgcg                                              >  4:21043/1‑149 (MQ=255)                                                                                                          agcagTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTc                                           <  5:291715/149‑1 (MQ=255)                                                                                                            cagTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGc                                         >  2:5659/1‑149 (MQ=255)                                                                                                                        gcAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGt                                                                                <  5:75361/99‑1 (MQ=255)                                                                                                                        gcAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCAAGCCCTCGATGACGAAACTGGCCAGGGCCGt                                                                                >  6:75361/1‑99 (MQ=255)                                                                                                                           gtggtgACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACgcagca                          <  5:117938/149‑1 (MQ=255)                                                                                                                           gtggtgACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACgcagca                          <  5:112823/149‑1 (MQ=255)                                                                                                                           gtggtgACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACgcagca                          <  3:21043/149‑1 (MQ=255)                                                                                                                                 aCCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGc                     <  1:5659/148‑1 (MQ=255)                                                                                                                                  ccGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGTGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTg                   <  8:51711/149‑1 (MQ=255)                                                                                                                                  ccGGCACGCACGGCAATGCACTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTg                   <  5:296155/149‑1 (MQ=255)                                                                                                                                  ccGGCACGCACGGCAATGCACTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTg                   <  5:301009/149‑1 (MQ=255)                                                                                                                                   cGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTgg                  >  7:189468/1‑149 (MQ=255)                                                                                                                                   cGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTgg                  >  5:2868/1‑149 (MQ=255)                                                                                                                                   cGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGAACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTgg                  >  1:252000/1‑149 (MQ=255)                                                                                                                                        cgcacgGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCTCCGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCa              <  1:212629/149‑1 (MQ=255)                                                                                                                                         gcacgGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAgg            >  2:267346/1‑149 (MQ=255)                                                                                                                                           acgGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAggct          <  1:4956/149‑1 (MQ=255)                                                                                                                                             ggCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACG‑CCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAggctggc       >  3:74818/1‑149 (MQ=255)                                                                                                                                                  tGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAGGCTGGCTTCa   <  3:194621/149‑1 (MQ=255)                                                                                                                                                    cccTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAGGCTGGCTTCAt  <  6:167103/148‑1 (MQ=255)                                                                                                                                                    cccTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAGGCTGGCTTCAt  <  6:2868/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGGTGCGCCCAACACCCGCACGCCCT‑CGCGGTCGCCCAGGCCACGCAGCAGGTGCTGGTGCAGGCTGGCTTCAT  >  NC_002947/1735328‑1735622

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