Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,330	G→C	5.6%	A66A (GCC→GCG)	pyrC ←	dihydroorotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,245,330	0	G	C	5.6%	227.8 / 3.6	71	A66A (GCC→GCG)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base C (2/2);  total (37/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 = 2.93e-01

TGGTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGC  >  NC_002947/1245182‑1245477                                                                                                                                                     |                                                                                                                                                      tGGTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTgccg                                                                                                                                                       >  6:435872/1‑149 (MQ=255)           ggAATCGGAGTTGGTGGTGGCGCCGGCCGGGTCCAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGGGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcggg                                                                                                                                             <  1:110785/149‑1 (MQ=255)              aTCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGTGGCTGCCCGCCCGGCCGGCAg                                                                                                                                          <  6:353786/149‑1 (MQ=255)               tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGTCGGGCAGc                                                                                                                                         <  7:255317/149‑1 (MQ=255)                 ggAGGTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCa                                                                                                                                       <  4:192307/149‑1 (MQ=255)                  gAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAg                                                                                                                                      <  3:35331/149‑1 (MQ=255)                  gAGTTGGTGGGGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCGGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGccgc                                                                                                                                      <  8:152494/149‑3 (MQ=255)                     ttGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGaa                                                                                                                                    <  2:149295/148‑1 (MQ=255)                        gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                                <  7:342324/149‑1 (MQ=255)                         tggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                               <  4:209261/149‑1 (MQ=255)                         tggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                               <  1:322889/149‑1 (MQ=255)                             ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCAc                                                                                                                           >  1:8987/1‑149 (MQ=255)                              gcgcCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACg                                                                                                                          >  6:1862/1‑149 (MQ=255)                                gccggccgGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACg‑g                                                                                                                        >  6:46208/1‑148 (MQ=255)                                 ccggccggGTACAACTTGGCGGTGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACg‑g                                                                                                                        >  4:42500/1‑147 (MQ=255)                                   ggccggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAgg                                                                                                                    <  5:427627/149‑1 (MQ=255)                                        ggTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCg                                                                                                               <  3:42277/149‑1 (MQ=255)                                        ggTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCg                                                                                                               <  2:8987/149‑1 (MQ=255)                                             aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTc                                                                                                          >  4:128183/1‑149 (MQ=255)                                             aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTc                                                                                                          <  5:338057/149‑1 (MQ=255)                                              agttGGCTGCGGACACGAATCCGCTGGGCTTGGCCGGTCGGACGTCCTCGGGGCTGGGGCGGTCGGTGTGTTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCTGGCGGGCAGCCAGAATGAGCTCACG‑GTAGGCGCCGGCCTCg                                                                                                         <  8:73175/147‑1 (MQ=255)                                              aCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCg                                                                                                         <  5:435872/149‑1 (MQ=255)                                                ttGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATGAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGt                                                                                                       <  3:187080/149‑1 (MQ=255)                                                        gTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCCCACG‑GTAGGCGCCGGCCTCGGTGACATTgc                                                                                               <  4:323189/149‑1 (MQ=255)                                                         tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTgcg                                                                                              >  5:235017/1‑149 (MQ=255)                                                         tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTgcg                                                                                              <  3:42500/149‑1 (MQ=255)                                                           cacGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCa                                                                                            >  3:68920/1‑149 (MQ=255)                                                            acGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTgc                                                                                               <  5:163006/145‑1 (MQ=255)                                                            acGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTgc                                                                                               >  6:163006/1‑145 (MQ=255)                                                            acGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCAc                                                                                           >  4:60915/1‑149 (MQ=255)                                                              gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                         >  7:293283/1‑149 (MQ=255)                                                              gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                         >  1:30135/1‑149 (MQ=255)                                                              gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                         >  2:73753/1‑149 (MQ=255)                                                              gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATT‑CGCACggg                                                                                        >  2:364328/1‑149 (MQ=255)                                                                  ccGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                     >  2:107914/1‑149 (MQ=255)                                                                   cGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGggcggc                                                                                    >  5:136103/1‑149 (MQ=255)                                                                   cGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                     >  2:134379/1‑148 (MQ=255)                                                                   cGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                     >  2:139711/1‑148 (MQ=255)                                                                     cTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCAc                                                                                  <  6:235017/149‑1 (MQ=255)                                                                      tGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCAc                                                                                  >  3:75965/1‑148 (MQ=255)                                                                        gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAg                                                                               <  6:136103/149‑1 (MQ=255)                                                                             ggCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                           <  1:73753/148‑1 (MQ=255)                                                                              gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                         <  3:128183/149‑1 (MQ=255)                                                                              cccGCGCGGACGTCCTCGGGGCTCGTGCGGTCCGTGAGGTACAGCAGGATCAGCCGCTCGAAGCGGCTGCCGGCCGGCCGCGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGCTGACATTGCGCACTGGCGGCACGAGGTTcgg                                                                         >  4:300426/2‑149 (MQ=255)                                                                                    cgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑atgat                                                                  >  2:157178/1‑149 (MQ=255)                                                                                      gACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑atgatga                                                                <  5:397604/149‑1 (MQ=255)                                                                                         gTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATgg                                                             >  1:69223/1‑149 (MQ=255)                                                                                           ccTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGTCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCa                                                           >  1:115334/1‑149 (MQ=255)                                                                                            cacGGGGCTGGTTCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCAGGCCGGGCGGGCAGCCAGAATGAGCTCACG‑GTAGGCGCCGGCCTCGGTGACAATGCGCACGGGCGCCACGAGGTTCGGC‑ATGATGATGGCAc                                                          >  2:414285/3‑149 (MQ=255)                                                                                             tgggggCTGGTGCGGTCGGTGAGGTACAGGACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGAGCTCTCT‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCAAGAGGTTGGGC‑ATGATGATGGCACg                                                         >  7:280248/3‑149 (MQ=255)                                                                                              cGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACgg                                                        >  3:118088/1‑149 (MQ=255)                                                                                               ggggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACggg                                                       <  4:118088/149‑1 (MQ=255)                                                                                                  gCTGGTGCGGTCGGTCAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGa                                                    <  8:103417/149‑1 (MQ=255)                                                                                                          ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTAcgc                                            <  4:68920/149‑1 (MQ=255)                                                                                                           gtaggtGAGGTACAGCACCATAAGCGGCTCGAATCGTATCCCGGCAGGGCGGGCAGACAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGAGACATAGCGCACGGGCGGCAGGAGTTGGGGC‑ATGATGATGGCACGGGCGAATGTAcgcg                                           >  2:306372/4‑149 (MQ=255)                                                                                                                      aCAGGACGATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCa                                <  3:307739/149‑1 (MQ=255)                                                                                                                      aCAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                 >  6:210202/1‑148 (MQ=255)                                                                                                                        aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAAc                              <  1:134379/149‑1 (MQ=255)                                                                                                                        aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAAc                              <  1:139711/149‑1 (MQ=255)                                                                                                                          cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGAATAGGCGCCGGCCTCGGTGCCATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCCCGTCGCCCACg                             >  1:125465/1‑149 (MQ=255)                                                                                                                          cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                            >  2:92107/1‑149 (MQ=255)                                                                                                                          cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                            >  5:186158/1‑149 (MQ=255)                                                                                                                             cATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtgt                         >  5:346613/1‑149 (MQ=255)                                                                                                                                  gCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgggg                     >  8:391607/1‑148 (MQ=255)                                                                                                                                   cGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGATGTTGGGCTATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGc                    <  1:107914/149‑1 (MQ=255)                                                                                                                                    ggCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAg                  <  4:75965/149‑1 (MQ=255)                                                                                                                                       tCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGAc               >  6:190424/1‑149 (MQ=255)                                                                                                                                         gAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg             >  6:284347/1‑149 (MQ=255)                                                                                                                                         gAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg             <  3:60915/149‑1 (MQ=255)                                                                                                                                               gCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCTCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGTCACCAt       >  6:196869/1‑149 (MQ=255)                                                                                                                                               gCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCAt       >  1:40004/1‑149 (MQ=255)                                                                                                                                                 tGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcg     >  6:123340/1‑149 (MQ=255)                                                                                                                                                  gccggccgGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcgc    >  4:49269/1‑149 (MQ=255)                                                                                                                                                   ccggccggGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcg     <  1:92107/147‑1 (MQ=255)                                                                                                                                                    cggccggGCGGTCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGTCAGGACGGCACCATCgcgc  <  1:414285/149‑1 (MQ=255)                                                                                                                                                    cggccggGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCgcgc  <  5:1862/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TGGTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGC  >  NC_002947/1245182‑1245477

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