Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	425,348	A→C	7.3%	intergenic (‑39/+18)	PP_0349 ← / ← PP_0350	membrane protein/ferrichrome‑iron receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	425,348	0	A	C	7.3%	172.0 / 2.0	62	intergenic (‑39/+18)	PP_0349/PP_0350	membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (3/2);  total (32/30)
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

ATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCC  >  NC_002947/425206‑425495                                                                                                                                               |                                                                                                                                                     aTAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcgg                                                                                                                                                <  8:167972/149‑1 (MQ=255)  tAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCTCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                               <  8:133239/149‑1 (MQ=255)  tAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                               >  3:152400/1‑149 (MQ=255)   aaGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                              >  7:154409/1‑149 (MQ=255)    aGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGccc                                                                                                                                             <  4:240672/149‑1 (MQ=255)     ggggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcccc                                                                                                                                            >  7:147050/1‑149 (MQ=255)      gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCACGGCGGCCCCa                                                                                                                                           >  3:179493/1‑149 (MQ=255)       ggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAg                                                                                                                                          >  2:279997/1‑149 (MQ=255)              cGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctccct                                                                                                                                    >  6:166258/1‑148 (MQ=255)               ggTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcact                                                                                                                                    >  4:52140/1‑147 (MQ=255)                   aCTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGaa                                                                                                                               <  8:154409/148‑1 (MQ=255)                       ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                              <  4:315153/129‑1 (MQ=255)                       ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                              >  3:315153/1‑129 (MQ=255)                       ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGATGTAGGTTGGAACCCGGGGCCCCAACGCGGCCCCAGCTCACTCACAACTTg                                                                                                                           >  8:17563/1‑148 (MQ=255)                          gCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAg                                                                                                                       >  6:120758/1‑149 (MQ=255)                             aCAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTgg                                                                                                                    <  8:44531/149‑1 (MQ=255)                               agcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccggg                                                                                                                                         >  3:78010/1‑123 (MQ=255)                               agcagGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccggg                                                                                                                                         <  4:78010/126‑4 (MQ=255)                                   ggAACAGGGTACAGACGAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTgg                                                                                                              <  5:189589/149‑1 (MQ=255)                                    gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTTGTTCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCATAACTTGTAATTGGCGCTgac                                                                                                             >  6:281346/1‑147 (MQ=255)                                    gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGGCGCCACGGGGCCCCCGCGCCCGCAGGACGTGTATGTGGGGCgggg                                                                                                             >  4:8190/1‑148 (MQ=255)                                              cAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCt                                                                                                   <  4:152400/148‑1 (MQ=255)                                              cAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGGTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCt                                                                                                   <  5:281346/148‑1 (MQ=255)                                               aGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcact                                                                                                                                    <  1:329208/115‑1 (MQ=255)                                                  ccAGGCTGCTCCAGGTGTGGATCTAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTgcgc                                                                                               <  7:126695/148‑1 (MQ=255)                                                     ggCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCAt                                                                                           >  4:84252/1‑149 (MQ=255)                                                                gtgtGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGtagtag                                                                                <  7:17563/149‑1 (MQ=255)                                                                      ataaaGGACCAGAGGCGGATGGTGGGTATCTTCATCACGGTACCGATTGATGTCGGATGGAACCCGGGGCCCCAAGGCCGCCCCAGATCACACACAAAATGTAAGTTTCGCATGCGACA‑ACGATGCACAGATCGCCGTAGTATCagtag                                                                          >  4:125536/4‑149 (MQ=255)                                                                      aTCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                          >  4:51442/1‑149 (MQ=255)                                                                      aTCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                          >  7:150077/1‑149 (MQ=255)                                                                       tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGACCCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACGTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                          >  5:169332/1‑148 (MQ=255)                                                                       tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTGATCGCCGTAGTAGCagtag                                                                          >  7:14675/1‑148 (MQ=255)                                                                       tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCa                                                                                            <  5:291941/130‑1 (MQ=255)                                                                       tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCa                                                                                            >  6:291941/1‑130 (MQ=255)                                                                       tCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCagtagt                                                                         <  3:52140/149‑1 (MQ=255)                                                                         aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAgtagta                                                                        <  5:120758/148‑1 (MQ=255)                                                                          aGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAgg                                                                      >  3:17664/1‑149 (MQ=255)                                                                             aCCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGt                                                                   <  3:84252/149‑1 (MQ=255)                                                                               cAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTc                                                                  >  7:166155/1‑148 (MQ=255)                                                                                           gtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                      <  8:147050/148‑1 (MQ=255)                                                                                               ggCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAAt                                                 >  5:342861/1‑149 (MQ=255)                                                                                               ggCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAAt                                                 >  1:5821/1‑149 (MQ=255)                                                                                               ggCTCTTCATCACGGTACCGATTGTTGGCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGGCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAAt                                                 <  8:102784/149‑1 (MQ=255)                                                                                                      catcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCtt                                             <  3:8190/146‑1 (MQ=255)                                                                                                        tcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCCCTCACTCAGAACATGTAAGTGGCGCTGGCGACA‑AAGCTGCGCGCATAGCCGTAGTAGCAGTAGTAGGAGTAACAAGTAGAGAGGTAATACTTGtt                                          >  4:185463/1‑147 (MQ=255)                                                                                                                      gttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGc                           <  6:342861/148‑1 (MQ=255)                                                                                                                       ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGcc                                                                                       >  7:315798/1‑87 (MQ=255)                                                                                                                       ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGcc                                                                                       <  8:315798/87‑1 (MQ=255)                                                                                                                       ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCgg                         <  8:150077/149‑1 (MQ=255)                                                                                                                       ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCg                          <  3:185463/148‑1 (MQ=255)                                                                                                                                 tgagctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggtg                             <  2:254391/129‑1 (MQ=255)                                                                                                                                 tgagctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggtg                             >  1:254391/7‑135 (MQ=255)                                                                                                                                 gAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGGGGCGGTGACCGcga                <  4:179493/148‑1 (MQ=255)                                                                                                                                     ccGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgaccgac           >  1:86933/1‑149 (MQ=255)                                                                                                                                       ggggCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGAcgc         <  4:17664/149‑1 (MQ=255)                                                                                                                                            cGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCt    >  1:223171/1‑149 (MQ=255)                                                                                                                                             gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACAGACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCt    >  1:87064/1‑149 (MQ=255)                                                                                                                                             gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGAcc                    <  2:136615/132‑1 (MQ=255)                                                                                                                                             gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGAcc                    >  1:136615/1‑132 (MQ=255)                                                                                                                                             gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc   >  6:202096/1‑149 (MQ=255)                                                                                                                                             gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc   <  5:166258/149‑1 (MQ=255)                                                                                                                                              cAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTcc  <  2:5821/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     ATAAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGCTGGCGACA‑ACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTCC  >  NC_002947/425206‑425495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑