Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,334	G→C	5.7%	P65R (CCG→CGG)	pyrC ←	dihydroorotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,245,334	0	G	C	5.7%	223.0 / 3.7	70	P65R (CCG→CGG)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base C (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245192‑1245481                                                                                                                                               |                                                                                                                                                      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                                                                                                              <  5:338057/149‑1 (MQ=255)                                   aaCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTc                                                                                                              >  4:128183/1‑149 (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                                                                                             >  2:73753/1‑149 (MQ=255)                                                    gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                             >  7:293283/1‑149 (MQ=255)                                                    gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                             >  1:30135/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:139711/1‑148 (MQ=255)                                                         cGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                         >  2:134379/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:139711/149‑1 (MQ=255)                                                                                                              aGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAAc                                  <  1:134379/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                 <  3:60915/149‑1 (MQ=255)                                                                                                                               gAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg                 >  6:284347/1‑149 (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)                                                                                                                                         cccgcccgGCCGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcg         <  2:40004/144‑1 (MQ=255)                                                                                                                                          cggccggGCGGTCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGTCAGGACGGCACCATCgcgc      <  1:414285/149‑1 (MQ=255)                                                                                                                                          cggccggGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCgcgc      <  5:1862/149‑1 (MQ=255)                                                                                                                                              cgggcgggCAGCCAGATTGCGCTCACG‑GTAGGCGCCGGCCTCGGTTACATTGCGCACGGGCTGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  <  8:280248/149‑1 (MQ=255)                                                                                                                                              cgggcgggCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  >  4:405304/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      GGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACG‑GTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGC‑ATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245192‑1245481

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