Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,547,395	A→G	100%	intergenic (‑226/‑6)	yddG ← / → fdnG	aromatic amino acid exporter/formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,547,395	0	A	G	99.0%	321.6 / ‑4.7	97	intergenic (‑226/‑6)	yddG/fdnG	aromatic amino acid exporter/formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (48/48);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

TTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGAAGCAA  >  NC_000913/1547300‑1547490                                                                                                |                                                                                                tttGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGc                                                                                             >  1:2816424/1‑100 (MQ=255)  ttGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCa                                                                                            >  2:2320520/1‑100 (MQ=255)  tccgatctaaCAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCTCAGTAATACCCCTGAAAAAAGAGGGAAGCa                                                                                            <  1:71207/92‑1 (MQ=255)      ggCAAACAGATGTTCTTTTTGATTTAGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATgg                                                                                        >  2:527780/1‑100 (MQ=255)        cAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATg                                                                                         >  2:1244569/1‑97 (MQ=255)        cAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGAc                                                                                      >  1:1966201/1‑100 (MQ=255)        cAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGAc                                                                                      >  2:2999435/1‑100 (MQ=255)        cAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGAc                                                                                      >  1:1671643/1‑100 (MQ=255)         aaaTAGATGTTCTTTTTGATTTCGCGCAAAAAGATTTAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACg                                                                                     <  1:1108987/100‑1 (MQ=255)          aaCAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCTCAGTAATACCCCTGAAAAAAGAGGGAAGCaagatcgga                                                                                    >  2:71207/1‑93 (MQ=255)           aCAGATGTTCTTTTTGATTTCCCGCAAACAGATTCAGAAGTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACgac                                                                                   <  2:1676531/100‑3 (MQ=255)              gATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGt                                                                                <  2:372802/100‑1 (MQ=255)                tGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTTAAAAAAGAGGGAAGCAATGGACGTCAGTCg                                                                              <  2:189778/100‑1 (MQ=255)                 gttgttTTTGATTTCGCACAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGc                                                                             <  2:1540477/100‑1 (MQ=255)                 gTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGc                                                                             >  1:294688/1‑100 (MQ=255)                  ttCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTGTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCa                                                                            <  2:1239283/100‑1 (MQ=255)                   tCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAg                                                                           >  1:2489332/1‑100 (MQ=255)                     tttttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGAc                                                                         >  2:795479/1‑100 (MQ=255)                     tttttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGAc                                                                         <  1:2183015/100‑1 (MQ=255)                       tttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACaa                                                                       >  1:1476842/1‑100 (MQ=255)                       tttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACaa                                                                       <  2:2518865/100‑1 (MQ=255)                        ttGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAt                                                                      <  1:2999468/100‑1 (MQ=255)                         tGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAtt                                                                     >  1:729599/1‑100 (MQ=255)                          gATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCATGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAAttt                                                                    <  1:1637740/100‑1 (MQ=255)                            tttCCCGCAAAAAGTTTCAGAATTTTACTGTTAGTTTCCTCGAGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGTCGTCAGTCGCAGACAAttttt                                                                  <  2:69633/100‑1 (MQ=255)                              tCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTa                                                                <  1:1228187/100‑1 (MQ=255)                               cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  2:2067542/100‑1 (MQ=255)                               cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  2:861278/100‑1 (MQ=255)                               cgcgCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  2:3013412/100‑1 (MQ=255)                                gcgcAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaaa                                                              <  1:330252/100‑1 (MQ=255)                                 cgcAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaaaa                                                             >  2:2218562/1‑100 (MQ=255)                                  gcAAAAAGATTCAGAATTTTATTGTTAGTTTCCTCGCTCAGTAATACAAATAAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTaa                                                               <  2:932498/97‑1 (MQ=255)                                   cAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATc                                                           >  2:404445/1‑100 (MQ=255)                                     aaaaGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTg                                                         <  2:1084200/100‑1 (MQ=255)                                       aaGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATTCCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCtn                                                         <  2:1082423/98‑2 (MQ=255)                                       aaGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTgcg                                                       <  2:164544/100‑1 (MQ=255)                                           ttCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcggg                                                   <  2:889369/100‑1 (MQ=255)                                           ttCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcgg                                                    >  1:518933/1‑99 (MQ=255)                                             cAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTTAAATCTGCgcgggcg                                                 >  2:906448/1‑100 (MQ=255)                                             cAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcgggcg                                                 >  2:2348016/1‑100 (MQ=255)                                             cAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcgggcg                                                 <  1:732839/100‑1 (MQ=255)                                              aGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCgcgggcgg                                                <  2:2485863/100‑1 (MQ=255)                                               gAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGt                                               >  2:2036062/1‑100 (MQ=255)                                                aaTTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTa                                              <  1:2036085/100‑1 (MQ=255)                                                aaTTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTa                                              <  1:903828/100‑1 (MQ=255)                                                  ttttACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATg                                            >  2:2441901/1‑100 (MQ=255)                                                   tttACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTAATTAAAATCTGCGCGGGCGGTATgg                                           >  1:2571361/1‑100 (MQ=255)                                                    ttACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGc                                          <  1:1731413/100‑1 (MQ=255)                                                    ttACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAACAATGGACGTCAATCGCAGACAATTTTTTAAAATCTGCGCGGTCGGTa                                              >  1:1890965/1‑96 (MQ=255)                                                      aCTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTg                                        <  1:1443021/100‑1 (MQ=255)                                                      aCTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTg                                        <  2:1602413/100‑1 (MQ=255)                                                       cTGTTAGTTTCCTCGCTCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTgg                                       <  2:2395690/100‑1 (MQ=255)                                                       cTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTgg                                       <  2:880163/100‑1 (MQ=255)                                                       cTGTTAGTTTCCTCGCGCACTAATACCCCTGAAAAAACAAGGAATCAATGGACGTCAGTCGCAGCCAATTTTTTAACATCTGCGCGGGCGGTATGGCTgg                                       <  1:2424721/100‑1 (MQ=255)                                                        tcttAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGa                                      <  1:2640204/98‑1 (MQ=255)                                                        tGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGa                                      <  1:350727/100‑1 (MQ=255)                                                        tGTTAGTTTCCTCGCACAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGa                                      >  2:601691/1‑100 (MQ=255)                                                         gTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaa                                     >  2:631476/1‑100 (MQ=255)                                                         gTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACATCATTCGCAGACAA‑TTTTTAAAATCTGCGCGGGCGGTATGGCTGGaa                                     >  2:905136/1‑99 (MQ=255)                                                          ttAGTTTCCTCGCGCAGTAATACCTCTGAAAAAAGGGGGAAGCTATGGACGTCTGTCGCAGACAATTTTTTAATATCTGCGCGGGCGGTATGGCTGGaac                                    <  2:1192614/100‑1 (MQ=255)                                                          ttAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaac                                    <  1:557161/100‑1 (MQ=255)                                                           tAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaaga                                   >  2:2640172/1‑98 (MQ=255)                                                            aGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaacaa                                  >  1:2149950/1‑100 (MQ=255)                                                             ntttCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAAATTTTTAAAATCTACGCGGGCGGCATGGCTGGaacaac                                 >  1:2453055/2‑100 (MQ=255)                                                             gTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaacaac                                 <  2:2448109/100‑1 (MQ=255)                                                             gTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGaacaac                                 >  1:297897/1‑100 (MQ=255)                                                                 ccTCTCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTa                             >  2:1282696/1‑100 (MQ=255)                                                                 ccTCGCTCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTa                             >  1:2140483/1‑100 (MQ=255)                                                                  cTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTa                             >  2:1564141/1‑99 (MQ=255)                                                                  cTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAg                            >  2:2530386/1‑100 (MQ=255)                                                                  cTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGAGCGGGCGGTATGGCTGGAACAACAGTAg                            >  1:175040/1‑100 (MQ=255)                                                                  cTCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGAAATGGACGTCAGTCGCAGACAACTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTa                             >  2:2585077/1‑99 (MQ=255)                                                                   tCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGc                           >  2:627540/1‑100 (MQ=255)                                                                   tCGCGCAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCAGGCGGTATGGCTAGAACAACAGTAGc                           >  1:2634633/1‑100 (MQ=255)                                                                     gcgcAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCgg                         <  2:80716/100‑1 (MQ=255)                                                                      cgcAGTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGc                        >  2:1879944/1‑100 (MQ=255)                                                                          gTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTg                    <  2:1613115/100‑1 (MQ=255)                                                                          gTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTg                    >  2:1520040/1‑100 (MQ=255)                                                                          gTAATACCCCTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTg                    >  2:1545111/1‑100 (MQ=255)                                                                                 cccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTg             <  1:1310203/100‑1 (MQ=255)                                                                                 cccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTg             >  1:1638014/1‑100 (MQ=255)                                                                                 cccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTg             >  1:1207836/1‑100 (MQ=255)                                                                                 cccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGATGGAACAACAGTAGCGGCATTAGGCTTTg             >  1:674408/1‑100 (MQ=255)                                                                                  ccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGc            <  2:297894/100‑1 (MQ=255)                                                                                  ccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGc            >  1:2226790/1‑100 (MQ=255)                                                                                  ccTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGAAGCGGCATTGGGCTTTGc            >  1:381154/1‑100 (MQ=255)                                                                                   cTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcc           <  1:161891/100‑1 (MQ=255)                                                                                   cTGAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTGTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcc           >  1:1356720/1‑100 (MQ=255)                                                                                    tGAAAAAAGAGTGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGccc          <  2:733512/100‑1 (MQ=255)                                                                                     gAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcccc         <  2:1671623/100‑1 (MQ=255)                                                                                     gAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcccc         <  1:837570/100‑1 (MQ=255)                                                                                     gAAAAAAGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGcccc         <  1:285063/100‑1 (MQ=255)                                                                                     gAAAAAAGAGGGAAGCAATGGAAGTCAGTCGCAGACAATTTTTTAAAATCTGCGAGGGCGGTATGGCTGGACCAACAGTAGCGGCATTGGGCTTTGccc          >  2:350165/1‑99 (MQ=255)                                                                                        aaaaGATGGAATCAATTGACGTCAGTCGTAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGa       <  1:596118/99‑1 (MQ=255)                                                                                         aaaGAGGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGaag     <  1:2325900/100‑1 (MQ=255)                                                                                           agagGGAAGCAATGGACGTCAGTCGCAGACAATTTTGTAAAAGCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCTGCATTGGGCTTTGCCTCGgagc    <  2:1804912/99‑1 (MQ=255)                                                                                            gagGGAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGaagcaa  >  1:466242/1‑100 (MQ=255)                                                                                                |                                                                                                TTTGTGGCAAACAGATGTTCTTTTTGATTTCGCGCAAAAAGATTCAGAATTTTACTGTTAGTTTCCTCGCGCAGTAATACCCCTGAAAAAAGAGGAAAGCAATGGACGTCAGTCGCAGACAATTTTTTAAAATCTGCGCGGGCGGTATGGCTGGAACAACAGTAGCGGCATTGGGCTTTGCCCCGAAGCAA  >  NC_000913/1547300‑1547490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑