Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	217,044	A→G	7.1%	intergenic (‑41/+13)	yaeF ← / ← proS	putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	217,044	0	A	G	7.1%	313.7 / 10.9	98	intergenic (‑41/+13)	yaeF/proS	putative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base G (3/4);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

TACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCGG  >  NC_000913/216899‑217191                                                                                                                                                  |                                                                                                                                                    tACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGcc                                                                                                                                                   <  1:577833/148‑1 (MQ=255)     gctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttt                                                                                                                                               >  2:769751/1‑148 (MQ=255)      ctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACACCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                                                                             <  4:872200/149‑1 (MQ=255)      ctgTGGCTGATGGGTCGGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                                                                                              <  2:540900/148‑1 (MQ=255)       tgtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttt                                                                                                                                             <  1:769751/148‑3 (MQ=255)       tgtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                                                                            <  3:500637/149‑1 (MQ=255)           gCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCAt                                                                                                                                        <  3:289802/149‑1 (MQ=255)             tGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCa                                                                                                                                      >  3:420030/1‑149 (MQ=255)              gATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCa                                                                                                                                      >  1:595207/1‑148 (MQ=255)                     ctggctggCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTa                                                                                                                               <  2:224071/148‑1 (MQ=255)                       ggctggctGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                                                                                             <  2:880517/148‑1 (MQ=255)                           ggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                                         <  1:215602/148‑1 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                                          <  2:366426/143‑1 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                                          >  1:366426/1‑143 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                                     >  2:317922/1‑148 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                                     >  1:644774/1‑148 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                                     >  1:598605/1‑148 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                                     >  2:524364/1‑148 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                                    <  3:839327/149‑1 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATATAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                                          >  1:754060/1‑143 (MQ=255)                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATATAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                                          <  2:754060/143‑1 (MQ=255)                                    cAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                                                                               >  3:761692/1‑149 (MQ=255)                                           aCAAGCCGAAAGTATAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCACCGAAATTCAGGACGCTTTCTGGCAGCAGCCTTTAATCTGTTTGACCAGATATTcga                                                                                                         >  2:704182/1‑148 (MQ=255)                                            cAAGCCGAAAGTAAAAAAAAGTTGGGGAAAAACAGGCGGAAGAACCCCTTTGG‑TTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                      <  4:17468/149‑1 (MQ=255)                                              aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                      >  2:479949/1‑148 (MQ=255)                                                          aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtt                                                                                         >  3:632038/1‑149 (MQ=255)                                                           agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtt                                                                                         >  1:251349/1‑148 (MQ=255)                                                           agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtt                                                                                         >  2:460072/1‑148 (MQ=255)                                                            gaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                                                                       <  4:420030/149‑1 (MQ=255)                                                             aaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAGCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                                                                      <  4:805017/149‑1 (MQ=255)                                                                gCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatta                                                                                    <  1:479949/148‑1 (MQ=255)                                                                 cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                                   <  2:575084/148‑1 (MQ=255)                                                                 cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                                   <  1:460072/148‑1 (MQ=255)                                                                  tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                                                                  <  2:598605/148‑1 (MQ=255)                                                                   ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACt                                                                                 <  2:644774/148‑1 (MQ=255)                                                                   ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACt                                                                                 >  1:353921/1‑148 (MQ=255)                                                                   ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACt                                                                                 <  1:831884/148‑1 (MQ=255)                                                                     ggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATATCACCAGTCTTAATTAACTGt                                                                               <  1:552842/148‑1 (MQ=255)                                                                        gaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                                                                            <  2:162941/148‑1 (MQ=255)                                                                           aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                                                                        >  3:485589/1‑149 (MQ=255)                                                                            cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                                                                        >  1:180965/1‑148 (MQ=255)                                                                              ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGt                                                                     <  4:761692/149‑1 (MQ=255)                                                                                       aCGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGat                                                            <  4:637075/149‑1 (MQ=255)                                                                                           cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAttt                                                        >  4:38635/1‑149 (MQ=255)                                                                                              tGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtata                                                     <  4:632038/149‑1 (MQ=255)                                                                                                   tATCCATTTTATACAATCGATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTCCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCg                                                 <  1:704182/148‑1 (MQ=255)                                                                                                     tCCATTTTATACAATCCATGTAAAAAAAGGGCCATGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                                                                         <  3:674279/90‑1 (MQ=255)                                                                                                     tCCATTTTATACAATCCATGTAAAAAAAGGGCCATGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                                                                         >  4:674279/1‑90 (MQ=255)                                                                                                        aTTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtc                                            >  1:161851/1‑148 (MQ=255)                                                                                                        aTTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtc                                            >  2:356296/1‑148 (MQ=255)                                                                                                         ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgt                                          >  4:588454/1‑149 (MQ=255)                                                                                                         ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgt                                          >  4:270949/1‑149 (MQ=255)                                                                                                         ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgt                                          >  3:317265/1‑149 (MQ=255)                                                                                                         ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcg                                           >  1:659495/1‑148 (MQ=255)                                                                                                          tttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                                         <  4:627811/149‑1 (MQ=255)                                                                                                             ataCAAGCCATGGAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATcgtcgt                                       <  1:524364/148‑1 (MQ=255)                                                                                                              taCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGtt                                      >  1:88796/1‑148 (MQ=255)                                                                                                               aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                                         >  2:183357/1‑144 (MQ=255)                                                                                                               aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                                         <  1:183357/144‑1 (MQ=255)                                                                                                                 aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTc                                   >  3:552424/1‑148 (MQ=255)                                                                                                                 aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTc                                   <  4:552424/148‑1 (MQ=255)                                                                                                                   tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAg                                <  3:38635/149‑1 (MQ=255)                                                                                                                       tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                                <  2:145667/97‑1 (MQ=255)                                                                                                                       tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                                >  1:145667/1‑97 (MQ=255)                                                                                                                       tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCCGGCATCAGCCTTTAATCTGTGTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTa                            >  4:545875/1‑149 (MQ=255)                                                                                                                         tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTAc                           >  2:634705/1‑148 (MQ=255)                                                                                                                         tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg                          >  4:667853/1‑149 (MQ=255)                                                                                                                          aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg                          <  2:596734/148‑1 (MQ=255)                                                                                                                          aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg                          <  2:251349/148‑1 (MQ=255)                                                                                                                           aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACgg                         <  2:595207/148‑1 (MQ=255)                                                                                                                             aaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCg                      <  3:813065/149‑1 (MQ=255)                                                                                                                              aaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg                          >  2:711292/1‑144 (MQ=255)                                                                                                                              aaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg                          <  1:711292/144‑1 (MQ=255)                                                                                                                              aaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCg                      >  1:834855/1‑148 (MQ=255)                                                                                                                               aaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGcc                    <  3:262673/149‑1 (MQ=255)                                                                                                                               aaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGc                     >  1:61262/1‑148 (MQ=255)                                                                                                                                   gtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCg                                  <  1:351334/129‑1 (MQ=255)                                                                                                                                   gtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCg                                  >  2:351334/3‑131 (MQ=255)                                                                                                                                      cTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcac              >  2:335552/1‑148 (MQ=255)                                                                                                                                       tGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcaca             >  1:153756/1‑148 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                                   >  2:672469/1‑77 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                                   <  1:672469/77‑1 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACgg                         >  2:752354/1‑135 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACgg                         <  1:752354/135‑1 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACaa            >  1:842673/1‑148 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGATTCTCGCCGTTGCGACGATATTTATAtt                                                   <  2:424156/109‑1 (MQ=255)                                                                                                                                        gAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGATTCTCGCCGTTGCGACGATATTTATAtt                                                   >  1:424156/1‑109 (MQ=255)                                                                                                                                         aaaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAAt           <  2:353921/148‑1 (MQ=255)                                                                                                                                         aaaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAAt           <  1:335552/148‑1 (MQ=255)                                                                                                                                         aaaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAAt           <  1:317922/148‑1 (MQ=255)                                                                                                                                          atttCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAg         >  4:333820/3‑149 (MQ=255)                                                                                                                                           aTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAgt        >  4:849175/1‑149 (MQ=255)                                                                                                                                               aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGtgtg     <  1:356296/148‑1 (MQ=255)                                                                                                                                               aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGc    >  3:250442/1‑149 (MQ=255)                                                                                                                                                 gACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtat                                                      <  4:513988/97‑1 (MQ=255)                                                                                                                                                 gACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCgg  >  4:135119/1‑149 (MQ=255)                                                                                                                                                 gACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCg   <  2:180965/148‑1 (MQ=255)                                                                                                                                                 gACCCTTTCTGACATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtat                                                      >  3:513988/1‑97 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCGG  >  NC_000913/216899‑217191

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: ‑