Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,656,820	C→G	5.9%	Q88H (CAG→CAC)	fdx ←	[2Fe‑2S] ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,656,820	0	C	G	5.9%	278.1 / 6.1	81	Q88H (CAG→CAC)	fdx	[2Fe‑2S] ferredoxin
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base G (2/3);  total (36/45)
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

AGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGT  >  NC_000913/2656672‑2656965                                                                                                                                                     |                                                                                                                                                  aGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATAGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTAGACTACTAAATCTTCGTCGGTAACGCGCGcc                                                                                                                                                   >  7:175936/1‑149 (MQ=255)  gATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCt                                                                                                                                                  <  5:96214/149‑1 (MQ=255)  gATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCt                                                                                                                                                  <  3:105533/149‑1 (MQ=255)   aTCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTg                                                                                                                                                 <  4:198012/149‑1 (MQ=255)    tCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTgg                                                                                                                                                <  7:261580/149‑1 (MQ=255)      gggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCa                                                                                                                                              <  1:236314/149‑1 (MQ=255)      gggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCa                                                                                                                                              <  7:59313/149‑1 (MQ=255)       ggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAg                                                                                                                                             <  7:376901/149‑1 (MQ=255)          tACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCtt                                                                                                                                          <  3:60224/149‑1 (MQ=255)             gCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTaaa                                                                                                                                       <  8:233790/149‑1 (MQ=255)                tCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACgg                                                                                                                                    >  2:34201/1‑149 (MQ=255)                   tACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCtt                                                                                                                                 <  6:150701/149‑1 (MQ=255)                    aCAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCttt                                                                                                                                >  7:233559/1‑149 (MQ=255)                     cAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               <  8:205224/149‑1 (MQ=255)                      aGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               <  7:315256/148‑1 (MQ=255)                      aGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               >  8:315256/1‑148 (MQ=255)                              gCCAATTTCGCGGCTATCGGTCGACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAg                                                                                                                      >  8:73976/1‑149 (MQ=255)                              gCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCa                                                                                                                       <  8:39444/148‑1 (MQ=255)                               ccAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGt                                                                                                                     >  3:301752/1‑149 (MQ=255)                               ccAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGt                                                                                                                     >  3:296808/1‑149 (MQ=255)                                cAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTc                                                                                                                    >  3:90859/1‑149 (MQ=255)                                  aTTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         <  6:202713/142‑1 (MQ=255)                                  aTTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         >  5:202713/1‑142 (MQ=255)                                   tttCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         <  4:19062/141‑1 (MQ=255)                                   tttCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         >  3:19062/1‑141 (MQ=255)                                      cgcgGCTATAGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGATTTGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAg                                                                                                               <  1:64981/148‑1 (MQ=255)                                       gcgGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAgg                                                                                                              <  4:90859/148‑1 (MQ=255)                                          gCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGGTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCttt                                                                                                          <  7:353917/149‑1 (MQ=255)                                           cTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTg                                                                                                         <  7:126949/149‑1 (MQ=255)                                              tCGGTCCACTTAAGTCCCATACTAACCACTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTc                                                                                                       >  7:78347/1‑148 (MQ=255)                                                 gTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAg                                                                                                    >  8:231285/1‑148 (MQ=255)                                                 gTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAg                                                                                                    >  1:202505/1‑148 (MQ=255)                                                 gTCCACTTAAGACCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGc                                                                                                   <  5:317068/149‑1 (MQ=255)                                                    cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                <  5:201612/149‑1 (MQ=255)                                                    cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                <  8:78347/149‑1 (MQ=255)                                                    cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                >  1:258114/1‑149 (MQ=255)                                                    cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                >  3:131955/1‑149 (MQ=255)                                                      cTTAAGTCCAATACTAACCTCTGTTAATGCTCAATCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtc                                                                                              <  7:341257/149‑1 (MQ=255)                                                       ttAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtcg                                                                                             >  4:168309/1‑149 (MQ=255)                                                          aGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtcgtc                                                                                           <  1:213403/148‑1 (MQ=255)                                                          aGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCt                                                                                          <  8:233559/149‑1 (MQ=255)                                                           gTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCtt                                                                                         >  1:319776/1‑149 (MQ=255)                                                           gTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCtt                                                                                         >  1:185973/1‑149 (MQ=255)                                                              ccATAGTAACGTCTGTTAATGCTCACGCGTATGGTTGATAGTGTAACGCGGGATTTCGACTACTATATCTTCGTAGGTAACGCGCGCCTGGCAGCTTTAACGTCTTTCCGGCTCCAGTCCCGTGGCTTTGTCCAGCATGTGGTCTTcc                                                                                       <  1:267572/148‑1 (MQ=255)                                                               cATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTg                                                                                     >  1:231112/1‑149 (MQ=255)                                                                 tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGct                                                                                   <  2:319776/149‑1 (MQ=255)                                                                 tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGc                                                                                    >  6:114001/1‑148 (MQ=255)                                                                 tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGGCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGct                                                                                   >  5:239609/1‑149 (MQ=255)                                                                      aCCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTg                                                                                <  7:329977/147‑1 (MQ=255)                                                                      aCCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAg                                                                              <  7:73976/149‑1 (MQ=255)                                                                         tctGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCtt                                                                           >  8:17077/1‑149 (MQ=255)                                                                            gTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGa                                                                               >  3:753/1‑142 (MQ=255)                                                                            gTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGa                                                                               <  4:753/142‑1 (MQ=255)                                                                            gTTAATGCTCACGCGCATGGTTGATAGCGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGGAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTcc                                                                        >  8:183245/1‑149 (MQ=255)                                                                               aaTGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCgg                                                                      >  8:50952/1‑148 (MQ=255)                                                                                aTGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCCGAGCTTTCCGGc                                                                     <  2:258114/148‑1 (MQ=255)                                                                                 tGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAg                                                                   <  4:102969/149‑1 (MQ=255)                                                                                   cTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGt                                                                  <  1:210620/148‑1 (MQ=255)                                                                                           cATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCaaaa                                                         <  7:231285/149‑1 (MQ=255)                                                                                            aTGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAAc                                                        <  5:272953/149‑1 (MQ=255)                                                                                               gTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCt                                                      >  4:114409/1‑148 (MQ=255)                                                                                                    tAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacg                                                 >  3:137922/1‑148 (MQ=255)                                                                                                         tAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGAt                                           <  1:34201/149‑1 (MQ=255)                                                                                                            cgcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATg                                          >  8:262929/1‑147 (MQ=255)                                                                                                            cgcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCa                                        <  2:202505/149‑1 (MQ=255)                                                                                                             gcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  4:296808/149‑1 (MQ=255)                                                                                                             gcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  4:301752/149‑1 (MQ=255)                                                                                                             gcgGGATTTCGACTACTAAATCATGGTCGGTAACGCGCGCGTGGCAGCTTAAACGGATTTCAGGCTACAGTCCCCAGGCTTTGTAAAGAATGGCGTCTGCCTGCTCTGAGCTTTCCGGGAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  1:58713/149‑1 (MQ=255)                                                                                                                        actactAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtggt                            <  7:309933/149‑1 (MQ=255)                                                                                                                        actactAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGCCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtggt                            >  7:367850/1‑149 (MQ=255)                                                                                                                         ctactaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGgtggtg                           >  1:198498/1‑149 (MQ=255)                                                                                                                            ctaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGc                          <  7:262929/147‑1 (MQ=255)                                                                                                                            ctaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGc                          >  7:129176/1‑147 (MQ=255)                                                                                                                             taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTc                                                                         >  6:84567/1‑99 (MQ=255)                                                                                                                             taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTc                                                                         <  5:84567/99‑1 (MQ=255)                                                                                                                             taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTACGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAg                       >  2:183648/1‑149 (MQ=255)                                                                                                                                 tCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGca                     >  2:290196/1‑147 (MQ=255)                                                                                                                                         ggTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt           >  6:68387/1‑149 (MQ=255)                                                                                                                                           tAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGTCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCTGCAGTGAGTCAAAACCTTCACGAAAGATGAAGTGGAAGGTGGTGCAAGCACAGGATTTTTca         >  3:33947/1‑149 (MQ=255)                                                                                                                                            aaCGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTca         >  2:258510/1‑148 (MQ=255)                                                                                                                                                cgcgcCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGc    <  2:231112/149‑1 (MQ=255)                                                                                                                                                  cgcCTGTCAGCTTAAACGGCTTTCCGGCTCGAGTCACCAGGATTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCgt  <  8:367850/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGT  >  NC_000913/2656672‑2656965

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