Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,786,746	Δ1 bp	pseudogene (350/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,786,746	0	C	.	100.0%	391.5 / NA	92	pseudogene (350/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (45/47);  total (45/47)

TTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAC  >  NC_000913/2786616‑2786873                                                                                                                                   |                                                                                                                                 ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                                         >  1:469733/1‑139 (MQ=255) ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                                         >  1:840827/1‑139 (MQ=255)   aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                                       <  2:474527/139‑1 (MQ=255)   aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                                       <  1:375086/139‑1 (MQ=255)                aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  2:469733/139‑1 (MQ=255)                aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  1:1416368/139‑1 (MQ=255)                     ttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                                     <  1:459729/139‑1 (MQ=255)                      tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtgt                                                                                                    >  1:1005127/1‑139 (MQ=255)                       cccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                                   <  2:840827/139‑1 (MQ=255)                        ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  <  2:657853/139‑1 (MQ=255)                           ttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATc                                                                                               <  2:628845/139‑1 (MQ=255)                                  gATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                        <  2:1005127/139‑1 (MQ=255)                                  gATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                        <  2:1446883/139‑1 (MQ=255)                                      cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGt                                                                                    <  1:1028140/139‑1 (MQ=255)                                       aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                                   <  1:721830/139‑1 (MQ=255)                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                                <  1:248106/121‑1 (MQ=255)                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                                >  2:248106/1‑121 (MQ=255)                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                              >  1:629885/1‑139 (MQ=255)                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                              >  1:1073707/1‑139 (MQ=255)                                            aCCAAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  >  2:640634/1‑119 (MQ=255)                                            aCCAAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  <  1:640634/119‑1 (MQ=255)                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                                        <  2:1159226/112‑1 (MQ=255)                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                                        >  1:1159226/1‑112 (MQ=255)                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  <  1:784931/114‑1 (MQ=255)                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  >  2:784931/1‑114 (MQ=255)                                                 gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                         >  2:569306/1‑139 (MQ=255)                                                    ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                      <  1:520864/139‑1 (MQ=255)                                                    ttGAATAAAGAAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATATCag                                                                      >  1:64771/1‑138 (MQ=255)                                                     tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                     <  2:981276/139‑1 (MQ=255)                                                         tAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                                 <  2:1441223/139‑1 (MQ=255)                                                                 cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt                                                         >  2:801367/1‑139 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          >  2:1180328/1‑85 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  1:1180328/85‑1 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                                    >  1:909100/1‑139 (MQ=255)                                                                           cATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCtt                                               >  2:609411/1‑139 (MQ=255)                                                                           cATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCtt                                               >  2:853754/1‑139 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAaga                                                                                                                     <  1:520294/64‑1 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAaga                                                                                                                     >  2:520294/1‑64 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  2:1026811/75‑1 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          >  1:1026811/1‑75 (MQ=255)                                                                                  gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCaaa                                       >  2:119120/1‑139 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  >  1:732601/1‑80 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACAtt                                                               <  2:1372857/115‑1 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACAtt                                                               >  1:1372857/1‑115 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCaaaa                                      <  1:609411/139‑1 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTATATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  <  2:732601/80‑1 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGATTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCTCATTATATCTGACAAACCTT‑CCCCaaaa                                      >  2:1298312/1‑139 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  1:1047802/70‑1 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          >  2:1047802/1‑70 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGa                                    <  1:593727/139‑1 (MQ=255)                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATaa                                 >  2:828856/1‑139 (MQ=255)                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATaa                                 <  2:1064659/139‑1 (MQ=255)                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTA‑CCCCAAAAGATaa                                 <  2:629885/139‑1 (MQ=255)                                                                                            ggATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAatat                             <  2:1073707/139‑1 (MQ=255)                                                                                             gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCCAAAAGATAAat                               <  2:710196/137‑1 (MQ=255)                                                                                             gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCCAAAAGATAAat                               >  1:710196/1‑137 (MQ=255)                                                                                             gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATg                            >  1:808910/1‑139 (MQ=255)                                                                                               tAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAata                              <  2:655321/135‑1 (MQ=255)                                                                                               tAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAata                              >  1:655321/1‑135 (MQ=255)                                                                                                  gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                        >  2:391815/1‑75 (MQ=255)                                                                                                  gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                                        <  1:391815/75‑1 (MQ=255)                                                                                                    gaaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCata                     <  1:512664/139‑1 (MQ=255)                                                                                                    gaaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCata                     <  1:119120/139‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          >  2:136966/1‑52 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  1:670264/52‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          >  2:670264/1‑52 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                                          <  1:136966/52‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                         >  1:1024406/1‑85 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                                         <  2:1024406/85‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAt                                                                           <  2:284287/83‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAt                                                                           >  1:284287/1‑83 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAatat                             <  2:13820/128‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAatat                             >  1:13820/1‑128 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCat                      <  2:458776/135‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCat                      >  1:458776/1‑135 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGc                  >  1:586369/1‑139 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGc                  >  2:873801/1‑139 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGc                  >  1:285673/1‑139 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGc                  >  2:1375146/1‑139 (MQ=255)                                                                                                        gaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCt                 <  2:909100/139‑1 (MQ=255)                                                                                                          aaaaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCTAt               <  2:586369/139‑1 (MQ=255)                                                                                                                   gTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCCT‑CCCCAAAAGATAAATATGTTGCATATGCTATTAACATTcc      >  2:1194601/1‑139 (MQ=255)                                                                                                                     ttATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTg    >  2:1282173/1‑139 (MQ=255)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  >  1:1265170/1‑44 (MQ=38)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                  <  2:1265170/44‑1 (MQ=255)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCata                     <  2:686418/120‑1 (MQ=255)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCata                     >  1:686418/1‑120 (MQ=255)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAc  >  2:798747/1‑139 (MQ=255)                                                                                                                         cAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCaaa                                       <  2:64031/100‑1 (MQ=255)                                                                                                                         cAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCaaa                                       >  1:64031/1‑100 (MQ=255)                                                                                                                             gagaTA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGt                           >  2:1008940/1‑108 (MQ=255)                                                                                                                             gagaTA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGt                           <  1:1008940/108‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 TTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTT‑CCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAC  >  NC_000913/2786616‑2786873

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