Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,409,206	C→G	5.6%	E748Q (GAG→CAG)	acnA‑I ←	aconitate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,409,206	0	C	G	5.6%	239.4 / 4.3	72	E748Q (GAG→CAG)	acnA‑I	aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base G (2/2);  total (38/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.43e-01

GAACCGGTGCCGTACTCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGATGTTG  >  NC_002947/2409064‑2409348                                                                                                                                               |                                                                                                                                               gAACCGGTGCCGTACTCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGc                                                                                                                                           <  4:98252/148‑1 (MQ=255)            gTACTCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGt                                                                                                                               >  1:351740/1‑149 (MQ=255)             tACTCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGtt                                                                                                                              <  4:120822/149‑1 (MQ=255)               cTCCTGGCCGTCAATCACCAGCAGCGGGGTGCCTTCCTGCTTGTAACTCATGGCCGCGTCGTAGATCGACAGCTTGTCGACAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTTCCAGCATCTCGTTGCGAATACGGATGTTg                                                                                                                             <  3:267189/148‑1 (MQ=255)                tCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGc                                                                                                                           >  5:150202/1‑149 (MQ=255)                            aTCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGTATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCcgcg                                                                                                               >  6:151905/1‑149 (MQ=255)                            aTCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCcgcg                                                                                                               <  7:225766/149‑1 (MQ=255)                             tCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCcgcg                                                                                                               >  8:198495/1‑148 (MQ=255)                                 cagcagCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGa                                                                                                                                  >  3:3097/1‑125 (MQ=255)                                 cagcagCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGa                                                                                                                                  <  4:3097/125‑1 (MQ=255)                                        ggggTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTCGCGAAGGTGCCGCGCATCATCACTTc                                                                                                   <  5:336973/149‑1 (MQ=255)                                          ggTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGt                                                                                                 <  7:198495/149‑1 (MQ=255)                                           gTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTg                                                                                                >  1:397230/1‑149 (MQ=255)                                            tGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTgg                                                                                               >  4:38370/1‑149 (MQ=255)                                               cTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTg                                                                                            >  4:286651/1‑149 (MQ=255)                                                ttCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGc                                                                                           >  6:327407/1‑149 (MQ=255)                                                  ccTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGc                                                                                           <  2:397230/147‑1 (MQ=255)                                                         gTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCAcgccgc                                                                                   >  2:9008/1‑148 (MQ=255)                                                         gTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCAcgccgc                                                                                   <  3:38370/148‑1 (MQ=255)                                                            aCGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGaa                                                                                <  6:210729/148‑1 (MQ=255)                                                            aCGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAAc                                                                               >  4:111609/1‑149 (MQ=255)                                                              gCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGt                                                                                                 >  1:143904/1‑129 (MQ=255)                                                              gCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGt                                                                                                 <  2:143904/129‑1 (MQ=255)                                                              gCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCg                                                                             >  3:356509/1‑149 (MQ=255)                                                                aTGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTa                                                                           <  2:351740/149‑1 (MQ=255)                                                                        gtcgtAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGa                                                                   <  5:151905/149‑1 (MQ=255)                                                                         tcgtAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTg                                                                    >  5:252261/1‑147 (MQ=255)                                                                         tcgtAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGaa                                                                  <  3:286651/149‑1 (MQ=255)                                                                           gtAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGt                                                                >  4:80782/1‑149 (MQ=255)                                                                            tAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTc                                                               >  6:208241/1‑149 (MQ=255)                                                                              gATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTcgc                                                             <  6:252261/149‑1 (MQ=255)                                                                                tCGACAGCTTCTCGGCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGt                                                           <  1:378104/149‑1 (MQ=255)                                                                                tCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGaa                                                                  >  5:183525/1‑142 (MQ=255)                                                                                tCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGaa                                                                  <  6:183525/142‑1 (MQ=255)                                                                                tCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGt                                                           <  5:327407/149‑1 (MQ=255)                                                                                tCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGt                                                           <  3:111609/149‑1 (MQ=255)                                                                                   aCAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTcgc                                                             <  3:90055/144‑1 (MQ=255)                                                                                   aCAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTcgc                                                             >  4:90055/1‑144 (MQ=255)                                                                                   aCAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGc                                                        <  6:408608/149‑1 (MQ=255)                                                                                   aCAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGc                                                        >  6:172481/1‑149 (MQ=255)                                                                                        ttCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCca                                                    >  2:250386/1‑148 (MQ=255)                                                                                         tctcGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCaaca                                                  >  3:332727/1‑149 (MQ=255)                                                                                           tcGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGAGAAGGTGCCGCTCATCATCACTTAGTGGTTGCCACGCCGCGAACAGTAGGAGTTGAAGTCGAGTGGCTCCacaca                                                >  8:266772/1‑148 (MQ=255)                                                                                             gCCAGTGGGCATATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGt                                                                                              <  2:400136/101‑1 (MQ=255)                                                                                             gCCAGTGGGCATATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGt                                                                                              >  1:400136/1‑101 (MQ=255)                                                                                             gCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCtt                                              >  6:299519/1‑149 (MQ=255)                                                                                                aGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGtt                                           <  1:9008/149‑1 (MQ=255)                                                                                                 gTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGt                                                                      >  4:123707/1‑121 (MQ=255)                                                                                                 gTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGt                                                                      <  3:123707/121‑1 (MQ=255)                                                                                                     gCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCg                                       >  4:21357/1‑148 (MQ=255)                                                                                                         aTGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGa                                                                        >  3:49658/1‑111 (MQ=255)                                                                                                         aTGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGa                                                                        <  4:49658/111‑1 (MQ=255)                                                                                                         aTGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGt                                                                      <  8:262236/113‑1 (MQ=255)                                                                                                         aTGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGt                                                                      >  7:262236/1‑113 (MQ=255)                                                                                                             aaGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGtt                                                                                             <  6:228081/86‑1 (MQ=255)                                                                                                             aaGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGtt                                                                                             >  5:228081/1‑86 (MQ=255)                                                                                                                   ttGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCccgg                         <  5:208241/148‑1 (MQ=255)                                                                                                                      ccgccTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCccggccgg                     >  5:62128/1‑149 (MQ=255)                                                                                                                            tcttcACCTGGCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATGATCGCTTTGTGGTTGCCACGCGGCGAGCCGTAGGAGTTGAAGTCGCGTGGGTCCAGGCCGTGTTCGCGCAGGTAGAGCCCGGCGGGGCTggc               >  4:164688/1‑147 (MQ=255)                                                                                                                                cACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCg              <  5:172481/146‑1 (MQ=255)                                                                                                                                cACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTc           >  7:356005/1‑149 (MQ=255)                                                                                                                                  ccTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACACCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCtt         >  3:165334/1‑149 (MQ=255)                                                                                                                                   cTGCCAGCATCTGGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTc                                                      <  3:162262/103‑1 (MQ=255)                                                                                                                                   cTGCCAGCATCTGGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTc                                                      <  3:156303/103‑1 (MQ=255)                                                                                                                                   cTGCCAGCATCTGGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTc                                                      >  4:156303/1‑103 (MQ=255)                                                                                                                                   cTGCCAGCATCTGGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTc                                                      >  4:162262/1‑103 (MQ=255)                                                                                                                                    tGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAAACGTAGGAGTTGAAGTCGCGAGGCTCCACACCTTGTTCCCGAAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGa       >  7:337785/1‑149 (MQ=255)                                                                                                                                      ccAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGa       <  3:21357/147‑1 (MQ=255)                                                                                                                                      ccAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGa       <  6:150202/147‑1 (MQ=255)                                                                                                                                      ccAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGa       <  1:349776/147‑1 (MQ=255)                                                                                                                                      ccAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGATg     >  1:406408/1‑149 (MQ=255)                                                                                                                                         gCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGATGTTg  >  7:70968/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                               GAACCGGTGCCGTACTCCTGGCCGGCAATCACCAGCAGCGGGGTGCCTTCCTGCTGGTAACGCATGGCCGCGTCGTAGATCGACAGCTTCTCGCCAGTGGGCACATGCAAGGTGTTGCCGCCTTCTTCACCTGCCAGCATCTCGTTGCGAATACGGATGTTGGCGAAGGTGCCGCGCATCATCACTTCGTGGTTGCCACGCCGCGAACCGTAGGAGTTGAAGTCGCGTGGCTCCACACCTTGTTCCCGCAGGTAGCGCCCGGCCGGGCTGTCGGTCTTGATGTTG  >  NC_002947/2409064‑2409348

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