Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,334	G→C	6.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	6.7%	191.1 / 4.3	60	P65R (CCG→CGG)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base C (2/2);  total (26/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.63e-01

GCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGC  >  NC_002947/1245189‑1245477                                                                                                                                                  |                                                                                                                                                gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggc                                                                                                                                              >  3:166948/1‑149 (MQ=255)   cGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcgg                                                                                                                                            <  8:194879/149‑1 (MQ=255)   cGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcgg                                                                                                                                            >  1:360196/1‑149 (MQ=255)        tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                       <  7:110488/149‑1 (MQ=255)        tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                       <  1:35833/149‑1 (MQ=255)           gAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGccgc                                                                                                                                    <  6:102824/149‑3 (MQ=255)              ttGGTGGTGGAGCCGGCCGTGTAGAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                 <  8:94002/149‑1 (MQ=255)                ggtggtggGGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACGGGTGGCGTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGGTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATg                                                                                                                                <  6:243152/148‑1 (MQ=255)                 gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                              <  2:360196/149‑1 (MQ=255)                 gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                               <  4:166948/148‑1 (MQ=255)                        cgcCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGTGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACgg                                                                                                                       <  1:95134/149‑1 (MQ=255)                        cgcCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACgg                                                                                                                       <  3:284230/149‑1 (MQ=255)                             gccggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGgc                                                                                                                  <  8:260206/149‑1 (MQ=255)                             gccggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGgc                                                                                                                  <  1:243277/149‑1 (MQ=255)                               cggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGgcgc                                                                                                                >  2:358958/1‑149 (MQ=255)                                        cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACgg                                                                                                                       >  7:85500/1‑133 (MQ=255)                                        cTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACgg                                                                                                                       <  8:85500/133‑1 (MQ=255)                                               gcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACAtt                                                                                                <  4:142466/149‑1 (MQ=255)                                                cgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTg                                                                                               <  8:120482/149‑1 (MQ=255)                                                     acGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCAc                                                                                          >  2:313934/1‑149 (MQ=255)                                                      cGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACg                                                                                         >  6:15629/1‑149 (MQ=255)                                                       gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCg                                                                                                              <  3:197446/127‑1 (MQ=255)                                                       gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCg                                                                                                              >  4:197446/1‑127 (MQ=255)                                                        aaaCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAg                                                                                                                                        <  5:347485/100‑1 (MQ=255)                                                        aaaCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAg                                                                                                                                        >  6:347485/1‑100 (MQ=255)                                                        aaaCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACggg                                                                                       >  2:108414/1‑149 (MQ=255)                                                          aCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcg                                                                                     >  6:122714/1‑149 (MQ=255)                                                               tGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCAc                                                                                 >  1:169812/1‑148 (MQ=255)                                                                       gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                         <  1:288777/148‑1 (MQ=255)                                                                        ccGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGc                                                                       <  2:169812/149‑1 (MQ=255)                                                                         cgcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCa                                                                      <  1:358958/149‑1 (MQ=255)                                                                                 cGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatgatg                                                              <  6:328476/149‑1 (MQ=255)                                                                                  gTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatgatg                                                              <  7:340738/148‑1 (MQ=255)                                                                                  gTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatgatg                                                              >  8:340738/1‑148 (MQ=255)                                                                                    ccTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                             >  4:301441/1‑83 (MQ=255)                                                                                    ccTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                             <  3:301441/83‑1 (MQ=255)                                                                                     cTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgggt                                                                       >  3:168814/1‑135 (MQ=255)                                                                                     cTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgggt                                                                       <  4:168814/136‑2 (MQ=255)                                                                                               gTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGt                                                <  5:45358/149‑1 (MQ=255)                                                                                               gTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGt                                                <  1:313934/149‑1 (MQ=255)                                                                                                    gtcggtGAGGTACAGCACCATTAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                        >  2:229338/1‑120 (MQ=255)                                                                                                    gtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                        <  1:229338/120‑1 (MQ=255)                                                                                                     tcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACgcgc                                          >  5:308746/1‑149 (MQ=255)                                                                                                        gtGAGGTACAGCACCATCAGCGGCTCGAAGTGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCAc                                       <  7:75651/149‑1 (MQ=255)                                                                                                               aCAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAg                                                  <  8:55867/131‑1 (MQ=255)                                                                                                               aCAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATCGCACGGGCGAAg                                                  >  7:55867/1‑131 (MQ=255)                                                                                                                   cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                            >  1:258283/1‑149 (MQ=255)                                                                                                                    aCCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgt                           >  1:163167/1‑149 (MQ=255)                                                                                                                       aTCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                        <  5:15629/149‑1 (MQ=255)                                                                                                                       aTCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                        <  5:122714/149‑1 (MQ=255)                                                                                                                        tCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                       >  8:42767/1‑149 (MQ=255)                                                                                                                              gCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAgg                 <  2:163167/149‑1 (MQ=255)                                                                                                                                  gAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg             >  7:169063/1‑149 (MQ=255)                                                                                                                                   aaGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGc            <  8:169063/149‑1 (MQ=255)                                                                                                                                   aaGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGG‑CTCGGTGACCTTGCGCACGGGCGGCCCCACGTTGGGCATGATGATGGCACGGGCGAAGGTACCCGCCACGTCGCCAACGGTGTGGGGCAGGACGGcc           >  5:105080/1‑148 (MQ=255)                                                                                                                                    aGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCa           <  6:220098/149‑1 (MQ=255)                                                                                                                                    aGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTAGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCa           >  4:53327/1‑149 (MQ=255)                                                                                                                                     ctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCAc          >  2:221731/3‑149 (MQ=255)                                                                                                                                      tggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                   <  7:236029/106‑1 (MQ=255)                                                                                                                                      tggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                   >  8:236029/2‑107 (MQ=255)                                                                                                                                        gCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCAt       >  7:60213/1‑149 (MQ=255)                                                                                                                                          tGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcg     >  5:226529/1‑149 (MQ=255)                                                                                                                                          tGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcc     >  4:337690/1‑148 (MQ=255)                                                                                                                                             cggccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCgcgc  <  2:258283/149‑1 (MQ=255)                                                                                                                                               gccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATTGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcgcg   >  8:79573/1‑146 (MQ=255)                                                                                                                                                  |                                                                                                                                                GCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGC  >  NC_002947/1245189‑1245477

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