Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,031	0	C	T	39.7%	80.5 / 36.2	63	intergenic (‑94/+109)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base C (29/9);  new base T (0/25);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 2.31e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTT  >  NC_000913/1824881‑1825170                                                                                                                                                         |                                                                                                                                             aCCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGc                                                                                                                                              >  2:867886/1‑151 (MQ=255)    ggCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCttt                                                                                                                                           >  1:7398/1‑151 (MQ=255)     gcagcaTTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCtttt                                                                                                                                          <  1:159657/151‑1 (MQ=255)                 cacaGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGtt                                                                                                                              >  1:1145847/1‑151 (MQ=255)                  acaGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGttt                                                                                                                             >  1:897956/1‑151 (MQ=255)                  acaGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGtt                                                                                                                              >  1:458728/1‑150 (MQ=255)                   caGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTg                                                                                                                            <  2:1014212/151‑1 (MQ=255)                     gATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTgag                                                                                                                          <  1:752307/151‑1 (MQ=255)                     gATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTga                                                                                                                           <  2:1218496/150‑1 (MQ=255)                        tCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTgagaa                                                                                                                        <  1:45339/150‑1 (MQ=255)                              aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                  >  2:1143455/1‑150 (MQ=255)                                   ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                             <  2:757885/150‑1 (MQ=255)                                   ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                             <  2:279572/150‑1 (MQ=255)                                     aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑cc                                                                                                         >  2:315776/1‑151 (MQ=255)                                       gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ccgc                                                                                                       <  1:347518/151‑1 (MQ=255)                                        cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ccgcc                                                                                                      >  1:212596/1‑151 (MQ=255)                                           atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGAt                                                                                                   >  1:562799/1‑151 (MQ=255)                                             cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATcc                                                                                                 >  1:797539/1‑151 (MQ=255)                                              aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCt                                                                                                >  1:503083/1‑151 (MQ=255)                                                  tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTcc                                                                                            >  1:12508/1‑151 (MQ=255)                                                   aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTcc                                                                                            >  2:493941/1‑150 (MQ=255)                                                    ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCCACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ACGCCGATCCTGTTCAc                                                                                          >  2:679303/1‑151 (MQ=255)                                                     tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCAcc                                                                                         >  2:254037/1‑151 (MQ=255)                                                     tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCAcc                                                                                         >  1:16385/1‑151 (MQ=255)                                                     cccATCCTCCCTCCCGGCCTTTTTTTTTCATTGTTATCTTTC‑GGTGACGGGGGGATGGTTATCGGCAATGCGGGTTATTTCTGTTGG‑TGGGGTTCAGGTTTTTTTTACAACTTTTTAAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCAcc                                                                                         <  2:557961/150‑1 (MQ=255)                                                       cATACTCGCCCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTACC‑CGGGGATGTTTATCGGCACCTCGGGTTTTTTACGTTAG‑TGGTGATCAGGTTTTTTTTACACCTGTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGc                                                                                       <  1:25305/151‑1 (MQ=255)                                                               cTCCCGGCCCGTTTTTTTCTTTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGCCAACTGGGGTGATTTACGTTAG‑TGGTGATCGGTTTTTTTTTACAACTTTTTGAAAAGAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACt                                                                                 <  1:130870/149‑1 (MQ=255)                                                               cTCCCGCCCGTTTCTTTTCATTGTTCTCTTACAGGAAAC‑CGGGGATTTTTCTCTGCAACCGGGGAAATTTACGGTAG‑GTGCAATCAGGTTTTTTTTCCACTTTTGGGAAAAGAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACt                                                                                 <  2:907875/149‑1 (MQ=255)                                                                 cccGGCTCTTTTTTTTGTTTTTTATCTTTAAGGAAAC‑CGTGGATTGTTTGCGGCAGCTGGGGTAATTTCCGTTGG‑TGGTCAGCGTGTTTTTTTAACAACTGTTGGAGAAAGAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCAt                                                                                     <  2:412126/143‑1 (MQ=255)                                                                  ccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCa                                                                              >  2:545075/1‑149 (MQ=255)                                                                   cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                           >  1:795819/1‑151 (MQ=255)                                                                   cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                           >  2:355071/1‑151 (MQ=255)                                                                    gcccggTTTTTTTTTTTTTTATAGTTCAGGAAAC‑GGGGGTTGTTTAGCGTAAACTGGGGTATTTTACGTTAG‑TGGTTATCAGGTTTTTTTTACAGCTTTTAGAGAAGAGAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                           <  2:702771/145‑1 (MQ=255)                                                                         gTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                     >  1:1004887/1‑151 (MQ=255)                                                                             ttttCTTTTTTTTCGTTCAGTTACC‑CGTTGATTGTTACCGCCAACTGGGTTGTTTTACTTTAG‑TGGTACTCAGGTTTTTTTTAAATTTTTATAAAAAAAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAATGTAGTGGAc                                                                   <  2:697794/149‑1 (MQ=255)                                                                                  aTTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                            >  2:162334/1‑151 (MQ=255)                                                                                   tttttATCGTCCAGGGAAC‑GGTGTATGTTTATCGGCAACTCGGGTGATTTACGTTGG‑TGGTGATCAGGTTTTTTTCAAAATTGTTTGAAAAAAGAAAAGAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                            <  2:1076058/147‑1 (MQ=255)                                                                                     gTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                         >  2:1167952/1‑151 (MQ=255)                                                                                      ttatcgtcggggcagc‑cGTAGTTGGTTATCGTCAACTTGGTTTATTTACTTGAG‑TGGTGTTAATGTTTTTTTACCATCTATTTGAAAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                         <  1:736249/135‑1 (MQ=255)                                                                                       tATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                       >  1:467748/1‑151 (MQ=255)                                                                                        aTCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                      >  1:57405/1‑151 (MQ=255)                                                                                             tCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                    >  2:1116331/1‑148 (MQ=255)                                                                                                gttacc‑cGTGGATGTTTATCGGCACCTCGGTTGATTTATTTTAG‑TGGTATTCAGTTTTTTTTTACAATTGTTTGAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTaa                                                  <  2:412965/142‑1 (MQ=255)                                                                                                ggTACC‑CGTGGATGTTTATCGGCAACGGGGGTATTTTACTTTAG‑GGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                              <  1:944344/151‑1 (MQ=255)                                                                                                     c‑cGTGGATGTTTATCGGCATCTCGGGTGTTTTACTTTAG‑TGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                         <  1:76043/151‑1 (MQ=255)                                                                                                      cggTGGATGTTTATCGGAAACGGGGGTGATTTACTTTAG‑TGGTGATCGGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                          <  1:1271287/148‑1 (MQ=255)                                                                                                          ggATGGTTTCCGGCATCTGGGGTTATTTTGTTTAG‑TGGTTATCAGTTTTTTTTTCCAACTTTTAGAGAAGAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                      <  2:731369/150‑1 (MQ=255)                                                                                                          ggATGGTTCCCGGAAGCTGGGGTATTTTGCTTGAG‑TGGTGATCGGCTTTTTTTTACATATGTTTGAGAAAAAAAAAAAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                            <  1:517150/144‑1 (MQ=255)                                                                                                                ttATCGGCAACTCGGGTGTTTTACGTTAG‑TGGTGATCAGTTTTTTTTTACATCTGTTGGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa                                 <  1:185490/149‑1 (MQ=255)                                                                                                                  atcgccaccggggggGATTTA‑GTTAGTTGGTGACCAGTTTTTTTTTACAACTGTTGGGAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                              <  2:513914/139‑1 (MQ=255)                                                                                                                     ggCACCTGGGGTGATTTACGTTAG‑TGGTGATCAGGTTTTTTTTACAACTGTTTAAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                           <  1:466218/150‑1 (MQ=255)                                                                                                                     gacaatttgggTGATTTACGTATG‑TGGTGATCGGGTTTTTTTTACATATTTTAGAGAAAAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat                             <  2:706031/140‑1 (MQ=255)                                                                                                                     gacAACTGGGGTGATGTAGGTTAG‑GGGTAGTCGTGTTTTTTTAACAGTTTTTTGGAAAGA‑AAAAGAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                              <  2:845721/145‑1 (MQ=255)                                                                                                                       cAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                        >  1:1092414/1‑151 (MQ=255)                                                                                                                        acctggggTGTTTTACTTTAG‑TGGTGAGCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                       <  1:975006/146‑1 (MQ=255)                                                                                                                               gTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                  >  1:824079/1‑149 (MQ=255)                                                                                                                                tGATTTACGTTAG‑GGGTGATCAGGCTTTTTTTACAACTGGTTGAGAAGAGAAAAGAACA‑CCGCCGATCCTGTCAACCGCATTACTGCAAGGTAGTGGACACGACCGGCGGTCTTTAGATTTTTGGCGGAAAGGTTATTCAGCCGTTGCAgg               >  2:602144/1‑151 (MQ=255)                                                                                                                                   tttagtttAG‑GGGTGATCAGTTTTTTTTTACATTTGTTGAAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                 <  1:233587/140‑1 (MQ=255)                                                                                                                                   tttACTTTAG‑TGGTGATCAGGTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                <  1:1203296/147‑1 (MQ=255)                                                                                                                                   tttACGTTAG‑TGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg               <  1:636270/148‑1 (MQ=255)                                                                                                                                        gTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAGCAGTTGCAGGCAtttt         >  2:1031293/1‑149 (MQ=255)                                                                                                                                         ttAG‑TGGTGATCGGGTTTTTTTAAAAACTGTTTGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt           <  1:723852/146‑1 (MQ=255)                                                                                                                                          tAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt     >  1:359856/1‑151 (MQ=255)                                                                                                                                               ggtgATCGTGCTTTTTGTCCATTTGTTAGAAAAGAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt  <  2:7397/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                             ACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTT  >  NC_000913/1824881‑1825170

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