Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,181,760	A→G	5.8%	intergenic (‑35/+23)	ymfA ← / ← potD	DUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,181,760	0	A	G	5.8%	208.9 / 2.7	69	intergenic (‑35/+23)	ymfA/potD	DUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base G (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACT  >  NC_000913/1181613‑1181909                                                                                                                                                    |                                                                                                                                                      aTCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAAt                                                                                                                                                     >  1:455330/1‑150 (MQ=255)  tCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATa                                                                                                                                                    >  2:172534/1‑150 (MQ=255)      gTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccac                                                                                                                                                <  1:1059112/150‑1 (MQ=255)      gTCTTTGCTATAGGTGAACTGGGTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccac                                                                                                                                                <  1:1064468/150‑1 (MQ=255)       tCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAACCGGGCGGTAAAaccacc                                                                                                                                               >  1:975287/1‑150 (MQ=255)       tCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccacc                                                                                                                                               <  1:111272/150‑1 (MQ=255)       tCTATGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATaccacc                                                                                                                                               <  2:425380/150‑1 (MQ=255)          ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcc                                                                                                                                            <  2:603668/150‑1 (MQ=255)          ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcc                                                                                                                                            >  1:1184160/1‑150 (MQ=255)          ttGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACAAccgcc                                                                                                                                            >  1:308836/1‑150 (MQ=255)           tGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACAGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgccc                                                                                                                                           <  2:198317/150‑1 (MQ=255)             cTATAGGTGAATTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCAccgcccgc                                                                                                                                         >  2:1128555/1‑150 (MQ=255)             cTATAGGTGAACGGGTAATAAAAGACAACAACATAAATGAGGATCAGTACCAGACCGATACCGAGGAATAAACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTGGTATTAccgcccgt                                                                                                                                         <  1:958275/150‑2 (MQ=255)                taGGGGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTg                                                                                                                                      <  1:136241/150‑1 (MQ=255)                 aGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGc                                                                                                                                     <  2:961552/150‑1 (MQ=255)                 aGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGACGAATATACGAAAGAAAACTTTACTGGCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGc                                                                                                                                     <  1:56359/150‑1 (MQ=255)                  ggTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAGAAACGGGCGGTGGTATTACCGCCCGtttttt                                                                                                                                    <  1:339921/150‑4 (MQ=255)                   gTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTg                                                                                                                                   <  1:552398/150‑1 (MQ=255)                     gAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGaa                                                                                                                                 >  2:5854/1‑150 (MQ=255)                           gTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACg                                                                                                                           <  2:1133924/150‑1 (MQ=255)                            ttATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGt                                                                                                                          <  1:295435/150‑1 (MQ=255)                             tataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTc                                                                                                                         >  2:323888/1‑150 (MQ=255)                              ataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCTATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTcc                                                                                                                        <  1:956627/150‑1 (MQ=255)                              ataAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTACCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGAcc                                                                                                                        <  1:956630/150‑1 (MQ=255)                               taAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCt                                                                                                                       >  2:70678/1‑150 (MQ=255)                                     acaacaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCTCCGCCCGCTTGCTGAATTAACGTCCTGCTTTc                                                                                                                 <  1:1195050/150‑1 (MQ=255)                                      caacaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCa                                                                                                                <  1:570818/150‑1 (MQ=255)                                       aacaacaGAAATTAGTATCAGTACCAGACCTATACCGAGGAATATACGAAAGTAAACTTTACTGTCCTTAATCCTTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAg                                                                                                               <  2:362251/150‑1 (MQ=255)                                        acaacaGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGc                                                                                                              <  1:1061712/150‑1 (MQ=255)                                               aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1231064/150‑1 (MQ=255)                                               aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1162932/150‑1 (MQ=255)                                               aaaTGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGat                                                                                                       <  2:1157737/150‑1 (MQ=255)                                                   taGGATCAGTACCAGACCTATTCCGAGGAATATACGAAAGTAAACTTAACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACTTGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                   <  1:78696/149‑1 (MQ=255)                                                      gATCAGTACCAGACCGATACAGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCCTGCTAACTAAAAAACGGGCGTTAATAACAACGCCCGCTTGCTCAATAAACGTCATACTTACAGCTTCGGAAAATActg                                                                                                >  2:819184/1‑149 (MQ=255)                                                        tCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCtt                                                                                              >  1:240146/1‑150 (MQ=255)                                                        tCAGTACCAGACCGATACCGAGGAATAAACTAAAAAAAACTTTAAAGTAATCACTCATTAACCACATCCTAGATAACTAAAAAACGGGCGGAAAAACAACCGCCCGCTTGCTGAATTAACGTACAGATATCACCATCTGATAATAcacat                                                                                              >  2:285760/1‑146 (MQ=255)                                                         cAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTc                                                                                             >  2:944652/1‑150 (MQ=255)                                                          aGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCa                                                                                            <  1:377820/150‑1 (MQ=255)                                                                aGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAAtg                                                                                      <  1:167460/150‑1 (MQ=255)                                                                           gAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTgc                                                                           >  1:544078/1‑150 (MQ=255)                                                                               aaTATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCTCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCGTCTGATAATTCTCTTCATAAATGCTGCAGGCTGCGccc                                                                       >  2:857539/1‑149 (MQ=255)                                                                                  ataCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACg                                                                    <  2:929846/150‑1 (MQ=255)                                                                                   taCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGt                                                                   >  1:1057936/1‑150 (MQ=255)                                                                                      gaaagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCAAACGTCAt                                                                <  2:1132971/150‑1 (MQ=255)                                                                                          gaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAATAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTg                                                            >  1:941085/1‑150 (MQ=255)                                                                                                   aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:134712/150‑1 (MQ=255)                                                                                                   aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:934257/150‑1 (MQ=255)                                                                                                   aCTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGcc                                                   <  2:31732/150‑1 (MQ=255)                                                                                                        ccTGACTCATTAGACACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGATTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCTTTCTGCAATTCGAAAtttt                                              >  2:818207/1‑150 (MQ=255)                                                                                                          tGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTa                                            >  1:609177/1‑150 (MQ=255)                                                                                                             cTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATg                                         <  2:478618/150‑1 (MQ=255)                                                                                                                aTTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGtt                                      >  2:718995/1‑150 (MQ=255)                                                                                                                aTTAGCCACATCCTAGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTTCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGtt                                      <  1:962616/150‑1 (MQ=255)                                                                                                                    gCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAGTTTTAATGGTTTCAg                                  >  2:675089/1‑150 (MQ=255)                                                                                                                     ccACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGc                                 <  1:121176/150‑1 (MQ=255)                                                                                                                      cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                                >  1:924551/1‑150 (MQ=255)                                                                                                                        caTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATc                              >  1:582445/1‑150 (MQ=255)                                                                                                                            cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                          >  1:353079/1‑150 (MQ=255)                                                                                                                            cTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                          >  1:358091/1‑150 (MQ=255)                                                                                                                                  aactaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgt                    <  1:323888/150‑1 (MQ=255)                                                                                                                                    ctaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAAGTAACGTCCTGCTTTCAGCTTCTGCTAATACTCTTCATAAGTGCTGCTGGCTGCGCCAACGTCATTCTGCCAATCGCCATTTTTTATGTTTTCAGCTTCCGGGTAGAgagt                  >  2:792349/1‑150 (MQ=255)                                                                                                                                         aaaCGGGAGGTAATACGACCGCCCGCTTGCTGAATTACCGTCCTGCATTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTAAt             <  1:500292/150‑1 (MQ=255)                                                                                                                                           gcGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATGCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCg           >  2:95065/2‑150 (MQ=255)                                                                                                                                            cGGGCGGTGGTATTACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGt          >  2:727260/1‑150 (MQ=255)                                                                                                                                            cGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGTGTGTTTTATCGt          <  1:383830/150‑1 (MQ=255)                                                                                                                                            cGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGt          <  1:423138/150‑1 (MQ=255)                                                                                                                                             gggCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATAGTCTTCATAAATGCTTCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTATTGGTTTCAGCATCCGGGTAGAGTGTTTTATCGtt         >  2:650873/1‑150 (MQ=255)                                                                                                                                             gggCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGtt         >  1:1696/1‑150 (MQ=255)                                                                                                                                                 ggTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc     >  1:135325/1‑150 (MQ=255)                                                                                                                                                    aaTACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACt  <  2:66287/150‑1 (MQ=255)                                                                                                                                                    aaTACCACCGCCCGCTTGCTGAATTAACGTCCTGCTATCAGCTTCTGATAATACTCTACATAAATGCTGCTGGCTGCGCCAACGACATACTGCCAATCCCCAGTTTTAATGGATTACGCATCCGGGTAGAGTGTTTTATCGTTCGCCACt  >  1:1203899/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ATCGCGTCTTTGCTATAGGTGAACTGGTTATAAAAGACAACAACAGAAATGAGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACT  >  NC_000913/1181613‑1181909

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