Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,330	G→C	9.3%	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	9.3%	226.2 / 8.3	72	A66A (GCC→GCG)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base C (3/4);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCG  >  NC_002947/1245191‑1245476                                                                                                                                              |                                                                                                                                                   cGGAATCGGAGTTGGAGGTGGCGCCGGCCGGGTACAAGCTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCcgggcg                                                                                                                                            >  1:69313/1‑149 (MQ=255)    aaTCGGAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCa                                                                                                                                        <  4:252188/149‑1 (MQ=255)      tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCAGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CCGCCCGGCCGGCAGc                                                                                                                                      <  5:389607/149‑1 (MQ=255)        ggAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCa                                                                                                                                    <  7:230153/149‑1 (MQ=255)         gAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGc                                                                                                                                      >  8:230543/1‑146 (MQ=255)           gTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGaa                                                                                                                                 <  1:385293/149‑1 (MQ=255)            ttGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CCGCCCGGCCGGCAGccgcttc                                                                                                                                <  7:202997/149‑6 (MQ=255)               gtggtggCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                             >  7:166146/1‑149 (MQ=255)               gtggtggCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTAGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                             >  2:144277/1‑149 (MQ=255)               gtggtggAGCCGGCCGTTTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGTACGTCCTCGGGGCTGGTGCGGTCTGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                             <  2:188252/149‑1 (MQ=255)                tggtggCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATgcgc                                                                                                                            <  5:100221/149‑1 (MQ=255)                   tggCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGc                                                                                                                                      >  8:280581/1‑136 (MQ=255)                   tggCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGc                                                                                                                                      <  7:280581/136‑1 (MQ=255)                     gcgcCGGCCGGGTACAA‑CTTGGCGTCGTACACGAAACCGCTGGCTTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGTCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACg                                                                                                                       <  3:15769/149‑1 (MQ=255)                              gggTACAA‑CTTCGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCTGTGATGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGcc                                                                                                              <  8:257395/149‑1 (MQ=255)                               ggTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTTAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCg                                                                                                             <  8:177913/149‑1 (MQ=255)                               ggTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCcgggcgg                                                                                                                                           >  1:274989/1‑119 (MQ=255)                               ggTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCcgggcgg                                                                                                                                           <  2:274989/119‑1 (MQ=255)                                 tACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGc                                                                                                           >  5:119551/1‑149 (MQ=255)                                   cAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCt                                                                                                         >  8:239325/1‑149 (MQ=255)                                    aa‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                                        <  2:69313/149‑1 (MQ=255)                                     a‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCg                                                                                                       >  7:330632/1‑149 (MQ=255)                                     a‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAAAGCGCTCACGGTAGGCGCCGGCCTCg                                                                                                       >  7:135356/1‑149 (MQ=255)                                         tGGCGGCGTACACGATAGCGCTGGCCTTGGAGGCGCGTACGTCGTCGGGGCTGTGGCGGGCGGGGAGGTACTGCAGCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCATCGATAATGCGCTCAGTGTAGGCGCCGGCCTCGGt                                                                                                     <  8:327667/148‑1 (MQ=255)                                          ggcggcGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGt                                                                                                     <  5:174225/147‑1 (MQ=255)                                           gcggcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGGCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGa                                                                                                   <  2:65219/148‑1 (MQ=255)                                                gTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                             >  4:391251/1‑149 (MQ=255)                                                 tACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgcg                                                                                            >  8:164745/1‑149 (MQ=255)                                                      gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                       <  3:65578/149‑1 (MQ=255)                                                      gAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACgg                                                                                       <  7:173306/149‑1 (MQ=255)                                                        aaCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggc                                                                                     <  1:144277/149‑1 (MQ=255)                                                         aCCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcg                                                                                    >  4:185871/1‑149 (MQ=255)                                                          ccGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                   <  7:337913/149‑1 (MQ=255)                                                          ccGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcgg                                                                                   <  7:230543/149‑1 (MQ=255)                                                             cTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGACTGCGCGCCCGGTAGGCGCCGGCCTCGGTGCCAGTGGGCAGGggcggccc                                                                                >  8:258218/1‑147 (MQ=255)                                                                gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAg                                                                             <  1:204710/149‑1 (MQ=255)                                                                   ttGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCCGTGACATTGCGCACGGGCGGCACGacgat                                                                          >  8:125926/1‑145 (MQ=255)                                                                      gCCGCGCGGACGTGCTCTGGGCTGGTTCGGTCGGTGAGGTACAGGAGCATGAGCGGCTCGAAGCGGCTGT‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                        <  2:241857/148‑1 (MQ=255)                                                                      gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                       <  8:83835/149‑1 (MQ=255)                                                                        cgcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGCCATTGCGCACGGGCGGCACGAGGTTGGGCa                                                                     >  6:27952/1‑149 (MQ=255)                                                                         gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCat                                                                    <  5:307477/149‑1 (MQ=255)                                                                           gcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                  <  2:183869/149‑1 (MQ=255)                                                                               aCGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggc                                                                                     >  6:23603/1‑126 (MQ=255)                                                                               aCGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggc                                                                                     <  5:23603/126‑1 (MQ=255)                                                                                 gTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATgg                                                            >  6:119446/1‑149 (MQ=255)                                                                                    cTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCAc                                                         >  5:26998/1‑149 (MQ=255)                                                                                    cTCGGGGCTGGTGCGGGCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGGACGAGGTTGGGCATGATGATGGCAc                                                         <  7:199217/149‑1 (MQ=255)                                                                                        gggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                         <  5:115645/129‑1 (MQ=255)                                                                                        gggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTg                                                                         >  6:115645/1‑129 (MQ=255)                                                                                                  ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                           <  7:164745/149‑1 (MQ=255)                                                                                                        tGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCGCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatga                                                               >  1:261841/1‑124 (MQ=255)                                                                                                        tGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatga                                                               <  2:261841/123‑1 (MQ=255)                                                                                                          aGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCCCGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTc                                   <  8:166146/149‑1 (MQ=255)                                                                                                           ggTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcggc                                                                                  >  6:124243/1‑101 (MQ=255)                                                                                                           ggTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGggcggc                                                                                  <  5:124243/101‑1 (MQ=255)                                                                                                            gTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAACCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGc                                 >  2:166619/1‑149 (MQ=255)                                                                                                              aCAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGTCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                >  1:279694/1‑148 (MQ=255)                                                                                                                    ccATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtg                         <  6:26998/149‑1 (MQ=255)                                                                                                                     cATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtgt                        >  6:140123/1‑149 (MQ=255)                                                                                                                     cATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtgt                        >  6:280243/1‑149 (MQ=255)                                                                                                                       tCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                      >  5:81067/1‑149 (MQ=255)                                                                                                                           cGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgggg                    >  7:213111/1‑147 (MQ=255)                                                                                                                             gCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                        <  1:208297/93‑1 (MQ=255)                                                                                                                             gCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTgg                                                                        >  2:208297/1‑93 (MQ=255)                                                                                                                                cGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGTTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGAACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACg             >  5:199629/1‑149 (MQ=255)                                                                                                                                 attctggCTGC‑CCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg            >  5:53193/6‑149 (MQ=255)                                                                                                                                   aGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCa          <  7:239325/149‑1 (MQ=255)                                                                                                                                    gCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGccc         >  8:147967/1‑147 (MQ=255)                                                                                                                                     tggCTGC‑CCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTggg                     <  6:66029/135‑1 (MQ=255)                                                                                                                                     tggCTGC‑CCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTggg                     >  5:66029/2‑136 (MQ=255)                                                                                                                                     tggCTGC‑CCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                      <  3:375242/134‑1 (MQ=255)                                                                                                                                     tggCTGC‑CCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                      >  4:375242/2‑135 (MQ=255)                                                                                                                                             cggccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATcgcg  >  6:372508/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                   CGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAA‑CTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGC‑CGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCG  >  NC_002947/1245191‑1245476

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