Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,030	1	.	G	78.3%	85.1 / 24.7	46	intergenic (+201/+535)	ycdU/serX	putative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (9/1);  new base G (11/25);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGC  >  NC_000913/1096884‑1097171                                                                                                                                                        |                                                                                                                                                    aTTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑g‑g‑‑gggg                                                                                                                                                >  2:650776/1‑151 (MQ=255) aTTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑gggg                                                                                                                                                <  2:1160061/151‑1 (MQ=255)  ttttGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑ggggg                                                                                                                                               >  1:380418/1‑151 (MQ=255)  ttttGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑ggggg                                                                                                                                               >  2:573875/1‑151 (MQ=255)     tGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑ggggggg                                                                                                                                             <  1:95514/150‑1 (MQ=255)      gATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑ggggggggg                                                                                                                                           >  1:643566/1‑149 (MQ=255)      gATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑gggggggg                                                                                                                                            >  2:429381/1‑149 (MQ=255)      gATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑ggggggggat                                                                                                                                          <  2:179476/151‑4 (MQ=255)       aTTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGAGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑G‑C‑‑gggggggggg                                                                                                                                          >  1:1149/1‑148 (MQ=255)                  aGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACt‑cgcttc                                                                                                                               <  1:814476/151‑1 (MQ=255)                             aGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑G‑G‑‑GGGGGGGGTTGACT‑CGCTTCGCTCGccc                                                                                                                       >  1:191529/1‑148 (MQ=255)                              gTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGGTTGACT‑CCCTTCGCTCcccc                                                                                                                       >  2:244973/1‑147 (MQ=255)                              gTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCg                                                                                                                   <  2:736886/151‑1 (MQ=255)                              gTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCg                                                                                                                   <  1:964650/151‑1 (MQ=255)                                  aTGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑G‑AG‑GGGGGGGGGTGACT‑CGCCTCCCTCGCCCCGCggg                                                                                                                 >  1:507525/1‑149 (MQ=18)                                    gtgCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAg                                                                                                              <  2:961413/150‑1 (MQ=255)                                     tgCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGAGTGAAT‑CGCTTCGCCCGCCCCGCGGGGAGcc                                                                                                            >  2:514328/1‑151 (MQ=18)                                     tgCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑G‑AG‑GGGGGGGATTGACT‑CGCTTCTCCCGCCCTCCGGGGAGcc                                                                                                            >  2:964619/1‑151 (MQ=18)                                        aCGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGG‑G‑A‑‑GGGGGGGGTTTAGC‑CGCTTCGCCCCCCCTGCGGGGCGccc                                                                                                           >  1:314094/1‑150 (MQ=18)                                         cGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGAAT‑CGCCTCGCACGCCCTGCgggaccccc                                                                                                           >  2:605508/1‑142 (MQ=18)                                           aTTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCa                                                                                                      <  2:661561/151‑1 (MQ=255)                                             ttCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑G‑AG‑GGGGGGGGGTGACT‑CGCCTCCCTCCCCCCGCGGGGGGccccccc                                                                                                       >  1:1006525/1‑144 (MQ=17)                                               ccccTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGGTTGACT‑CGCTTCGCTCGCCCTGCGGGCAGccc                                                                                                           >  2:1156962/1‑142 (MQ=255)                                                   tcccccTAACAGACGATGGCGGGGTCACCCGTTGTGTTCTATATACG‑AAAAAAAAGCCCGTACTTTCGTTCTATCTCTTCTGTAACGACCAC‑CGG‑G‑A‑‑‑GGGCCGATTGACG‑CGATTCCCTCGCCCTCCCGGCACCCCCCTTACttgatt                                                                                               <  1:860393/145‑6 (MQ=9)                                                     gCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGGTTTACT‑CCCCTCGCTCGCCCCGCGGGCAGCCCCCTCCCCGCGTTCAc                                                                                            >  2:609736/1‑151 (MQ=17)                                                      cccTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAc                                                                                            <  2:296741/150‑1 (MQ=255)                                                           aCAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTg                                                                                      <  2:616162/151‑1 (MQ=255)                                                                cGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCaa                                                                                 <  1:37421/151‑1 (MQ=255)                                                                     gCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AGGGGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑Tgg                                                                            <  2:486502/151‑1 (MQ=255)                                                                      cGGGGGCACCCGGTGTGTGCAATATACG‑‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTg                                                                         <  1:559331/151‑1 (MQ=255)                                                                       gggggCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTg                                                                         <  2:180427/151‑1 (MQ=21)                                                                                       tgCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGGTTGACT‑CGCTTCTCTCCCCCTTCGGGCAGCCCGCCCCCCGTGTTCACGGTCTGTCCCAC‑CGGCTGTCGCCC‑Ggtg‑ttggacc                                                      >  1:385474/1‑145 (MQ=1)                                                                                          aaTATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑Acccc                                                    <  1:868774/150‑1 (MQ=14)                                                                                              taCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAA‑‑GGGGGGGATTGAAT‑CGCCTTCCCCCCCCCGCGGGGAGCCCGCCCCCTGCGTTTACGGGCTTTCCAAC‑CTGCTGTCGCCC‑CTTT‑TCG‑Acaccc                                                   >  2:752579/1‑147 (MQ=2)                                                                                                 g‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCggg                                            <  1:503017/151‑1 (MQ=14)                                                                                                   aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCgggg                                           <  1:852492/151‑1 (MQ=14)                                                                                                    aaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGc                                          <  1:781183/151‑1 (MQ=14)                                                                                                        aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTtct                                      <  1:1026320/151‑1 (MQ=14)                                                                                                         aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGGTTTACC‑CGCTTCCCTCGCCCTCCGGGCCGCCCCCTCACTGCGTTTCCGGTCTGTCCCAC‑CGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCCGTCGGGGGtt                                        >  2:195174/1‑148 (MQ=11)                                                                                                          aGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGGTTTACT‑CGCCTCTCTCCCCCCCCGGGCAGCCCGCTCACTGCGGTCACGGTCTCTCCCAC‑CCGCCGTCGCCC‑GTTG‑Tcccaccccccgcggg                                             >  1:316744/1‑128 (MQ=0)                                                                                                                tACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑G‑A‑‑GGGGGGGATTTAGT‑CTGCTCTCCCGCTCTTCCGGTCGCCCGCTCCCTTCGGTTACGGTCTGTTCCAC‑CGGCGGGCCCCC‑CCTT‑TTG‑ACCCCCCGTCGGGGCTTCTCCTcccccc                             >  1:69891/1‑151 (MQ=2)                                                                                                                  cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGG‑G‑‑GGGGGGGGTTGACTCCCCTCCGCCCCCCCTG‑GGGGAGCCCCCCCCCCCCGTTTTCCGGCTGTCCAACATGGC‑GTCCCCACGTGT‑TCG‑ACCCCCGGCGGGGGCTT‑TCATccccccc                            >  2:1105735/1‑150 (MQ=2)                                                                                                                      cGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGgtg                         <  2:175014/150‑1 (MQ=14)                                                                                                                      cGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGAGTGACT‑CGCTTCCCTCGCCCCGCGGGGCCCCCGCCCACCGCGTTCTCCGGCTGTTCAAC‑TGGCTGTCCCCC‑GGTG‑TCG‑ACCCCCCGCCGGGGCTTTCCCTccccccc                            >  1:213752/1‑146 (MQ=2)                                                                                                                       gTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGGTTTACT‑CCCCCCCCTCGCCCTCCTGGCAACCCGCTCCCTGGGTTTACGGGCCGTCCCAC‑TGGGCGTTGGCA‑GCTT‑TCGCACCCCCGG‑CGGGGGCTCTCCTcccccccggg                         >  1:147692/1‑144 (MQ=0)                                                                                                                            aGCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCaa                  <  1:244920/151‑1 (MQ=14)                                                                                                                            aGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGGTTGACT‑CGCTTCGCTCCCCCCGCGGGGCGCCCGCCCCCCGCGGTTACGGGCTGTTCCAC‑CGGCTGGCGCCC‑GGTGTTGG‑ACCCCCGG‑CGGGGGGTCTCCTcccccccggg                         >  1:1177162/1‑139 (MQ=0)                                                                                                                             gCTCTTCTTTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAAt                 <  2:16445/151‑1 (MQ=14)                                                                                                                                  tcttTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATacac             <  1:118941/150‑1 (MQ=18)                                                                                                                                   cttTAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATac               <  2:126532/147‑1 (MQ=17)                                                                                                                                    tttAAATATGGC‑GGT‑G‑AG‑GGGGGGGGTTGACT‑CGCTTCGCCCGCCCTGCGGGCAGCCCGCCCACTGCGTTCACGGGCTGTCCAAC‑TGGCTGTCGCCC‑GTTG‑TCG‑ACCCCCCGGCGGGGCTTCTCCTCCCCCCCGGGGGTGGAAAaca              >  2:987002/1‑147 (MQ=11)                                                                                                                                     ttAAATATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACg            <  2:160346/148‑1 (MQ=18)                                                                                                                                       aaaTATGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGaaaaa       <  2:832558/151‑1 (MQ=18)                                                                                                                                           atGGC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGaaaaaaaaa   <  1:479084/151‑2 (MQ=18)                                                                                                                                             ggC‑GGT‑G‑AG‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGc  <  2:937232/150‑1 (MQ=18)                                                                                                                                                        |                                                                                                                                                    ATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑G‑A‑‑GGGGGGGATTGACT‑CGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC‑TGGCTGTCGCCA‑GTTG‑TCG‑ACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGC  >  NC_000913/1096884‑1097171

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