Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,491,268	T→A	7.7%	D228V (GAT→GTT)	oxc ←	oxalyl‑CoA decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,491,268	0	T	A	7.7%	273.8 / 9.0	90	D228V (GAT→GTT)	oxc	oxalyl‑CoA decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (46/37);  new base A (5/2);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCG  >  NC_000913/2491120‑2491417                                                                                                                                                     |                                                                                                                                                       aacaaGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTtcatc                                                                                                                                                       >  1:1056065/1‑150 (MQ=255)         tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                                                                              <  2:409581/150‑1 (MQ=255)         tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                                                                              <  1:55634/150‑1 (MQ=255)            cAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAAt                                                                                                                                           >  2:978990/1‑150 (MQ=255)              aCGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAc                                                                                                                                         <  1:1270753/150‑1 (MQ=255)                 cgttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                                                                      >  1:1001951/7‑150 (MQ=255)                 cgttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                                                                      >  1:31026/7‑150 (MQ=255)                  gttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgccg                                                                                                                                     >  2:1532077/6‑150 (MQ=255)                  cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgccg                                                                                                                                     >  2:1314087/1‑150 (MQ=255)                  cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                                                                                     <  2:353682/150‑7 (MQ=38)                  cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                                                                                     <  2:419229/150‑7 (MQ=38)                    gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCcgcg                                                                                                                                   <  1:82055/150‑1 (MQ=255)                    gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCcgcg                                                                                                                                   >  2:1237606/1‑150 (MQ=255)                         tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCttt                                                                                                                              >  1:593809/1‑150 (MQ=255)                          gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                                                                             <  2:1263709/150‑1 (MQ=255)                          gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                                                                             <  2:373270/150‑1 (MQ=255)                           ccAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGc                                                                                                                            >  1:1411995/1‑150 (MQ=255)                            cAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGcc                                                                                                                           <  2:300601/150‑1 (MQ=255)                              gggCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCTATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCaa                                                                                                                         >  2:181861/1‑150 (MQ=255)                                gCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAgg                                                                                                                       <  2:1179256/150‑1 (MQ=255)                                  aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAg                                                                                                                        >  2:634561/1‑147 (MQ=255)                                  aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAg                                                                                                                        <  1:634561/147‑1 (MQ=255)                                   cacgaacgCTCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGCCATTGGAAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGAGCCTTTGCCAAGGAta                                                                                                                    >  2:649690/2‑150 (MQ=255)                                   aacgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAta                                                                                                                    >  2:1393420/1‑150 (MQ=255)                                     cgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                  <  2:303958/150‑1 (MQ=255)                                     cgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                  >  2:1128609/1‑150 (MQ=255)                                          gcgcAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATg                                                                                                             <  2:1311147/150‑1 (MQ=255)                                           cgcAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATgg                                                                                                            >  2:1450900/1‑150 (MQ=255)                                             cAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                                                                          >  2:200338/1‑150 (MQ=255)                                                cTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGtt                                                                                                       >  1:1608699/1‑150 (MQ=255)                                                   ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                                                                                    >  2:512213/1‑150 (MQ=255)                                                   ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                                                                                    <  2:566492/150‑1 (MQ=255)                                                      cAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtt                                                                                                 >  1:131197/1‑150 (MQ=255)                                                           aGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCt                                                                                            <  2:1056065/150‑1 (MQ=255)                                                            gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                                   <  2:692941/142‑1 (MQ=255)                                                            gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                                   >  1:692941/1‑142 (MQ=255)                                                            gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                                                                           <  1:489381/150‑1 (MQ=255)                                                             tGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACATTCAATAAATTCACGAAGCTGTTCATCAGCGT‑GGGAATACGCCGCGCCTTTGCCAAGGCTAATTAAGGGCCGTGCAGCTTTTGCTaa                                                                                          >  2:329366/1‑150 (MQ=255)                                                               gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                        >  1:421541/1‑150 (MQ=255)                                                                aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAg                                                                                       <  2:1280218/150‑1 (MQ=255)                                                                  gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                                                                                     >  2:52193/1‑150 (MQ=255)                                                                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                    >  1:332261/1‑150 (MQ=255)                                                                   cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                    <  1:1237606/150‑1 (MQ=255)                                                                       tCTTCAAGGATCCCTTTCGCCATAGACATTGGAAGGAATGGAATCTGAGCACTTTCAATAAATTCAAGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAtt                                                                                >  1:1467513/1‑150 (MQ=255)                                                                        cTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCATTTGCCAAGGATAATTAATGGCCGTTAAGCTTTTGCTAAAAGCGAAATTg                                                                               >  1:24251/1‑150 (MQ=255)                                                                         ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                              >  2:1116276/1‑150 (MQ=255)                                                                         ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                              >  1:1331417/1‑150 (MQ=255)                                                                          tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcg                                                                             <  1:1532077/150‑1 (MQ=255)                                                                                 tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtga                                                                      >  2:530236/1‑150 (MQ=255)                                                                                 tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtga                                                                      >  2:1379708/1‑150 (MQ=255)                                                                                    cTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                                   >  2:1532073/1‑150 (MQ=255)                                                                                      ttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactgac                                                                 >  1:107395/1‑150 (MQ=255)                                                                                            aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg                                                           >  2:1115915/1‑150 (MQ=255)                                                                                            aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg                                                           >  2:1609113/1‑150 (MQ=255)                                                                                             tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc                                                          >  1:498072/1‑150 (MQ=255)                                                                                             tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc                                                          <  2:1428582/150‑1 (MQ=255)                                                                                              aGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCa                                                         <  2:107395/150‑1 (MQ=255)                                                                                               gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                      <  2:853441/120‑1 (MQ=255)                                                                                               gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                      >  1:853441/1‑120 (MQ=255)                                                                                               gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAt                                                        >  1:1415332/1‑150 (MQ=255)                                                                                                aCATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATg                                                       <  2:972853/150‑1 (MQ=255)                                                                                                  aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc                                                     <  2:360259/150‑1 (MQ=255)                                                                                                  aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc                                                     <  1:1128609/150‑1 (MQ=255)                                                                                                   ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCa                                                    <  1:542718/150‑1 (MQ=255)                                                                                                   ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCa                                                    <  2:421541/150‑1 (MQ=255)                                                                                                       cAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaataa                                                >  2:1052747/1‑150 (MQ=255)                                                                                                                aatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGa                                       >  2:1412230/1‑150 (MQ=255)                                                                                                                     gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAt                                  >  2:996130/1‑150 (MQ=255)                                                                                                                       gCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAttt                                <  2:827602/150‑1 (MQ=255)                                                                                                                         aCTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTc                              <  1:1052747/150‑1 (MQ=255)                                                                                                                            ttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAAc                           <  1:1532073/150‑1 (MQ=255)                                                                                                                             tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACt                          <  1:1393420/150‑1 (MQ=255)                                                                                                                             tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACt                          >  1:377456/1‑150 (MQ=255)                                                                                                                               ggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                        >  2:1172389/6‑150 (MQ=255)                                                                                                                               ggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                        >  1:749553/6‑150 (MQ=255)                                                                                                                               aaTAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                        >  2:342037/1‑150 (MQ=255)                                                                                                                                gcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTa                       >  1:317047/5‑150 (MQ=255)                                                                                                                                 tAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                      >  1:988604/1‑150 (MQ=255)                                                                                                                                  aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAAc                     <  1:181861/150‑1 (MQ=255)                                                                                                                                   aaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACa                    <  2:1608699/150‑1 (MQ=255)                                                                                                                                       cACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCg                >  2:959058/1‑150 (MQ=255)                                                                                                                                        aCGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGt               <  2:1497213/150‑1 (MQ=255)                                                                                                                                         cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg              >  1:554755/1‑150 (MQ=255)                                                                                                                                         cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg              >  2:1448306/1‑150 (MQ=255)                                                                                                                                         cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg              >  2:278372/1‑150 (MQ=255)                                                                                                                                         acaaGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGt               >  1:986489/3‑150 (MQ=255)                                                                                                                                         acaaGCTGATGAACAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg              >  2:453704/3‑150 (MQ=255)                                                                                                                                          gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg             <  1:959058/150‑1 (MQ=255)                                                                                                                                          gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg             <  2:1432602/150‑1 (MQ=255)                                                                                                                                            aGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGtt           <  2:332261/150‑1 (MQ=255)                                                                                                                                               tGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAc        <  1:329366/150‑1 (MQ=255)                                                                                                                                                   catcaGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCtt    >  1:204836/1‑150 (MQ=255)                                                                                                                                                     tcaGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCg  >  2:137402/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AACAAGCATGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCG  >  NC_000913/2491120‑2491417

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