BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,012	T→C	I116V (ATT→GTT)	b3691 ←	D‑galactonate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,012	0	T	C	99.8%	1334.9 / ‑5.6	417	I116V (ATT→GTT)	b3691	D‑galactonate transport
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (215/201); total (215/202)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 434 total aligned reads displayed.

TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101

tataAAAACCAACGTCAGGAGCGGGTTCATGTTCCGGGAACCAGCTGTTCACCATACGCTTATTGGTCGGCACAGCAGGCGCTT‑GAAAAtcac                                                                                        <  2:1226744/93‑4 (MQ=255)

  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       >  1:5335571/1‑93 (MQ=255)

    aaaaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGt                                                                                     <  1:6885264/93‑1 (MQ=255)

     aaaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGtt                                                                                    >  1:2099914/1‑93 (MQ=255)

      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   >  2:5909394/1‑93 (MQ=255)

       aCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTAt                                                                                  <  2:4041767/93‑1 (MQ=255)

       aCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTAt                                                                                  >  1:9971015/1‑93 (MQ=255)

         cAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcg                                                                                <  1:3842873/93‑1 (MQ=255)

           aCGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcgcg                                                                              <  1:6191464/93‑1 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCgcgcgc                                                                           >  1:3176321/1‑93 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCgcgcgc                                                                           <  2:2456940/93‑1 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCgcgcgc                                                                           >  1:4403531/1‑93 (MQ=255)

              gCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGGCACCATCCGGTTATTGGTCGGGAGAGAAGGCGCTTCGACAACACCGGTTATCgagcgc                                                                           >  1:958221/1‑93 (MQ=255)

               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          >  1:890225/1‑93 (MQ=255)

               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  1:8580615/93‑1 (MQ=255)

               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  1:7530258/93‑1 (MQ=255)

               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  2:8508736/93‑1 (MQ=255)

               cAGAAGCGCGTTCATGTGCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  2:1343944/93‑1 (MQ=255)

                aGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAg                                                                         >  1:9187911/1‑93 (MQ=255)

                    gcgcGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCa                                                                     >  1:2062012/1‑93 (MQ=255)

                    gcgcGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCa                                                                     <  2:4709116/93‑1 (MQ=255)

                      gcgTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAAt                                                                   <  1:4908010/93‑1 (MQ=255)

                      gcgGTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAAt                                                                   <  2:3278665/93‑1 (MQ=255)

                        gTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTa                                                                 >  1:5232852/1‑93 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  2:961489/1‑93 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  1:3692785/1‑93 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                <  1:5641834/93‑1 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  2:2663961/1‑93 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  2:169854/1‑93 (MQ=255)

                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  1:10031350/1‑93 (MQ=255)

                           cATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATg                                                              >  2:7163616/1‑93 (MQ=255)

                            aTGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGa                                                             >  1:25550/1‑93 (MQ=255)

                            aTGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGa                                                             >  1:3534838/1‑93 (MQ=255)

                             tGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGAc                                                            <  2:1088828/93‑1 (MQ=255)

                             agTTCCGGCAACCAGCTGGTCACCATCCGTTTAGTGGTCTGGAAGGCAGCCGCTTCGAAACCACCGCTTATCGCGCGCAGACCAATTAATGAc                                                            <  2:3829454/92‑1 (MQ=255)

                               ttCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAt                                                          >  1:4485683/1‑93 (MQ=255)

                               ttCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAt                                                          <  1:5813591/93‑1 (MQ=255)

                                tCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAtt                                                         >  2:2877821/1‑93 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  1:659534/93‑1 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       >  2:9009153/1‑93 (MQ=255)

                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  2:759127/93‑1 (MQ=255)

                                   gggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAg                                                      >  2:2023809/1‑93 (MQ=255)

                                   gggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAg                                                      >  1:4958/1‑93 (MQ=255)

                                   gggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAg                                                      >  1:1321564/1‑93 (MQ=255)

                                    ggAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGc                                                     >  2:1877499/1‑93 (MQ=255)

                                     gAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGcc                                                    <  1:7237435/93‑1 (MQ=255)

                                      aaCCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGccc                                                   <  2:6402632/93‑1 (MQ=255)

                                        ccAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGt                                                 >  2:7044273/1‑93 (MQ=255)

                                        ccAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGt                                                 <  2:7398598/93‑1 (MQ=255)

                                        ccAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGt                                                 <  1:3425035/93‑1 (MQ=255)

                                         cAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTg                                                <  2:5206548/93‑1 (MQ=255)

                                         cAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTg                                                <  1:1593359/93‑1 (MQ=255)

                                           gCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGc                                              >  2:8795124/1‑93 (MQ=255)

                                            cTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCa                                             >  1:2666007/1‑93 (MQ=255)

                                             tGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCaa                                            <  1:6700922/93‑1 (MQ=255)

                                              ggTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCaaa                                           <  2:5432819/93‑1 (MQ=255)

                                              ggTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCaaa                                           >  1:4629243/1‑93 (MQ=255)

                                                 cACCATCCTGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGcc                                        <  2:4630678/93‑1 (MQ=255)

                                                    cATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTg                                     >  2:3427955/1‑93 (MQ=255)

                                                     aTCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTgg                                    <  2:9342206/93‑1 (MQ=255)

                                                        cGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAAt                                 >  2:671745/1‑93 (MQ=255)

                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              <  1:9945494/93‑1 (MQ=255)

                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              <  1:6251939/93‑1 (MQ=255)

                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              >  2:7839904/1‑93 (MQ=255)

                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              >  1:2461825/1‑93 (MQ=255)

                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              <  2:4430810/93‑1 (MQ=255)

                                                              ttGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTc                           <  1:4096002/93‑1 (MQ=255)

                                                                tgTCGGGAAAGCAGGCGCTTCGAAAACACCGGTGATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGc                         <  1:9722689/92‑1 (MQ=255)

                                                                ggTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGc                         <  2:1512543/93‑1 (MQ=255)

                                                                   cGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCAc                      <  2:9757520/93‑1 (MQ=255)

                                                                    gggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACt                     <  2:9408258/93‑1 (MQ=255)

                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTg                    >  2:3043030/1‑93 (MQ=255)

                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTg                    >  1:6300657/1‑93 (MQ=255)

                                                                      gAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGa                   <  2:8330517/93‑1 (MQ=255)

                                                                       aaaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAc                  <  2:80750/93‑1 (MQ=255)

                                                                        aaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAcc                 >  2:3880583/1‑93 (MQ=255)

                                                                        aaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAcc                 <  1:5475837/93‑1 (MQ=255)

                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                >  2:3641322/1‑93 (MQ=255)

                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                <  1:1904175/93‑1 (MQ=255)

                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGACCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                >  2:2839384/1‑93 (MQ=255)

                                                                          gCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAg               <  2:2278865/93‑1 (MQ=255)

                                                                          gCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAg               >  2:6134432/1‑93 (MQ=255)

                                                                          gCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCAACTGACCAg               >  1:3269129/1‑93 (MQ=255)

                                                                           cAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGc              >  1:8683316/1‑93 (MQ=255)

                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGAGATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  1:7544328/1‑93 (MQ=255)

                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            <  2:9260845/93‑1 (MQ=255)

                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  2:6273395/1‑93 (MQ=255)

                                                                              gcgcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCaa           <  2:6615372/93‑1 (MQ=255)

                                                                               cgcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCaaa          <  1:300789/93‑1 (MQ=255)

                                                                                gcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAat         >  2:5032999/1‑93 (MQ=255)

                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCCAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        >  1:8343489/1‑93 (MQ=255)

                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        >  2:7962185/1‑93 (MQ=255)

                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        <  2:1026374/93‑1 (MQ=255)

                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       >  1:6037773/1‑93 (MQ=255)

                                                                                   tCGAAAACACCGGTTATCGCGCGCAGACCAATTTATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATc      <  1:7835169/93‑1 (MQ=255)

                                                                                   tCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATc      >  1:8143329/1‑93 (MQ=255)

                                                                                   tCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATc      >  1:4784987/1‑93 (MQ=255)

                                                                                     gAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGc    >  2:6924392/1‑93 (MQ=255)

                                                                                     gAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGc    >  2:1907657/1‑93 (MQ=255)

                                                                                       aaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCaa  <  1:2885380/93‑1 (MQ=255)

TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101

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