Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	532,886	A→C	N146H (AAC→CAC)	allA →	ureidoglycolate lyase, releasing urea

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	532,886	0	A	C	100.0%	312.1 / NA	76	N146H (AAC→CAC)	allA	ureidoglycolate lyase, releasing urea
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (36/40);  total (36/40)

GACGACAAACCAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGAC  >  NC_000913/532739‑533026                                                                                                                                                    |                                                                                                                                             gacgacAAACCAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCAc                                                                                                                                            <  2:651061/150‑1 (MQ=255) gacgacAAACCAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACt                                                                                                                                           <  2:272546/151‑1 (MQ=255)      cAAACCAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGAt                                                                                                                                      <  1:951586/151‑1 (MQ=255)          ccAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTg                                                                                                                                  >  1:363308/1‑151 (MQ=255)            aGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGaa                                                                                                                                <  1:53615/151‑1 (MQ=255)                cTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGt                                                                                                                             >  1:1062239/1‑150 (MQ=255)                  gTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTAtt                                                                                                                          <  1:694485/151‑1 (MQ=255)                  gTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTAtt                                                                                                                          <  1:327138/151‑1 (MQ=255)                    cAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGTATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCCGCAGTGACCACTGTGATGTTGAAAGTATTc                                                                                                                         <  1:808554/150‑1 (MQ=255)                    cAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTcc                                                                                                                        <  2:590385/151‑1 (MQ=255)                    cAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTc                                                                                                                         <  1:759895/150‑1 (MQ=255)                       cGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGa                                                                                                                     <  2:606230/151‑1 (MQ=255)                        gCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGa                                                                                                                     >  2:179297/1‑150 (MQ=255)                         cTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAAc                                                                                                                   >  2:505902/1‑151 (MQ=255)                           gCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACa                                                                                                                  <  2:986029/150‑1 (MQ=255)                              ggCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGaa                                                                                                              <  2:1099268/151‑1 (MQ=255)                                 gTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTTACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAAct                                                                                                            >  2:314878/1‑150 (MQ=255)                                 gTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAActc                                                                                                           <  2:367362/151‑1 (MQ=255)                                   ttATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTg                                                                                                         >  1:332927/1‑151 (MQ=255)                                   ttATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTg                                                                                                         >  2:748027/1‑151 (MQ=255)                                    tATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGt                                                                                                        >  1:1038776/1‑151 (MQ=255)                                    tATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACGTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAGctct                                                                                                          <  2:216078/149‑1 (MQ=255)                                       cACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGtttt                                                                                                     <  2:519922/151‑1 (MQ=255)                                       cACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGttt                                                                                                      <  1:159213/150‑1 (MQ=255)                                       cACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGttt                                                                                                      <  1:787836/150‑1 (MQ=255)                                       cACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGttt                                                                                                      <  1:830755/150‑1 (MQ=255)                                         ccAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTg                                                                                                    <  2:877960/150‑1 (MQ=255)                                             cGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGa                                                                                               <  1:512403/151‑1 (MQ=255)                                              ggCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGAc                                                                                              <  1:117935/151‑1 (MQ=255)                                                cGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGTCAGTGACCACTGTGATGGTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGa                                                                                               >  2:490670/1‑148 (MQ=255)                                                 gAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGcc                                                                                           >  1:643622/1‑151 (MQ=255)                                                  aaCAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCt                                                                                          <  1:797715/151‑1 (MQ=255)                                                  aaCAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCt                                                                                          >  1:898307/1‑151 (MQ=255)                                                     aGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGc                                                                                        >  2:1055754/1‑150 (MQ=255)                                                          gTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCg                                                                                   >  2:1060697/1‑150 (MQ=255)                                                            gAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGAc                                                                                 <  1:418356/150‑1 (MQ=255)                                                            gAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACt                                                                                >  1:61491/1‑151 (MQ=255)                                                              aTTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTg                                                                              <  1:672706/151‑1 (MQ=255)                                                                tACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGc                                                                             <  2:5818/150‑1 (MQ=255)                                                                  ccATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCAt                                                                           <  2:345960/150‑1 (MQ=255)                                                                  ccATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATa                                                                          <  2:344490/151‑1 (MQ=255)                                                                           cGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACt                                                                 >  1:739005/1‑151 (MQ=255)                                                                               tGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAAtt                                                             >  1:698083/1‑151 (MQ=255)                                                                                ggCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAAtt                                                             <  2:483988/150‑1 (MQ=255)                                                                                ggCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTg                                                            <  2:202134/151‑1 (MQ=255)                                                                                  cATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGtt                                                          <  2:8431/151‑1 (MQ=255)                                                                                    tCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTc                                                         >  1:1030153/1‑150 (MQ=255)                                                                                    tCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTc                                                         >  1:682841/1‑150 (MQ=255)                                                                                      aCCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCAtt                                                      >  2:479894/1‑151 (MQ=255)                                                                                       cccACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTg                                                     <  2:1042539/151‑1 (MQ=255)                                                                                           cTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGttt                                                  >  1:321310/1‑150 (MQ=255)                                                                                             tttCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAATATAAACTAATTGTTCATTGTTTATg                                               >  2:482495/1‑151 (MQ=255)                                                                                              ttCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGc                                              >  1:724968/1‑151 (MQ=255)                                                                                                 gCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTc                                            >  1:670666/1‑150 (MQ=255)                                                                                                     ggCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACt                                         <  1:1035412/149‑1 (MQ=255)                                                                                                     ggCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTg                                       >  1:747844/1‑151 (MQ=255)                                                                                                         gcgcGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAgg                                   <  2:455473/151‑1 (MQ=255)                                                                                                           gcgTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTc                                 >  1:480412/1‑151 (MQ=255)                                                                                                             gTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCgg                               >  2:583731/1‑151 (MQ=255)                                                                                                             gTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCgg                               >  1:528600/1‑151 (MQ=255)                                                                                                             gTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCgg                               >  2:1041043/1‑151 (MQ=255)                                                                                                             gTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCgg                               >  2:1127082/1‑151 (MQ=255)                                                                                                               cACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGa                              <  1:414069/150‑1 (MQ=255)                                                                                                                aCCGATTTTCTGACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGt                            >  2:871151/1‑151 (MQ=255)                                                                                                                           gACCATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGt                  <  2:171137/150‑1 (MQ=255)                                                                                                                             ccATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATg               <  2:298999/151‑1 (MQ=255)                                                                                                                             ccATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATg               <  2:289897/151‑1 (MQ=255)                                                                                                                              cATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATg               >  2:458329/1‑150 (MQ=255)                                                                                                                              cATCGATCGCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATg               >  2:420503/1‑150 (MQ=255)                                                                                                                                 cgatcgCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACgg           >  1:383659/1‑151 (MQ=255)                                                                                                                                 cgatcgCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACg            <  2:430619/150‑1 (MQ=255)                                                                                                                                  gatcgCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACgg           >  1:181013/1‑150 (MQ=255)                                                                                                                                  gatcgCGGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACgg           >  2:1048165/1‑150 (MQ=255)                                                                                                                                  gatcgCGGCGGCAGGGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACgg           >  2:638657/1‑150 (MQ=255)                                                                                                                                     cgcgGCGGCAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAg        <  2:323429/150‑1 (MQ=255)                                                                                                                                          cggcAGTGACCACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGAc  <  2:129700/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             GACGACAAACCAGACCTGTCAACGCTGCGGGCGTTTATCACCAACGGCGAACAGGGAGTGAATTACCATCGTAACGTCTGGCATCACCCACTTTTCGCCTGGCAGCGCGTCACCGATTTTCTGACCATCGATCGCGGCGGCAGTGACAACTGTGATGTTGAAAGTATTCCTGAACAGGAACTCTGTTTTGCGTGACGCCTGCAACCGACTTGCATAAGATAAACTAATTGTTCATTGTTTATGCTCACTTGTAGGTCGGAGTTAACGTAGGTATGACGGAAGTTAGAC  >  NC_000913/532739‑533026

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