Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,061,818	G→C	8.3%	W323C (TGG→TGC) ‡	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,818	0	G	C	8.3%	252.4 / 8.1	88	W323C (TGG→TGC) ‡	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base C (6/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCC  >  NC_000913/3061674‑3061966                                                                                                                                                 |                                                                                                                                                     cttcttcTGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGcc                                                                                                                                                  <  2:222680/149‑1 (MQ=255) attcttcGGCATCGGCATGGATCTGTTTATGACCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGcc                                                                                                                                                  <  2:660056/148‑1 (MQ=255)  ttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCTGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTTGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCct                                                                                                                                                 <  1:1474464/149‑1 (MQ=255)   tcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctg                                                                                                                                                >  2:1256471/1‑149 (MQ=255)    cttcGTCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                                                                                                                               <  2:398641/149‑1 (MQ=255)    cttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                                                                                                                               <  1:1124287/149‑1 (MQ=255)     ttcGGCATCGGCATGGAGCTGTTTATGAACGTCGCCATGTTGCGGGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                              <  2:932647/149‑2 (MQ=255)      tcggcatcggcaTGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCTCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                                                                                                                                             <  2:533145/149‑1 (MQ=255)      tcggcatcggcaTGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                                                                                                                                             <  1:622010/149‑1 (MQ=255)           atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgaggtctgg                                                                                                                                        <  1:734735/149‑8 (MQ=255)              ggcatGGATCTGATTATGACCGTGGCCAGGGTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                                                                                                                                     <  2:992793/149‑1 (MQ=255)               gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGCCCCGAAATCACTGGCGCTGCTTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGCACAGGAt                                                                                                                                    >  1:306766/1‑149 (MQ=255)               gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTCCCCACTGCCAGACCTCAGGCGGGAGCCTGACTGAACCGGAt                                                                                                                                    >  2:1031401/1‑149 (MQ=255)               gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                    <  2:1331100/149‑1 (MQ=255)               gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                    >  2:254064/1‑149 (MQ=255)               gcatGGATCTCTTTATGAACGTCGCCATGTTGCGTGCGGCACGGTATTTATTGAGCGAAGCTGTCAGTGGAAATTGCGCACAGGACCCGAAAACACTGGCGCTGCGTACCCACTCCCAGACCTCAGGCTGCAGCCTGACTGAACAGGAt                                                                                                                                    >  1:280195/1‑149 (MQ=255)                 atggtgcTGATTATGATCGTCGCCATTGTGCGTCCGGCACTTTATATATTGAGCGAAGCGGTGAGTGGATTTGGCGCAGAGGACCCGATATCACTGGCTCTTCGTACCCTCTGCCAGACCTCAGGGGGGAGCCTGACTGAACAGGATcc                                                                                                                                  <  2:245446/143‑1 (MQ=255)                     aTCTGGTTATGACCGTCGCCATGTTGCGCGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                              <  2:159553/149‑1 (MQ=255)                        ggTTTATGAACGTCGCGATGTTGCGGGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaac                                                                                                                           <  1:280201/148‑1 (MQ=255)                         gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                          >  1:1484706/1‑149 (MQ=255)                         gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                          >  1:242131/1‑149 (MQ=255)                         gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                          >  2:1517538/1‑149 (MQ=255)                           ttATGTACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGCAGCGGGGACTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaac                                                                                                                        <  2:855087/149‑1 (MQ=255)                            tATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACg                                                                                                                       <  1:1256471/149‑1 (MQ=255)                               gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa                                                                                                                    >  1:1167214/1‑149 (MQ=255)                                aaCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTAt                                                                                                                   >  1:40442/1‑149 (MQ=255)                                 aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc                                                                                                                  <  1:34847/149‑1 (MQ=255)                                      gCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcac                                                                                                             <  2:1303366/149‑1 (MQ=255)                                        cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGCGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacca                                                                                                           >  1:302616/1‑149 (MQ=255)                                        cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcaccc                                                                                                           >  1:79464/1‑148 (MQ=255)                                           gTTGCGTGCGGCACGGTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCCAAACCACTGCCAGACCTCAGACTGCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacaacaa                                                                                                        >  2:163833/1‑147 (MQ=255)                                                    ggCACGTTATTTATCGAGCGAAGCTGGCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGACgcgc                                                                                               <  2:1376689/149‑1 (MQ=255)                                                     gctcGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCGAGCCGGACTGAACAGGATCCGTATAACCACGTTATCCGCACCCCCATGGCAgcgcg                                                                                              >  2:1033585/4‑148 (MQ=255)                                                     gCACGTTATTTATGGAGCGAAGAGGTCAGTGGATTTGGAGCACAGGACCCGAAATCACTGGCGCGGCGTACCCACTGCCAGACCTCAGGCTGTAGCCTGACTGATCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgct                                                                                              <  2:487867/149‑1 (MQ=255)                                                          ttatttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggctg                                                                                         >  1:1118632/1‑149 (MQ=255)                                                          ttatttaaGGAGCGAAGCGGTGAGTGGATTTGGCGCACAGGACCCGAAATCACTTGCGATGCGTACCCTCTGCCATACATCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATACGCACCACCATTGAAGCgctggctg                                                                                         <  2:223590/149‑1 (MQ=255)                                                            atttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCg                                                                                       <  1:254064/149‑1 (MQ=255)                                                                  ggAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGACCCTGACTGAAAAGGATCCGTATAACAACGTTAACCACACCCGCATTGAAGCGCTGGCAGCGACTCTg                                                                                 >  1:1037667/1‑149 (MQ=255)                                                                     gCGAAGCGGTCAGTGGTTTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAACTCAGGCTGGAGCCTGACTTAACAGGATCCGTATAACAACGTTATCCGCACCACCGTTGAAGCGCTGGCTGCGACGCTGGGc                                                                              <  2:1399559/149‑1 (MQ=255)                                                                       gAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGTCGCGGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCgg                                                                            <  1:139339/149‑1 (MQ=255)                                                                        aaGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCTTGACTGAACAGGATCCGTATAACAAGGTTATCCGGAACACCATTGAAGCGCTGGCTGCGACGCTGGGCGGt                                                                           >  2:863952/1‑149 (MQ=255)                                                                         aGCGGTCAGTGGATTTGGCGCACTGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGTGCCTGACTGAACAGGATCCGTATTACAACGTTATCCGCACCACCAGTGAAGCGCTGGCTGCGCCGCTGGGCGGTa                                                                          >  1:266409/1‑149 (MQ=255)                                                                         aGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTa                                                                          >  1:729818/1‑149 (MQ=255)                                                                             gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCTCTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCa                                                                      <  2:1341812/149‑1 (MQ=255)                                                                               cAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGt                                                                    <  1:461734/149‑1 (MQ=255)                                                                                  tGGATTTGGCGCACAGGACCCGAAATCACTGGCTCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGGTGCGACGCTGGGCGGTACTCAGTCAc                                                                 >  1:1352043/1‑149 (MQ=255)                                                                                   ggATTTGGTGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                <  2:779471/149‑1 (MQ=255)                                                                                   ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGACCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                <  1:646521/149‑1 (MQ=255)                                                                                   ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                <  2:779215/149‑1 (MQ=255)                                                                                    gATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTg                                                               >  2:167716/1‑149 (MQ=255)                                                                                    gATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTg                                                               >  1:435594/1‑149 (MQ=255)                                                                                      tttGGCGCACAGGACCCGAAATCACTGGCGCTGCGTCCCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAAAGTTATCCGCCCCACCATTGAAGCGCTGGCGGCGACGCTGGGCGGTACTCAGTCGCGGCa                                                             >  2:1041050/1‑149 (MQ=255)                                                                                      tttGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCAGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGCAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCa                                                             >  1:414533/1‑149 (MQ=255)                                                                                       ttGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCAt                                                            <  1:1183327/149‑1 (MQ=255)                                                                                            taacaGGACCCGAAATCACTGGCGCTGCGTACCAACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCaa                                                       <  1:1026175/147‑1 (MQ=255)                                                                                             cacaGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAc                                                      <  1:689905/149‑1 (MQ=255)                                                                                             cacaGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCCGGGCGGTACTCAGTCACTGCATACCAAc                                                      >  1:142817/1‑149 (MQ=255)                                                                                              acaGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGCCTGAACAGGATCCGTATAACAACGTAATCCGCACCACCATTGAAGCGCGGGCTGCGACGCTGGGCGGTCAGCAGTCACTGCCTACCAACg                                                     >  2:1019321/1‑149 (MQ=255)                                                                                              acaGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACg                                                     <  2:40442/149‑1 (MQ=255)                                                                                               caGGACCCGAAATCACTGGCGCTGCGTACCCACAGCCAGACCTCAGGCGGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCGGGCGGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGc                                                    <  1:925148/149‑1 (MQ=255)                                                                                                     ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                              <  2:752954/149‑1 (MQ=255)                                                                                                     ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                              <  1:1197254/149‑1 (MQ=255)                                                                                                      cGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGAc                                             >  1:1210739/1‑149 (MQ=255)                                                                                                       gAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACg                                            >  1:817378/1‑149 (MQ=255)                                                                                                        aaaTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGa                                           >  2:254625/1‑149 (MQ=255)                                                                                                             aCTGGCGCTGCGTACCCACTGCCAGACATCAGGCTGGAGACTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGAGGCAgcga                                      >  2:923032/1‑148 (MQ=255)                                                                                                                 gcgcTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGCAGCGCTTgg                                  >  1:724324/1‑149 (MQ=255)                                                                                                                     tGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTg                              <  2:729818/149‑1 (MQ=255)                                                                                                                      gCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGc                             <  2:1167214/149‑1 (MQ=255)                                                                                                                          aCCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACCTTATCCGCCCCCCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCACCGCCTTTGACGAAGCGCAGGGTTTGCCTag                         >  1:298805/1‑148 (MQ=255)                                                                                                                           cccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTAcc                        >  1:1101860/1‑149 (MQ=255)                                                                                                                            ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                       >  2:1499843/1‑149 (MQ=255)                                                                                                                              aCTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAt                     >  1:426049/1‑149 (MQ=255)                                                                                                                                tGCCAGACCTGAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGTAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                   <  1:536533/149‑1 (MQ=255)                                                                                                                                tGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTTAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                   <  1:1384485/149‑1 (MQ=255)                                                                                                                                  ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGCCGCTGGGCGGTACTCAGTCACTGCATACCAACGCGTTTGACGAAGCGCTTGGTTTGCGTACCGATTtct                 >  2:562808/1‑149 (MQ=255)                                                                                                                                  ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtct                 <  2:1210739/149‑1 (MQ=255)                                                                                                                                  acAGACCTCAGTCTGGAGCCTGACTGAACAGGATCAGTATAACAACTTTATCCGCACCACCATTGAAGCGCTGGCTTCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtct                 <  1:382539/148‑1 (MQ=255)                                                                                                                                   cAGACCTCAGGCGGGAGCCGGACTGTACAGGATCCGGATTACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCTTACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtctc                <  2:1037667/149‑1 (MQ=255)                                                                                                                                    ttcagtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCTCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa               >  2:482820/6‑149 (MQ=255)                                                                                                                                      aCCTCAGGCTGGAGCCTGACTGCACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCGTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGc             >  1:848505/1‑149 (MQ=255)                                                                                                                                      aCCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTGTGACGCAGCGCTTGGTTTGCCTACCGATTTCTCAGc             >  1:44267/1‑149 (MQ=255)                                                                                                                                       agtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGCAGCCCTGGCTGCGACGCAGGGCGCTACTCAGTCACTGCATACCAACGCGTATGACGACGCGCATGGTGTCCCGACCGATTACTCAGcc            >  1:875686/3‑148 (MQ=255)                                                                                                                                        gtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATTCCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc           >  2:630862/2‑149 (MQ=255)                                                                                                                                        gtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCGTTGACGAAGCGCTTGGTTTGCCTACCGAGTTCTCAGCAc           >  2:1154649/2‑149 (MQ=255)                                                                                                                                        cTCCGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc           <  1:1031401/149‑1 (MQ=255)                                                                                                                                        cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAc           >  2:462707/1‑149 (MQ=255)                                                                                                                                           aGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTAGACGAAGCGCTTGGTTTGCCTCCCGATTTCTCAGCACGCa        >  2:915979/1‑149 (MQ=255)                                                                                                                                                ggAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGcc   <  1:1499843/149‑1 (MQ=255)                                                                                                                                                 gAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGccc  <  2:104980/149‑1 (MQ=255)                                                                                                                                                 gAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCGGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGccc  <  2:1395938/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCC  >  NC_000913/3061674‑3061966

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