BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,898	Δ1 bp	coding (212/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,897	0	T	.	98.2%	515.1 / 0.2	114	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base . (47/64); minor base C (0/2); total (48/66)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACA  >  NC_000913/1793749‑1794032

 cgcgACTGCGGTCGGCGCGGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTGTAATCTCTTCATCGCTGCCCGGTTCtg                                                                                                                                        >  2:1052559/1‑148 (MQ=255)

   cgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCCTCATCGCTGCCCGGTTCtgtt                                                                                                                                      <  1:1100345/149‑4 (MQ=255)

   cgACTGCGGGCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCCATCTTCCTGCACATCGCGTACGTCACCCCCCATGTGGTGGTACAGTAAGTCTTAATCTCTTCATCGCTGCCCGGTTCtgtt                                                                                                                                      <  2:53071/149‑4 (MQ=255)

    gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                     <  2:1393819/149‑5 (MQ=255)

    gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                     >  1:220958/1‑145 (MQ=255)

       tGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                 >  2:1394597/1‑149 (MQ=255)

       tGCGGTCGGCGCTGCGGCAATCAATTTTGGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                 <  1:789529/149‑1 (MQ=255)

       ggCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGATCGTCTCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                 <  2:795899/148‑1 (MQ=255)

          ggTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACTTCGGGAACGTCACCCCCCATGTGGTGGGACAGTAAGGTTTAATCTCTTCATAGCTGCCCGGTTC‑TGTTCCAGcaa                                                                                                                              >  2:1558912/1‑149 (MQ=255)

              ggCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTgg                                                                                                                          <  2:1432841/149‑1 (MQ=255)

                cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGtt                                                                                                                        <  2:1281195/149‑1 (MQ=255)

                      ggCAATGTATTGTCGATACAGCGGATGTGCCCCTTCGCCAGTAACTTCAATCTTACTGAACACCGGAGACGTCATCCCTCAGGTGGTGGTACAGGAAGTGTTAATCCCTCCGTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                  <  2:500224/149‑1 (MQ=255)

                         aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGaa                                                                                                               >  2:1768763/1‑149 (MQ=255)

                            caatTTTTGATACAGCGGATGGCGTCCTGCGCCATAAACTTCAATCTTACTGAACATCGGGAAAGTCACCCACCGTAAGGCGTTACCGTAACTTTTAACCTCTTCATCGCTGACCGGGT‑ATGTTCCCGAAACTCGTTGCACGGGAATcc                                                                                                            >  2:304393/1‑149 (MQ=255)

                              atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCa                                                                                                          >  1:1522444/1‑149 (MQ=255)

                                 tttGAGACAGCGGATGGCGTCCTTCGCCATTAAGTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTCTAATCTCTTCATCGATACCCGCTTC‑TGTTCAAGAAACCGGTTGCACTGTAATCCCAGCa                                                                                                       <  2:1491071/149‑1 (MQ=255)

                                 attGAGACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCATACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGGTCCAGAAACTGGTTGCACGGGAATCCCAGCa                                                                                                       <  1:1671905/148‑1 (MQ=255)

                                     aTACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTTGTTGCACGGGAATCCCAGCACCAt                                                                                                   <  1:1325500/149‑1 (MQ=255)

                                         aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                               >  2:1480895/1‑149 (MQ=255)

                                         aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCCCCCCCGCTTCCCTTTTC‑TGTTCCAGAAACTGGTTGCCCGGGAATCCCAGCACCATaaaa                                                                                               >  1:51370/1‑149 (MQ=255)

                                           cGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                                                                             >  1:1454430/1‑149 (MQ=255)

                                             gATGGCCTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCtc                                                                                           <  2:312812/149‑1 (MQ=255)

                                               tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                                                         <  1:467585/149‑1 (MQ=255)

                                                tgCGTCATTCGTAATTAACTTCAATCTTACTGAACATTGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                        <  2:535244/148‑1 (MQ=255)

                                                ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                        <  1:1729500/149‑1 (MQ=255)

                                                      cTTCGCCGTTCACCTCAATCCTACGGAACATGGTGACGTTCACCCCCCCTGTGTTGGTACAGTAAGTTTTACTCTCTTCTTCGCTTCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCGCGATCTAccc                                                                                  <  2:335753/149‑1 (MQ=255)

                                                      cTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAccc                                                                                  <  1:1536189/149‑1 (MQ=255)

                                                       ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGATC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                 >  2:881907/1‑149 (MQ=255)

                                                       ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGGACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCa                                                                                 >  1:792697/1‑149 (MQ=255)

                                                         cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAgg                                                                               <  1:395338/149‑1 (MQ=255)

                                                           ccATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                                                                             >  1:701883/1‑149 (MQ=255)

                                                           ccATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                                                                             >  2:561062/1‑149 (MQ=255)

                                                            taTTAAGTTCAATCTTACCGAACATCGGGAACGTCACCACCCCTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCtt                                                                            <  2:281114/148‑1 (MQ=255)

                                                             aTTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATCAAACCTCGATCGACCCAGGCttt                                                                           <  1:943576/149‑1 (MQ=255)

                                                             aTTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCttt                                                                           <  1:1448574/149‑1 (MQ=255)

                                                               tAACTTCAATCTTACTGGACATCGGGAACGTCACCCCCCATGTGGTTGTACAGGAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACGGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                         <  2:238840/149‑1 (MQ=255)

                                                               tAACTTCAATCTTACTGCACATCGGGAACGTCACCCCCCAGGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                         <  1:1144990/149‑1 (MQ=255)

                                                                 aCTTCAATCTTACTGAACATCGGGAACGTCGCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGa                                                                       <  2:572849/149‑1 (MQ=255)

                                                                   ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                     >  1:970937/1‑149 (MQ=255)

                                                                      aaTCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATAtt                                                                  <  1:272532/149‑1 (MQ=255)

                                                                         cTTACGGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                               <  1:1697633/149‑1 (MQ=255)

                                                                          ttACTGAACATCGGGAACGTGCCCCCCCGTGGGGTGGTACAGTAAGTTTTAATCTCTTCTTCGCTGGCCGGTTC‑TGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                                              <  2:1660515/149‑1 (MQ=255)

                                                                          ttACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                                              <  2:1731906/149‑1 (MQ=255)

                                                                              tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                          >  2:160574/1‑149 (MQ=255)

                                                                              tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                          >  1:1013646/1‑149 (MQ=255)

                                                                              tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGAGCGACCCAGGCTTTCTGAATATTCTCCAACt                                                          >  1:165791/1‑149 (MQ=255)

                                                                                aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGACACTGGCTGCACGGGAACCCCCGCACCATAAACCCTAGATCGCCACAGGCTCGCTGAACATGCTCCATCTGc                                                        <  2:1313027/149‑1 (MQ=255)

                                                                                aaCATCGGGAACGTCACCCCCCATGTGGCGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCCCGGGACTCCCAGCCCCATACAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACtac                                                        >  2:679744/1‑147 (MQ=255)

                                                                                 aCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                                       <  1:721673/149‑1 (MQ=255)

                                                                                   aTCGGGACCGTAGCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                                     <  1:1557757/149‑1 (MQ=255)

                                                                                       ggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGACACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGCTCGACCCAGGCTTCCTGCATATTCTCCAACTGCCCATAtt                                                 >  1:616849/1‑149 (MQ=255)

                                                                                         aaCGTCGCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                               <  1:1028638/149‑1 (MQ=255)

                                                                                           cgccACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcgg                                             <  2:82067/146‑1 (MQ=255)

                                                                                            gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                            <  2:1293036/149‑1 (MQ=255)

                                                                                             tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCACCTGCTCATATTgcggcg                                           >  1:1421529/1‑149 (MQ=255)

                                                                                               aCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGtt                                         >  1:14211/1‑149 (MQ=255)

                                                                                                 cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                       <  1:83810/149‑1 (MQ=255)

                                                                                                 cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTGGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCAAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGCCGTTaa                                       >  2:231646/1‑149 (MQ=255)

                                                                                                  ccccATTTCGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAATACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                      <  2:1611158/149‑1 (MQ=255)

                                                                                                  ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                      >  1:378536/1‑149 (MQ=255)

                                                                                                     cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                                   <  2:796310/149‑1 (MQ=255)

                                                                                                        gtggtgGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCacac                                <  2:1479329/149‑1 (MQ=255)

                                                                                                          ggtgggACAGTAAGCTTTAATCTCTTCATCGCTCCCCGGTTC‑TGGTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCGGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACtt                              <  1:578195/149‑1 (MQ=255)

                                                                                                           gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                             >  1:1559695/1‑149 (MQ=255)

                                                                                                           gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                             <  2:690479/149‑1 (MQ=255)

                                                                                                            tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTg                            <  1:623244/149‑1 (MQ=255)

                                                                                                             ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                           <  1:1152953/149‑1 (MQ=255)

                                                                                                               tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                         <  1:497353/149‑1 (MQ=255)

                                                                                                                aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTCCGGCGTTAAGCCACACTTTGAGGc                        >  2:1563665/1‑149 (MQ=255)

                                                                                                                aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTACACGGGAATACCAGCACCATAAAACCTCGATCGACCCACGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCGCACTTTGAGGc                        >  1:920178/1‑149 (MQ=255)

                                                                                                                 cAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACCCTTTGAGGCg                       >  2:1385164/1‑149 (MQ=255)

                                                                                                                  aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                      >  2:151285/1‑149 (MQ=255)

                                                                                                                  aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                      <  2:721210/149‑1 (MQ=255)

                                                                                                                  aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                      >  2:217128/1‑149 (MQ=255)

                                                                                                                  aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGAACGACCCAGGCTTTCTGAATATTATCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                      >  2:784126/1‑149 (MQ=255)

                                                                                                                  aGTAAGTTTTAATCTCTGCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                      <  2:694552/149‑1 (MQ=255)

                                                                                                                   gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCCGGCTTTCTGAATATTCTCCAACTGCTCATCTTGCGGCGTTATGCCAACTTTTGAGGCCAc                     >  1:531329/1‑149 (MQ=255)

                                                                                                                   gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAc                     <  2:1218462/149‑1 (MQ=255)

                                                                                                                   gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAc                     <  2:1121332/149‑1 (MQ=255)

                                                                                                                      aGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTTAGGCGACAtt                  >  1:740242/1‑149 (MQ=255)

                                                                                                                        ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                <  2:1454430/149‑1 (MQ=255)

                                                                                                                        ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                >  2:1354453/1‑149 (MQ=255)

                                                                                                                        ttttAATCTCTTCATCGCGGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                <  2:1031722/149‑1 (MQ=255)

                                                                                                                         tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc               >  1:1433535/1‑149 (MQ=255)

                                                                                                                         tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc               <  1:1362951/149‑1 (MQ=255)

                                                                                                                          ttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa              <  2:701883/149‑1 (MQ=255)

                                                                                                                              tctcTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCa          >  1:665522/1‑149 (MQ=255)

                                                                                                                              tctcTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCa          >  1:577097/1‑149 (MQ=255)

                                                                                                                              tctcTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGACGACATTGACACTcc          >  2:1326333/1‑148 (MQ=255)

                                                                                                                               ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa         >  2:14201/1‑149 (MQ=255)

                                                                                                                               ctctTCATCGCTGCCCGG‑TCCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa         <  1:296619/149‑1 (MQ=255)

                                                                                                                                tctTCATCGCTGCCCGGTTC‑TGTTTCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc        >  1:1328053/1‑149 (MQ=255)

                                                                                                                                tctTCATCGCTGCCCGGTTC‑TGTTCCAGACGTTGGTGGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc        >  1:874084/1‑149 (MQ=255)

                                                                                                                                tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc        <  1:253074/149‑1 (MQ=255)

                                                                                                                                 ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa       <  1:4955/149‑1 (MQ=255)

                                                                                                                                  ttCATCGCTGCCCGGTTC‑TCTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg      <  1:1511711/149‑1 (MQ=255)

                                                                                                                                  ctCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg      <  1:1257638/148‑1 (MQ=255)

                                                                                                                                    cATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATGGACAATCAACAGca    <  2:345487/149‑1 (MQ=255)

                                                                                                                                     aTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                  >  1:127715/1‑38 (MQ=255)

                                                                                                                                     aTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCAACGCACCATAAAACCTCGGTCGACCCAGACTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGcac   >  1:866976/1‑149 (MQ=255)

GCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACA  >  NC_000913/1793749‑1794032

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