BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,841,922	T→C	Y325C (TAC→TGC)	b1762 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,841,922	0	T	C	99.5%	722.7 / ‑4.6	205	Y325C (TAC→TGC)	b1762	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (103/101); total (103/102)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 206 total aligned reads displayed.

GCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAA  >  NC_000913/1841832‑1842011

 cccccGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGt                                                                                        <  1:3248447/93‑1 (MQ=255)

      ggCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTcc                                                                                   <  1:4694580/93‑1 (MQ=255)

       gCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTccc                                                                                  <  1:3219746/93‑1 (MQ=255)

       gCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTccc                                                                                  <  2:6246309/93‑1 (MQ=255)

       gCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTccc                                                                                  <  1:2158834/93‑1 (MQ=255)

        cAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCa                                                                                 <  2:456888/93‑1 (MQ=255)

         aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCaa                                                                                >  1:1386198/1‑93 (MQ=255)

         aaaaTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCaa                                                                                >  1:1825844/1‑93 (MQ=255)

            aTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAacac                                                                             >  2:2387439/1‑93 (MQ=255)

            aTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAAc                                                                             <  2:1527909/93‑1 (MQ=255)

            aTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAAc                                                                             <  2:4937958/93‑1 (MQ=255)

             tttCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACt                                                                            <  2:5628518/93‑1 (MQ=255)

               tCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGc                                                                          <  2:476380/93‑1 (MQ=255)

                cAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCt                                                                         <  1:1168709/93‑1 (MQ=255)

                 aGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCtt                                                                        <  1:2167121/93‑1 (MQ=255)

                 aGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCtt                                                                        >  2:654231/1‑93 (MQ=255)

                  ggCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTg                                                                       >  2:5703213/1‑93 (MQ=255)

                  ggCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTg                                                                       >  2:4388929/1‑93 (MQ=255)

                     gTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCg                                                                    >  1:2522411/1‑93 (MQ=255)

                      tttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGt                                                                   >  2:2910262/1‑93 (MQ=255)

                      tttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGt                                                                   <  2:5892018/93‑1 (MQ=255)

                       ttATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTGTGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTa                                                                  >  1:4523647/1‑93 (MQ=255)

                        tATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTaa                                                                 >  2:6367022/1‑93 (MQ=255)

                         aTGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAAc                                                                >  1:5751390/1‑93 (MQ=255)

                          tGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACg                                                               >  1:557671/1‑93 (MQ=255)

                           gAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGa                                                              <  2:6430527/93‑1 (MQ=255)

                           gAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGa                                                              <  1:4018005/93‑1 (MQ=255)

                             gTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATa                                                            >  2:32125/1‑93 (MQ=255)

                              tATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATaa                                                           >  2:6172582/1‑93 (MQ=255)

                              tATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATaa                                                           >  1:1961032/1‑93 (MQ=255)

                               aTTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAac                                                          <  2:3020147/93‑1 (MQ=255)

                                  tAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACg                                                       >  2:777161/1‑93 (MQ=255)

                                    aCGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCt                                                     <  2:3326129/93‑1 (MQ=255)

                                      ggATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTc                                                   >  2:6087925/1‑93 (MQ=255)

                                       gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                  >  2:4265330/1‑93 (MQ=255)

                                       gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                  >  2:369401/1‑93 (MQ=255)

                                       gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                  <  1:1650820/93‑1 (MQ=255)

                                       gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                  <  1:67636/93‑1 (MQ=255)

                                       gatgatGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCa                                                  >  2:2602753/1‑93 (MQ=255)

                                          gatgCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGt                                               >  2:3725758/1‑93 (MQ=255)

                                              cTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCaa                                           >  2:1639971/1‑92 (MQ=255)

                                              cTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAg                                           >  1:5740410/1‑93 (MQ=255)

                                               tCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGc                                          <  1:4537378/93‑1 (MQ=255)

                                               tCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGc                                          >  2:156253/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  2:5434511/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  1:1692478/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  2:235335/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  2:1001183/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  1:5900242/1‑93 (MQ=255)

                                                ccccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCa                                         >  1:4911197/1‑93 (MQ=255)

                                                 cccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                        <  1:3004893/93‑1 (MQ=255)

                                                 cccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                        >  2:1664177/1‑93 (MQ=255)

                                                 cccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCAt                                        >  1:4230778/1‑93 (MQ=255)

                                                  ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATc                                       <  2:1416204/93‑1 (MQ=255)

                                                  ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATc                                       >  1:3546229/1‑93 (MQ=255)

                                                    aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCg                                     <  2:1849900/93‑1 (MQ=255)

                                                    aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCg                                     <  1:4392520/93‑1 (MQ=255)

                                                    aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCg                                     <  2:2929888/93‑1 (MQ=255)

                                                     cGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGt                                    >  2:4060297/1‑93 (MQ=255)

                                                     cGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGt                                    >  1:2270666/1‑93 (MQ=255)

                                                      ggAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                   >  1:4476906/1‑93 (MQ=255)

                                                      ggAACATTTCTGATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGtt                                   <  2:3187401/93‑1 (MQ=255)

                                                       gAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTa                                  >  1:1768488/1‑93 (MQ=255)

                                                       gAACATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTa                                  >  1:5956640/1‑93 (MQ=255)

                                                        aaCATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTaa                                 <  2:5570844/93‑1 (MQ=255)

                                                        aaCATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTaa                                 >  2:2543821/1‑93 (MQ=255)

                                                        aaCATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTaa                                 >  1:3961358/1‑93 (MQ=255)

                                                         aCATTTCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAAc                                <  2:1681305/93‑1 (MQ=255)

                                                            tttCTTATGGGCCAACAGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACttt                             >  1:6308102/1‑93 (MQ=255)

                                                             ttCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTc                            <  1:3786090/93‑1 (MQ=255)

                                                             ttCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTc                            <  1:363527/93‑1 (MQ=255)

                                                             ttCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTc                            <  2:3066723/93‑1 (MQ=255)

                                                              tctcATGGGCCAACGGCATTTCTTACTGCAGTGCCCCCAAACCTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCt                           <  1:1419977/93‑1 (MQ=255)

                                                              tCTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCt                           >  2:6126645/1‑93 (MQ=255)

                                                               cTTATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTg                          <  1:1458327/93‑1 (MQ=255)

                                                                 tATGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCg                        >  1:729865/1‑93 (MQ=255)

                                                                  aTGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGa                       >  1:5879115/1‑93 (MQ=255)

                                                                   tGGGCCAACGGCATTTCTTACTGCGGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                      <  1:6342299/93‑1 (MQ=255)

                                                                   tGGGCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                      >  2:2441670/1‑93 (MQ=255)

                                                                    gggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                     <  1:3339870/93‑1 (MQ=255)

                                                                    gggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                     >  1:5570981/1‑93 (MQ=255)

                                                                    gggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATa                     >  1:1698327/1‑93 (MQ=255)

                                                                     ggCCAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATag                    <  2:2731806/93‑1 (MQ=255)

                                                                       ccAACGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagca                  <  2:1103499/93‑1 (MQ=255)

                                                                         aaCGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagcagc                >  1:824712/1‑93 (MQ=255)

                                                                          aCGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAgcagca               >  2:1546334/1‑93 (MQ=255)

                                                                          aCGGCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAgcagca               >  2:2122361/1‑93 (MQ=255)

                                                                            ggCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGa             <  2:2023829/93‑1 (MQ=255)

                                                                             gCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat            >  2:5162612/1‑93 (MQ=255)

                                                                             gCATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGat            <  1:3746413/93‑1 (MQ=255)

                                                                              cATTTCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata           >  1:5362098/1‑93 (MQ=255)

                                                                                tttCTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg         >  1:2723382/1‑93 (MQ=255)

                                                                                   cTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCa      <  2:3654938/93‑1 (MQ=255)

                                                                                   cTTACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCa      <  2:3311288/93‑1 (MQ=255)

                                                                                     tACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGt    >  2:5789445/1‑93 (MQ=255)

                                                                                     tACTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGt    <  2:2102810/93‑1 (MQ=255)

                                                                                      tcTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTa   <  1:4604263/92‑1 (MQ=255)

                                                                                      aCTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTa   <  2:3891550/93‑1 (MQ=255)

                                                                                      aCTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTa   <  1:6205656/93‑1 (MQ=255)

                                                                                       cTGCAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTaa  <  1:1197301/93‑1 (MQ=255)

GCCCCCGGCAAAATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAA  >  NC_000913/1841832‑1842011

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