BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,769,478	(G)_7→8	coding (405/1113 nt)	ydiK →	UPF0118 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,769,471	1	.	G	89.7%	250.2 / 18.7	78	coding (398/1113 nt)	ydiK	UPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (3/5); new base G (36/34); total (39/39)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769337‑1769596

tATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑g                                                                                                                             >  2:729260/1‑139 (MQ=255)

  tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑g                                                                                                                          >  1:1973466/1‑139 (MQ=255)

  tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑g                                                                                                                          >  2:1745965/1‑139 (MQ=255)

        aCGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑gggacgg                                                                                                                    <  2:1973466/139‑6 (MQ=255)

          ggCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCg                                                                                                                  <  1:1745965/139‑1 (MQ=255)

            cAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGAt                                                                                                                >  2:32783/1‑139 (MQ=255)

             aGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGTGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGG‑GG‑GG‑ACGGCGAt                                                                                                                <  1:359974/139‑1 (MQ=255)

              gCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCa                                                                                                              <  1:729260/139‑1 (MQ=255)

                        aaTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGt                                                                                                    <  1:615899/139‑1 (MQ=255)

                         aTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATTGCGAAAGTc                                                                                                   <  2:1379900/139‑1 (MQ=255)

                                   tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 >  2:1512343/1‑131 (MQ=255)

                                   tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 <  1:1512343/131‑1 (MQ=255)

                                   tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTtata                                                                                         >  1:471547/1‑139 (MQ=255)

                                       cAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTgg                                                                                     <  2:471547/139‑1 (MQ=255)

                                          cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcac                                                                                  <  1:257120/139‑1 (MQ=255)

                                          cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCATTATAtt                                                                                       >  2:54345/1‑134 (MQ=255)

                                          cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GGTACGGCGATCATGGCGAAAGTCCGCCCTTATAtt                                                                                       <  1:54345/134‑1 (MQ=255)

                                             tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCg                                                                                                                  <  1:955012/104‑1 (MQ=255)

                                             tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCg                                                                                                                  >  2:955012/1‑104 (MQ=255)

                                             tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   >  2:1292669/1‑119 (MQ=255)

                                             tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   <  1:1292669/119‑1 (MQ=255)

                                               aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                             >  1:679656/1‑139 (MQ=255)

                                                cATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccac                                                                            <  2:679656/139‑1 (MQ=255)

                                                    aCGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCaccaccaccacc                                                                        >  2:1201357/1‑139 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ATGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                           >  1:437360/1‑131 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 >  1:905161/1‑109 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 <  2:905161/109‑1 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc                                                                           <  2:437360/131‑1 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacc                                                                              <  2:1674202/128‑1 (MQ=255)

                                                         aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacc                                                                              >  1:1674202/1‑128 (MQ=255)

                                                               tCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                          <  2:754653/110‑1 (MQ=255)

                                                               tCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                          >  1:754653/1‑110 (MQ=255)

                                                               tCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTgg                                                             >  1:770300/1‑139 (MQ=255)

                                                                cTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑g                                                                                                                          >  1:496648/1‑77 (MQ=255)

                                                                cTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑g                                                                                                                          <  2:496648/77‑1 (MQ=255)

                                                                   gCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGTGG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAg                                                         <  1:1201357/139‑1 (MQ=255)

                                                                   gCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        >  1:1471254/1‑139 (MQ=255)

                                                                    cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   >  2:1915775/1‑96 (MQ=255)

                                                                    cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   <  1:1915775/96‑1 (MQ=255)

                                                                    cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        >  2:689134/1‑139 (MQ=255)

                                                                        gCTTAATACCATTCCGGTGTTTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCCTGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGATCGTTGGGCAggcggc                                                    >  2:88039/1‑139 (MQ=255)

                                                                        gCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAggcggc                                                    >  2:10209/1‑139 (MQ=255)

                                                                         cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        <  2:1036174/134‑1 (MQ=255)

                                                                         cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        >  1:1036174/1‑134 (MQ=255)

                                                                         cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGgcggcg                                                   >  2:485354/1‑139 (MQ=255)

                                                                         cTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGgcgc                                                  >  2:640672/1‑139 (MQ=255)

                                                                           tAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gg                                                                                                                                <  1:565757/62‑1 (MQ=255)

                                                                           tAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gg                                                                                                                                >  2:565757/1‑62 (MQ=255)

                                                                                  aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                          >  2:2002520/1‑139 (MQ=255)

                                                                                   ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGt                                                                                                    >  1:606378/1‑80 (MQ=255)

                                                                                   ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGt                                                                                                    <  2:606378/80‑1 (MQ=255)

                                                                                   ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgc                                         >  2:1602348/1‑139 (MQ=255)

                                                                                     ccGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      <  2:1471254/139‑1 (MQ=255)

                                                                                      cGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      <  1:32783/139‑1 (MQ=255)

                                                                                        gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                   <  1:1708882/92‑1 (MQ=255)

                                                                                        gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                   >  2:1708882/1‑92 (MQ=255)

                                                                                        gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑gggacgg                                                                                                                    <  2:259917/59‑6 (MQ=255)

                                                                                        gTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg‑gg‑gggacgg                                                                                                                    >  1:259917/1‑54 (MQ=255)

                                                                                             tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 <  2:428458/73‑1 (MQ=255)

                                                                                             tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 >  1:428458/1‑73 (MQ=255)

                                                                                                cgcgAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCa                            <  1:1442360/139‑1 (MQ=255)

                                                                                                     aGCTGTATGCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACGACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATtgtg                       >  2:1096013/1‑139 (MQ=255)

                                                                                                             gCAGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATg               <  1:689134/139‑1 (MQ=255)

                                                                                                               aGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGaaa                                                                                                      >  1:487279/1‑50 (MQ=255)

                                                                                                               aGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGaaa                                                                                                      <  2:487279/50‑1 (MQ=255)

                                                                                                               aGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        >  1:1317583/1‑96 (MQ=255)

                                                                                                               aGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        <  2:1317583/96‑1 (MQ=255)

                                                                                                               aGGCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCt             >  2:1672332/1‑139 (MQ=255)

                                                                                                                ggCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCtt            <  1:2002520/139‑1 (MQ=255)

                                                                                                                 gCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 >  1:915821/1‑133 (MQ=255)

                                                                                                                 gCTGGCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 <  2:915821/133‑1 (MQ=255)

                                                                                                                    ggCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTcc                                                                                                  >  2:987852/1‑49 (MQ=255)

                                                                                                                    ggCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTcc                                                                                                  <  1:987852/49‑1 (MQ=255)

                                                                                                                    ggCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 <  2:1299015/50‑1 (MQ=255)

                                                                                                                    ggCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCg                                                                                                 >  1:1299015/1‑50 (MQ=255)

                                                                                                                    ggCACAACTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct        >  1:1631790/1‑139 (MQ=255)

                                                                                                                          tcTTGTTGGATAT‑‑GGG‑GGTGG‑AGGGCGGTCATGGCGTAAGTCCGCCCTTATATGGGCACCACCACCACCTGATTCGGTGGGCAGGCCGCGCATATCGCGCGCGTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGta  <  1:88039/138‑2 (MQ=255)

                                                                                                                          aCTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg  <  1:10209/139‑1 (MQ=255)

                                                                                                                           cTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   >  2:606448/1‑41 (MQ=255)

                                                                                                                           cTTGCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTc                                                                                                   <  1:606448/41‑1 (MQ=255)

                                                                                                                           cTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGc                                                                                                >  1:429490/1‑43 (MQ=255)

                                                                                                                           cTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGc                                                                                                <  2:429490/43‑1 (MQ=255)

                                                                                                                              gCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGccc                                                                                              >  1:257926/1‑43 (MQ=255)

                                                                                                                              gCTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGccc                                                                                              <  2:257926/43‑1 (MQ=255)

                                                                                                                               cTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc                                                           <  2:181555/77‑1 (MQ=255)

                                                                                                                               cTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc                                                           >  1:181555/1‑77 (MQ=255)

                                                                                                                               cTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAggcgg                                                     <  2:1142176/83‑1 (MQ=255)

                                                                                                                               cTGGATATG‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAggcgg                                                     >  1:1142176/1‑83 (MQ=255)

TATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGG‑GG‑GG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769337‑1769596

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