BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,037,676	(C)_5→6	coding (177/1005 nt)	b1971 →	putative reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,037,671	1	.	C	98.6%	462.6 / ‑1.5	147	coding (172/1005 nt)	b1971	putative reductase
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (82/63); total (83/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 188 total aligned reads displayed.

AAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTT  >  NC_000913/2037579‑2037758

aaccACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATcgcccgcccg                                                                                          <  1:284951/90‑1 (MQ=255)

  gCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATcgcccgcccg‑cc                                                                                       <  2:3519585/93‑1 (MQ=255)

    aCTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑cccc                                                                                     <  1:89733/93‑1 (MQ=255)

     cTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑ccccc                                                                                    <  2:4211614/93‑1 (MQ=255)

     cTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccg‑ccccc                                                                                    <  1:2981902/93‑1 (MQ=255)

      tGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTANCTNGTTTAAAGGGCACGATCGCCCGcccg‑cccccc                                                                                   >  2:4833413/1‑92 (MQ=255)

       gggCATCAGTGGAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGATTAAAGGGAACGATCGCCCGcccg‑cgccctg                                                                                  >  2:3047275/1‑91 (MQ=255)

        ggCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGCGAACGATCGCCCGcccg‑ccccccgc                                                                                 >  1:266969/1‑90 (MQ=255)

           aTCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCgg                                                                               >  1:1899555/1‑93 (MQ=255)

           aTCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCgg                                                                               >  2:5754906/1‑93 (MQ=255)

            tCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGa                                                                              >  2:4318784/1‑93 (MQ=255)

            tCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGa                                                                              >  1:3616566/1‑93 (MQ=255)

              aGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaa                                                                            >  1:4538264/1‑93 (MQ=255)

               gTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGaaaa                                                                           <  1:1447476/93‑1 (MQ=255)

                tGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAg                                                                          <  1:2248629/93‑1 (MQ=255)

                  cAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAgcg                                                                        <  1:3244513/93‑1 (MQ=255)

                   aaCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAgcgc                                                                       <  1:1847083/93‑1 (MQ=255)

                   aaCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAgcgc                                                                       <  1:1681161/93‑1 (MQ=255)

                   aaCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAgcgc                                                                       >  1:5657150/1‑93 (MQ=255)

                    aCTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCt                                                                      >  1:430532/1‑93 (MQ=255)

                     cTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTg                                                                     >  2:3694277/1‑93 (MQ=255)

                     cTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCCCTg                                                                     >  2:3923318/1‑93 (MQ=255)

                      tGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTgg                                                                    >  2:2416693/1‑93 (MQ=255)

                       gCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGa                                                                   <  2:5032928/93‑1 (MQ=255)

                       gCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGa                                                                   <  2:1391069/93‑1 (MQ=255)

                       gCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGa                                                                   >  1:5737023/1‑93 (MQ=255)

                         aCTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGt                                                                 >  1:4979843/1‑93 (MQ=255)

                          ctttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCTCCCGCCCCCCGCCGGAAAAGCGCTGGAGtt                                                                >  1:1352493/1‑93 (MQ=255)

                          ctttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGtt                                                                >  2:2456845/1‑93 (MQ=255)

                           tttctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTc                                                               <  1:2045513/93‑1 (MQ=255)

                             tctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAg                                                             <  1:5703381/93‑1 (MQ=255)

                              ctttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGc                                                            <  1:2896301/93‑1 (MQ=255)

                               tttGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCa                                                           <  2:6294915/93‑1 (MQ=255)

                                  gCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGc                                                        <  2:5662547/93‑1 (MQ=255)

                                  gCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGc                                                        >  2:3286016/1‑93 (MQ=255)

                                   ccTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGcc                                                       <  1:1993634/93‑1 (MQ=255)

                                    cTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCt                                                      >  1:4623238/1‑93 (MQ=255)

                                     tCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTg                                                     >  1:2446548/1‑93 (MQ=255)

                                     tCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTg                                                     >  1:3300147/1‑93 (MQ=255)

                                       aCGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgcc                                                   >  1:2110267/1‑93 (MQ=255)

                                         gCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgccgc                                                 <  1:5622389/93‑1 (MQ=255)

                                          cTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgccgcc                                                >  2:725923/1‑93 (MQ=255)

                                            gcgcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                              >  1:1217195/1‑93 (MQ=255)

                                            gcgcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                              <  2:3940897/93‑1 (MQ=255)

                                            gcgcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                              <  2:1365065/93‑1 (MQ=255)

                                              gcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGc                                            >  1:4594717/1‑93 (MQ=255)

                                              gcATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGc                                            >  2:4581826/1‑93 (MQ=255)

                                                aTGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaa                                          >  1:1214376/1‑93 (MQ=255)

                                                 tGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaaa                                         >  1:3686851/1‑93 (MQ=255)

                                                 tGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaaa                                         >  1:6166792/1‑93 (MQ=255)

                                                  gCCTATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaaaa                                        <  1:917429/93‑1 (MQ=255)

                                                   ccGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAAt                                       <  2:4436187/93‑1 (MQ=255)

                                                    cGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATa                                      <  1:4853718/93‑1 (MQ=255)

                                                    cGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATa                                      >  2:472974/1‑93 (MQ=255)

                                                      aTCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAc                                    >  2:2109227/1‑93 (MQ=255)

                                                      aTCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAc                                    >  2:1736321/1‑93 (MQ=255)

                                                       tCTGCTTAGCTGTTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCCCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAcc                                   <  2:5371939/93‑1 (MQ=255)

                                                       tCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAct                                   >  2:4162909/1‑92 (MQ=255)

                                                       tCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAcc                                   >  2:2600602/1‑93 (MQ=255)

                                                       tCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAAcc                                   <  1:5854659/93‑1 (MQ=255)

                                                        cTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCt                                  <  2:922536/93‑1 (MQ=255)

                                                         tGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTg                                 >  2:5822275/1‑93 (MQ=255)

                                                         tGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTg                                 >  2:520410/1‑93 (MQ=255)

                                                         tGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTg                                 <  2:3815563/93‑1 (MQ=255)

                                                          gCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGc                                >  2:5746427/1‑93 (MQ=255)

                                                           cTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGcc                               <  1:2589160/93‑1 (MQ=255)

                                                           cTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGcc                               >  2:4793549/1‑93 (MQ=255)

                                                            ttAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCa                              <  2:2789802/93‑1 (MQ=255)

                                                            ttAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCa                              <  2:6135318/93‑1 (MQ=255)

                                                             tAGCTGGTTTAAAGGAAACGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCAc                             <  1:726731/93‑1 (MQ=255)

                                                               gCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCATCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTg                           >  1:6282290/1‑93 (MQ=255)

                                                                cgggTTTAAGGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGCAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGa                          <  1:3371871/91‑1 (MQ=255)

                                                                  ggTTTAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACg                        <  1:6038827/93‑1 (MQ=255)

                                                                   gTTTAAAGGAAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGc                       <  1:3810316/93‑1 (MQ=255)

                                                                    tttAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGcc                      >  2:845896/1‑93 (MQ=255)

                                                                    tttAAAGGGAACGACCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGcc                      <  2:4616822/93‑1 (MQ=255)

                                                                     ttAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCca                     >  1:5446746/1‑93 (MQ=255)

                                                                     ttAAAGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCca                     <  2:2683058/93‑1 (MQ=255)

                                                                        aaGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagca                  <  1:3556588/93‑1 (MQ=255)

                                                                         aGGGAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagcag                 >  1:5346347/1‑93 (MQ=255)

                                                                          gggAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGa                >  2:3168374/1‑93 (MQ=255)

                                                                           ggACCGAGCCCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGAt               <  1:2552795/93‑1 (MQ=255)

                                                                           ggAACGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGAt               >  1:2565315/1‑93 (MQ=255)

                                                                              aCGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATaaa            <  1:4160091/93‑1 (MQ=255)

                                                                              aCGATCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATaaa            <  2:5974695/93‑1 (MQ=255)

                                                                                gaacgccctcccgccCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGt          <  1:5172664/90‑1 (MQ=255)

                                                                                 aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCGGCAGATAAAGtt         >  1:5595312/1‑93 (MQ=255)

                                                                                 aTCGCCCGCCCGCCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGtt         >  1:2394668/1‑93 (MQ=255)

                                                                                   cccccccccc‑ccccCGCCGGAAAGGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      <  1:5835197/91‑1 (MQ=255)

                                                                                    gcccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      <  1:5228014/93‑1 (MQ=255)

                                                                                    gcccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      >  1:6070648/1‑93 (MQ=255)

                                                                                    ccccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTcc      <  1:1904420/92‑1 (MQ=255)

                                                                                     cccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCTCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCg     >  2:1796871/1‑93 (MQ=255)

                                                                                     cccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCg     <  1:1475135/93‑1 (MQ=255)

                                                                                     cccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCg     >  1:1234128/1‑93 (MQ=255)

                                                                                     cccgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCg     >  2:104988/1‑93 (MQ=255)

                                                                                       cgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGt   >  2:5686469/1‑93 (MQ=255)

                                                                                       cgcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGt   >  1:98413/1‑93 (MQ=255)

                                                                                        gcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGtt  >  1:4596704/1‑93 (MQ=255)

                                                                                        gcccgcccCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGtt  >  2:2045413/1‑93 (MQ=255)

AAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCG‑CCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTT  >  NC_000913/2037579‑2037758

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