Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,192	G→A	Q447Q (CAG→CAA)	ydgA →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,192	0	G	A	100.0%	305.8 / NA	93	Q447Q (CAG→CAA)	ydgA	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/51);  total (42/51)

GGATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCC  >  NC_000913/1691057‑1691317                                                                                                                                        |                                                                                                                              ggATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATg                                                                                                                            >  1:561193/1‑139 (MQ=255)     gCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTcc                                                                                                                        >  2:1123703/1‑139 (MQ=255)     gCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTcc                                                                                                                        >  2:177600/1‑139 (MQ=255)       gAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGt                                                                                                                      >  2:719904/1‑139 (MQ=255)              aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               >  1:731471/1‑139 (MQ=255)              aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  1:1361378/139‑1 (MQ=255)              aCCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  1:1634840/139‑1 (MQ=255)                 aTTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCt                                                                                                            <  2:731471/139‑1 (MQ=255)                   tCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                          <  2:202278/139‑1 (MQ=255)                       gTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  1:1099605/130‑1 (MQ=255)                       gTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               >  2:1099605/1‑130 (MQ=255)                       gTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAgg                                                                                                      <  1:177600/139‑1 (MQ=255)                             atGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                >  2:487403/1‑139 (MQ=255)                             atGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                >  1:47785/1‑139 (MQ=255)                             atGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                <  2:1760605/139‑1 (MQ=255)                               ggCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATAcc                                                                                              >  2:1765446/1‑139 (MQ=255)                                  aaCTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCACGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATc                                                                                           >  2:1347972/1‑139 (MQ=255)                                    cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                          <  2:1292121/122‑1 (MQ=255)                                    cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                          >  1:1292121/1‑122 (MQ=255)                                    cTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcac                                                                                         >  2:1707661/1‑139 (MQ=255)                                     tGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                                                            >  1:43495/1‑135 (MQ=255)                                     tGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                                                            <  2:43495/135‑1 (MQ=255)                                      gAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       <  1:530065/139‑1 (MQ=255)                                      gAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       <  2:47785/139‑1 (MQ=255)                                       aGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcaccac                                                                                      <  1:1123703/139‑1 (MQ=255)                                        gTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCaccacc                                                                                     <  2:561193/139‑1 (MQ=255)                                           tATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGc                                                                                  <  2:1110811/139‑1 (MQ=255)                                             tGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCt                                                                                <  2:267086/139‑1 (MQ=255)                                               aCTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGc                                                                              <  1:719904/139‑1 (MQ=255)                                                cTCAGGTAGCGAAGCTGGAATGTTATCAGTAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGTGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCa                                                                             <  1:487403/139‑1 (MQ=255)                                                 tCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCCAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                <  1:1115353/119‑1 (MQ=255)                                                 tCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCaa                                                                            <  2:871215/139‑1 (MQ=255)                                                 tCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGGTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                >  2:1115353/1‑119 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGt                                                                                                                      >  1:66825/1‑95 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGt                                                                                                                      <  2:66825/95‑1 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAata                                                                          >  2:839637/1‑139 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAata                                                                          >  2:768726/1‑139 (MQ=255)                                                   aGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAata                                                                          <  1:1804598/139‑1 (MQ=255)                                                     gTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       <  1:1664620/124‑1 (MQ=255)                                                     gTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       >  2:1664620/1‑124 (MQ=255)                                                        gCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCaa                                                                            <  1:1495086/132‑1 (MQ=255)                                                        gCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCaa                                                                            >  2:1495086/1‑132 (MQ=255)                                                         cGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCt                                                                      <  1:787984/137‑1 (MQ=255)                                                         cGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCt                                                                      >  2:787984/1‑137 (MQ=255)                                                         cGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTaa                                                                    <  1:744950/139‑1 (MQ=255)                                                         cGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTaa                                                                    <  1:1707661/139‑1 (MQ=255)                                                             gCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCCATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGcc                                                                                 >  1:820871/1‑122 (MQ=255)                                                             gCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGcc                                                                                 <  2:820871/122‑1 (MQ=255)                                                                 gAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGa                                                            <  1:768726/139‑1 (MQ=255)                                                                  aaGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                >  1:905446/1‑102 (MQ=255)                                                                  aaGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                <  2:905446/102‑1 (MQ=255)                                                                  aaGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGAt                                                           <  1:839637/139‑1 (MQ=255)                                                                  aaGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTCGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGAt                                                           <  1:1758953/139‑1 (MQ=255)                                                                        aTCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGTTAACGGTCAGATAACGtt                                                     >  2:1090516/1‑139 (MQ=255)                                                                               aGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               >  1:1635248/1‑74 (MQ=255)                                                                               aGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  2:1635248/74‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  2:964074/73‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               >  1:964074/1‑73 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                >  2:69766/1‑88 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATa                                                                                                <  1:69766/88‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGc                                                                                  >  1:1426992/1‑102 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGc                                                                                  <  2:1426992/102‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGt                                                                <  2:858469/120‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGt                                                                >  1:858469/1‑120 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTaaa                                                  >  2:561647/1‑134 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTaaa                                                  <  1:561647/134‑1 (MQ=255)                                                                                gATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGc                                             >  1:1764480/1‑139 (MQ=255)                                                                                    aaGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGt                                                                >  2:733612/1‑116 (MQ=255)                                                                                    aaGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGt                                                                <  1:733612/116‑1 (MQ=255)                                                                                     aGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGaaa                                        <  1:1765446/139‑1 (MQ=255)                                                                                       cGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                          >  2:125742/1‑71 (MQ=255)                                                                                       cGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTg                                                                                                          <  1:125742/71‑1 (MQ=255)                                                                                        gAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               >  2:1777669/1‑65 (MQ=255)                                                                                        gAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGacca                                                                                                               <  1:1777669/65‑1 (MQ=255)                                                                                          aaaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGAc                                                                                                    <  1:764635/74‑1 (MQ=255)                                                                                          aaaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGAc                                                                                                    >  2:764635/1‑74 (MQ=255)                                                                                           aaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCATATAACGTTAAATGGGCAGAAAATGTCg                                  >  2:1117708/1‑139 (MQ=255)                                                                                           aaaaCTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCg                                  >  2:838192/1‑139 (MQ=255)                                                                                                            aaGTTGAAGGTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTgg                 <  1:83130/139‑1 (MQ=255)                                                                                                                     gTGCATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGc        <  2:1351043/139‑1 (MQ=255)                                                                                                                        cATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGa                                                            <  2:1044209/84‑1 (MQ=255)                                                                                                                        cATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATgg                                               <  1:244716/97‑1 (MQ=255)                                                                                                                        cATCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATgg                                               >  2:244716/1‑97 (MQ=255)                                                                                                                        cATCAGCAATGGGGCAAACGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGa                                                            >  1:1044209/1‑84 (MQ=255)                                                                                                                         aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       <  1:180475/56‑1 (MQ=255)                                                                                                                         aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       >  1:191720/1‑56 (MQ=255)                                                                                                                         aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       <  2:191720/56‑1 (MQ=255)                                                                                                                         aTCAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATcacca                                                                                       >  2:180475/1‑56 (MQ=255)                                                                                                                           cAGCAATGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGcc  <  2:1764480/139‑1 (MQ=255)                                                                                                                               aaTGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGt                                    <  1:1406301/101‑1 (MQ=255)                                                                                                                               aaTGGGGCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGt                                    >  2:1406301/1‑101 (MQ=255)                                                                                                                                  ggggCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                                                            >  2:658317/1‑42 (MQ=255)                                                                                                                                  ggggCAAATGTTCCGTCTGACCACCTTGCAGGACAATACCAt                                                                                            <  1:658317/42‑1 (MQ=255)                                                                                                                                        |                                                                                                                              GGATGCGAAACTGACCATTCCGGTGGATATGGCAACTGAGTTTATGACTCAGGTAGCGAAGCTGGAAGGTTATCAGGAAGATCAAGCGAAAAAACTGGCGAAACAGCAAGTTGAAGGTGCATCAGCAATGGGGCAGATGTTCCGTCTGACCACCTTGCAGGACAATACCATCACCACCAGCCTGCAATATGCTAACGGTCAGATAACGTTAAATGGGCAGAAAATGTCGTTAGAAGATTTTGTTGGCATGTTTGCGATGCC  >  NC_000913/1691057‑1691317

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