Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,338	G→C	10.5%	R64G (CGC→GGC)	pyrC ←	dihydroorotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,245,338	0	G	C	10.5%	157.3 / 8.7	57	R64G (CGC→GGC)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base C (2/4);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

CCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245190‑1245481                                                                                                                                                      |                                                                                                                                                 ccGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcg                                                                                                                                                  <  4:316240/149‑1 (MQ=255)       tCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGTCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAg                                                                                                                                             >  1:272751/1‑149 (MQ=255)       tCGGAGTTGGTGGGGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGc                                                                                                                                            <  6:330082/149‑1 (MQ=255)        cGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CGGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCCCGGCCGGCAGcc                                                                                                                                           <  1:318029/149‑1 (MQ=255)           aGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcggg                                                                                                                                                >  1:317661/1‑141 (MQ=255)           aGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcggg                                                                                                                                                <  2:317661/141‑1 (MQ=255)            gTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGaa                                                                                                                                       <  5:60191/149‑1 (MQ=255)                gtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                                   >  2:132832/1‑149 (MQ=255)                     ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCgggcgggc                                                                                                                                               >  5:109428/1‑132 (MQ=255)                     ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCgggcgggc                                                                                                                                               <  6:109428/132‑1 (MQ=255)                     ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                            <  4:296795/135‑1 (MQ=255)                     ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                            >  3:296795/1‑135 (MQ=255)                      gcgcCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACg                                                                                                                             >  5:114607/1‑149 (MQ=255)                           ggccggGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAgg                                                                                                                        >  3:311144/1‑149 (MQ=255)                                   aCAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                          <  4:21152/139‑1 (MQ=255)                                   aCAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTa                                                                                                                          >  3:21152/1‑139 (MQ=255)                                     aaCTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                                              >  8:293017/1‑149 (MQ=255)                                     aaCTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                                              >  6:21731/1‑149 (MQ=255)                                          ggcggcGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                      >  7:113109/1‑120 (MQ=255)                                          ggcggcGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                                      <  8:113109/120‑1 (MQ=255)                                             ggcgTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACAt                                                                                                      >  5:167251/1‑149 (MQ=255)                                                 tACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                    >  7:296341/1‑115 (MQ=255)                                                 tACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                                    <  8:296341/115‑1 (MQ=255)                                                  acacGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgcg                                                                                                  >  1:5155/1‑148 (MQ=255)                                                        aaCCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGgg‑c                                                                                          <  1:132832/149‑1 (MQ=255)                                                                 gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAg                                                                                  >  2:323013/1‑149 (MQ=255)                                                                  ccTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAgg                                                                                 <  2:324116/149‑1 (MQ=255)                                                                   cTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGt                                                                                >  7:244215/1‑149 (MQ=255)                                                                    ttGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGtt                                                                               <  7:293017/149‑1 (MQ=255)                                                                       gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTggg                                                                            >  3:177226/1‑149 (MQ=255)                                                                       gCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTgg                                                                             <  5:21731/148‑1 (MQ=255)                                                                          gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                       >  6:235693/1‑120 (MQ=255)                                                                          gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                       <  5:235693/120‑1 (MQ=255)                                                                          gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCat                                                                         <  2:5155/149‑1 (MQ=255)                                                                           cgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatg                                                                        <  2:272751/149‑1 (MQ=255)                                                                           cgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCACCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatg                                                                        <  5:240088/149‑1 (MQ=255)                                                                            gcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTCCGCACGGG‑CGGCACGAGGTTGGGCatga                                                                       <  1:201815/149‑1 (MQ=255)                                                                               gACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGGGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatga                                                                    >  3:28904/1‑149 (MQ=255)                                                                                   tCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTg                                                                              >  1:24013/1‑135 (MQ=255)                                                                                   tCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTg                                                                              <  2:24013/135‑1 (MQ=255)                                                                                    ccTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCa                                                               >  3:24832/1‑149 (MQ=255)                                                                                                tGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                   <  8:244215/149‑1 (MQ=255)                                                                                                tGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                   >  1:345293/1‑149 (MQ=255)                                                                                                  cggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcg                                                 >  5:205077/1‑149 (MQ=255)                                                                                                   ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                <  4:311144/149‑1 (MQ=255)                                                                                                         tGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACg                                          >  8:244628/1‑149 (MQ=255)                                                                                                          gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatg                                                                     >  6:173811/1‑121 (MQ=255)                                                                                                          gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatg                                                                     <  5:173811/121‑1 (MQ=255)                                                                                                                cAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGTCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                     <  1:94162/148‑1 (MQ=255)                                                                                                                   cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatg                                                                     >  3:68788/1‑112 (MQ=255)                                                                                                                   cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatg                                                                     <  4:68788/112‑1 (MQ=255)                                                                                                                   cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACgg                                >  4:260132/1‑149 (MQ=255)                                                                                                                    aCCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCAGGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgt                               <  6:114607/149‑1 (MQ=255)                                                                                                                     ccATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtg                              <  4:24832/149‑1 (MQ=255)                                                                                                                        tCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                           <  6:205077/149‑1 (MQ=255)                                                                                                                         cAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                            <  7:201787/147‑1 (MQ=255)                                                                                                                         cAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                            >  8:201787/1‑147 (MQ=255)                                                                                                                                  attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg                 >  1:232799/6‑149 (MQ=255)                                                                                                                                  attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg                 <  2:232799/144‑1 (MQ=255)                                                                                                                                       ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg                          <  1:322852/135‑1 (MQ=255)                                                                                                                                       ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg                          >  2:322852/1‑135 (MQ=255)                                                                                                                                               gccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAg    <  3:260132/149‑1 (MQ=255)                                                                                                                                                  gggcgggcAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  >  3:105866/1‑148 (MQ=255)                                                                                                                                                  gggcgggcAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAgat  <  4:105866/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 CCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGAT  >  NC_002947/1245190‑1245481

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