Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,245,330	G→C	12.5%	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	12.5%	172.9 / 14.9	64	A66A (GCC→GCG)	pyrC	dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (3/5);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-02

GTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGG  >  NC_002947/1245184‑1245466                                                                                                                                                    |                                                                                                                                          gTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTgccggc                                                                                                                                         <  3:87141/149‑1 (MQ=255)    aCGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCACCAGCGGCTCGAAGCGGCTGccggccgg                                                                                                                                      >  6:240088/1‑149 (MQ=255)      gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggc                                                                                                                                    >  2:201815/1‑149 (MQ=255)       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                                                                                                                                <  6:109428/132‑1 (MQ=255)                           ggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCgggcgggc                                                                                                                                >  5:109428/1‑132 (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                                                                                               >  6:21731/1‑149 (MQ=255)                                           aaCTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTc                                                                                               >  8:293017/1‑149 (MQ=255)                                                ggcggcGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                       <  8:113109/120‑1 (MQ=255)                                                ggcggcGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAAt                                                                                                                       >  7:113109/1‑120 (MQ=255)                                                   ggcgTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACAt                                                                                       >  5:167251/1‑149 (MQ=255)                                                       tACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                     <  8:296341/115‑1 (MQ=255)                                                       tACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgc                                                                                                                     >  7:296341/1‑115 (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                                                                                        <  5:235693/120‑1 (MQ=255)                                                                                gcgcgGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                        >  6:235693/1‑120 (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)                                                                                              ggggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCccggc                                                                                                                                    >  8:133289/1‑62 (MQ=255)                                                                                              ggggCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCCGCccggc                                                                                                                                    <  7:133289/62‑1 (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                                                      <  5:173811/121‑1 (MQ=255)                                                                                                                gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCatgatg                                                      >  6:173811/1‑121 (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  <  2:232799/144‑1 (MQ=255)                                                                                                                                        attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACgg  >  1:232799/6‑149 (MQ=255)                                                                                                                                             ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg           >  2:322852/1‑135 (MQ=255)                                                                                                                                             ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACATCGCCAACGGTGTggg           <  1:322852/135‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                          GTCACGCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCG‑CTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGG‑CGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGG  >  NC_002947/1245184‑1245466

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