Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,061,817	G→C	7.7%	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.7%	242.5 / 10.7	88	W323S (TGG→TCG) ‡	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (39/42);  new base C (5/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

TCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCC  >  NC_000913/3061669‑3061965                                                                                                                                                     |                                                                                                                                                     tCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTg                                                                                                                                                      >  1:104980/1‑149 (MQ=255)      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                                                                                                                         >  2:1517538/1‑149 (MQ=255)                              gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                         >  1:1484706/1‑149 (MQ=255)                              gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                         >  1:242131/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                                             <  1:1197254/149‑1 (MQ=255)                                                                                                          ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                             <  2:752954/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)                                                                                                                                                     |                                                                                                                                                     TCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCC  >  NC_000913/3061669‑3061965

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