Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,061,817	G→C	7.1%	W323S (TGG→TCG) ‡	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,817	0	G	C	7.1%	333.8 / 14.3	98	W323S (TGG→TCG) ‡	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (47/44);  new base C (4/3);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCC  >  NC_000913/3061670‑3061965                                                                                                                                                      |                                                                                                                                                     cGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTgg                                                                                                                                                     <  3:5254/149‑1 (MQ=255)  gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGa                                                                                                                                                    <  1:752091/149‑1 (MQ=255)  gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGa                                                                                                                                                    >  2:85665/1‑149 (MQ=255)  gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCCGGa                                                                                                                                                    <  4:633618/149‑1 (MQ=255)   ttcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAg                                                                                                                                                   >  1:466206/1‑149 (MQ=255)    tcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGc                                                                                                                                                  <  4:276755/149‑1 (MQ=255)    tcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGc                                                                                                                                                  <  4:239479/149‑1 (MQ=255)      ttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCct                                                                                                                                                >  2:719521/1‑149 (MQ=255)      ttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCct                                                                                                                                                <  2:715451/149‑1 (MQ=255)        cttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCctga                                                                                                                                              >  4:239492/1‑149 (MQ=255)        cttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCctga                                                                                                                                              <  1:615772/149‑1 (MQ=255)         ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                             <  2:466206/149‑1 (MQ=255)         ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                             <  2:599296/149‑1 (MQ=255)         ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                             <  4:770419/149‑2 (MQ=255)         ttaggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                             <  4:524478/146‑1 (MQ=255)            ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTCCAGCctgaggtc                                                                                                                                          <  4:307718/149‑5 (MQ=255)             gcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGaa                                                                                                                                         >  3:769660/1‑149 (MQ=255)             gcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTAGTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGaa                                                                                                                                         >  4:617054/1‑149 (MQ=255)               atcggcttGGATCTGATTATGAACGTCGCAATGGTGCTTACGACAGGTTATTTATGTTCCGGAGCGGTCAGTAG‑ATTCGGCGCACTGGTTCCGACATCACTGGCGCTGCGTACCCAC‑TGCTAGACCTGAGACTGGAGTCCGGCGGAACa                                                                                                                                       <  2:588081/149‑1 (MQ=255)                    catGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                                                                                                                                  <  4:769660/149‑1 (MQ=255)                         aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                             >  2:138814/1‑149 (MQ=255)                          tCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGA‑CTAAGCGGTCAGGGGAATTTGCGGCAAGGGACCCGAAAT‑AATGGCCGTGGGGACCCACTTGCAGGACCTC‑GGGTGGAGGCCTACTGAACAGGAATCGTATaa                                                                                                                           >  2:150209/1‑149 (MQ=255)                              tttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaa                                                                                                                        <  1:818388/149‑1 (MQ=255)                              tttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaa                                                                                                                        <  3:312951/149‑1 (MQ=255)                               ttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaac                                                                                                                       <  1:96394/149‑1 (MQ=255)                                 aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGt                                                                                                                     >  2:115562/1‑149 (MQ=255)                                   gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                                                                                                                                    <  3:718832/132‑1 (MQ=255)                                   gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                                                                                                                                    >  4:718832/1‑132 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTCTCCGca                                                                                                             >  3:588943/1‑149 (MQ=255)                                         cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGca                                                                                                             <  3:25649/149‑1 (MQ=255)                                          gCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcac                                                                                                            <  1:138814/149‑1 (MQ=255)                                           ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATCACAACGTTATCCGcacc                                                                                                           >  2:626482/1‑149 (MQ=255)                                           ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacc                                                                                                           >  3:461551/1‑149 (MQ=255)                                             aTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcaccac                                                                                                         <  4:261999/149‑1 (MQ=255)                                             aTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcaccac                                                                                                         <  2:484780/149‑1 (MQ=255)                                                ttGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAt                                                                                                      <  1:160890/149‑1 (MQ=255)                                                 tgcgcgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                     <  2:329467/149‑1 (MQ=255)                                                    gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                  <  2:225519/149‑1 (MQ=255)                                                       cgGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAgcg                                                                                               >  3:658468/1‑149 (MQ=255)                                                         gCACGTTATTTATGGAGCGAAGCGGTCAGTAG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGAACTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACAACCATTGAAgcgaa                                                                                             <  2:82566/149‑3 (MQ=255)                                                          cACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctg                                                                                            >  3:665355/1‑149 (MQ=255)                                                           aCGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctgg                                                                                           >  1:894560/1‑149 (MQ=255)                                                           aCGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctgg                                                                                           >  1:358347/1‑149 (MQ=255)                                                            cGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggc                                                                                          <  1:118097/149‑1 (MQ=255)                                                             gTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggct                                                                                         >  1:583707/1‑149 (MQ=255)                                                                atttatGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCg                                                                                      <  4:606719/149‑1 (MQ=255)                                                                 tttatGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGa                                                                                     <  2:677294/149‑1 (MQ=255)                                                                    atGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGc                                                                                  <  2:382468/149‑1 (MQ=255)                                                                    atGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                             >  1:97210/1‑89 (MQ=255)                                                                    atGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                             <  2:97210/90‑2 (MQ=255)                                                                     tGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGATGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCAACACCATTGAAGCGCTGGCTGCGACGCt                                                                                 <  4:614813/149‑1 (MQ=255)                                                                      ggAGCGAAGAGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAGTGAAGCGCTGGCTGCGACGCTg                                                                                >  1:769252/1‑149 (MQ=255)                                                                         gCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGc                                                                             <  1:115562/149‑1 (MQ=255)                                                                          cGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCg                                                                            >  1:505517/1‑149 (MQ=255)                                                                            aaGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGt                                                                          >  4:626135/1‑149 (MQ=255)                                                                                 gTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCa                                                                     >  3:820924/1‑149 (MQ=255)                                                                                 gTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCa                                                                     <  2:719471/149‑1 (MQ=255)                                                                                    aGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTc                                                                  >  2:332640/1‑149 (MQ=255)                                                                                              ggCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATAc                                                         <  3:275911/149‑1 (MQ=255)                                                                                              ggCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATAc                                                         <  1:455496/149‑1 (MQ=255)                                                                                                  cacaGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAc                                                     >  4:611155/1‑149 (MQ=255)                                                                                                   acaGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACg                                                    >  3:794975/1‑149 (MQ=255)                                                                                                     aGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAAGGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATGCCATCGcc                                                  >  1:697763/1‑149 (MQ=255)                                                                                                       gACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCtt                                                >  3:165556/1‑149 (MQ=255)                                                                                                       gACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCtt                                                >  4:823826/1‑149 (MQ=255)                                                                                                        aCCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                               <  4:877482/149‑1 (MQ=255)                                                                                                          ccGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             <  2:583707/149‑1 (MQ=255)                                                                                                          ccGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             <  2:588878/149‑1 (MQ=255)                                                                                                          ccGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             >  3:286317/1‑149 (MQ=255)                                                                                                          ccGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCACTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             >  2:541321/1‑149 (MQ=255)                                                                                                           cGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGAc                                            <  4:658468/149‑1 (MQ=255)                                                                                                            gAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACg                                           >  1:648792/1‑149 (MQ=255)                                                                                                              aaTCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGaa                                         >  2:249762/1‑149 (MQ=255)                                                                                                              aaTCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACACCTTTGACGaa                                         <  3:692813/149‑1 (MQ=255)                                                                                                               aTCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATGCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             >  3:485581/1‑144 (MQ=255)                                                                                                               aTCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATGCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             <  4:485581/144‑1 (MQ=255)                                                                                                                    tGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCtt                                   <  4:588943/149‑1 (MQ=255)                                                                                                                     ggCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTg                                  <  1:719521/149‑1 (MQ=255)                                                                                                                     ggCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTg                                  >  3:33046/1‑149 (MQ=255)                                                                                                                          tGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTg                             >  1:514928/1‑149 (MQ=255)                                                                                                                                 ccAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                      >  1:38008/1‑149 (MQ=255)                                                                                                                                 ccAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                      >  1:226172/1‑149 (MQ=255)                                                                                                                                        ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtct                >  4:140173/1‑149 (MQ=255)                                                                                                                                          aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa              <  4:854008/149‑1 (MQ=255)                                                                                                                                          aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa              >  3:192416/1‑149 (MQ=255)                                                                                                                                          aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa              <  2:24457/149‑1 (MQ=255)                                                                                                                                          aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa              <  3:611155/149‑1 (MQ=255)                                                                                                                                          aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa              >  3:596602/1‑149 (MQ=255)                                                                                                                                              gtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:646925/2‑149 (MQ=255)                                                                                                                                              gtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:736892/2‑149 (MQ=255)                                                                                                                                              gtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:128400/2‑149 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:856650/1‑149 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:316533/1‑149 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          <  1:278943/149‑1 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  2:501221/1‑149 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  1:144609/1‑149 (MQ=255)                                                                                                                                              cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc          >  2:98038/1‑149 (MQ=255)                                                                                                                                               tCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg         >  1:548242/1‑149 (MQ=255)                                                                                                                                               tCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg         >  4:677062/1‑149 (MQ=255)                                                                                                                                                      ggAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGcc  <  2:894560/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGG‑ATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCAC‑TGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCC  >  NC_000913/3061670‑3061965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑