BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,211,165	A→G	I21T (ATA→ACA)	ycgW ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,211,165	0	A	G	86.4%	337.3 / 47.1	140	I21T (ATA→ACA)	ycgW	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/16); new base G (27/94); total (30/110)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 145 total aligned reads displayed.

GGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAAT  >  NC_000913/1211073‑1211255

 gATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCatat                                                                                           <  1:2189610/93‑1 (MQ=255)

  ataataCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCAtgtg                                                                                          >  1:480662/1‑93 (MQ=255)

    aataCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGc                                                                                        >  1:4139529/1‑93 (MQ=255)

    aataCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGc                                                                                        <  1:5509138/93‑1 (MQ=255)

    aataCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGc                                                                                        <  2:4879239/93‑1 (MQ=255)

      taCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTg                                                                                      >  1:1460901/1‑93 (MQ=255)

      taCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTg                                                                                      <  1:3816414/93‑1 (MQ=255)

       aCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGa                                                                                     <  1:3928301/93‑1 (MQ=255)

        ccAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGaa                                                                                    >  2:4063271/1‑93 (MQ=255)

          aGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGaaaa                                                                                  <  1:3026970/93‑1 (MQ=255)

          aGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGaaaa                                                                                  <  1:3504884/93‑1 (MQ=255)

           gACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGaaaaa                                                                                 >  1:1089220/1‑93 (MQ=255)

           gACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGaaaaa                                                                                 <  1:2661431/93‑1 (MQ=255)

           gACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGaaaaa                                                                                 <  1:4624107/93‑1 (MQ=255)

              ttAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGat                                                                              <  2:7392092/93‑1 (MQ=255)

               tAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGata                                                                             <  1:4225835/93‑1 (MQ=255)

               tAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGata                                                                             <  2:1676992/93‑1 (MQ=255)

               tAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGata                                                                             <  2:3532744/93‑1 (MQ=255)

                 aCCGGTGTAAATATACTGCCAGGCGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGAtata                                                                           <  2:6970453/93‑1 (MQ=255)

                 aCCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGAtata                                                                           >  1:6058466/1‑93 (MQ=255)

                 aCCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGAtata                                                                           >  1:2321997/1‑93 (MQ=255)

                  ccGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATAc                                                                          <  2:5545962/93‑1 (MQ=255)

                   cGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATAcc                                                                         >  1:4091128/1‑93 (MQ=255)

                   cGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATAcc                                                                         <  2:1688960/93‑1 (MQ=255)

                    ggTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACca                                                                        >  2:399475/1‑93 (MQ=255)

                    ggTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACca                                                                        <  1:3157043/93‑1 (MQ=255)

                    ggTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACca                                                                        >  1:5791981/1‑93 (MQ=255)

                     gtgtAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACcac                                                                       <  2:1750418/93‑1 (MQ=255)

                      tgtAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACcaca                                                                      <  2:6084034/93‑1 (MQ=255)

                       gtAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACAt                                                                     <  2:2650058/93‑1 (MQ=255)

                         aaaTATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGt                                                                   <  1:210375/93‑1 (MQ=255)

                          aaTATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTa                                                                  <  1:5249217/93‑1 (MQ=255)

                           atatACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAg                                                                 <  1:6274390/93‑1 (MQ=255)

                           atatACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAg                                                                 <  2:3143089/93‑1 (MQ=255)

                            tataCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAgg                                                                >  2:786242/1‑93 (MQ=255)

                             ataCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGt                                                               <  2:1116233/93‑1 (MQ=255)

                             ataCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGt                                                               <  2:824591/93‑1 (MQ=255)

                             ataCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGt                                                               <  1:586872/93‑1 (MQ=255)

                              taCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGtt                                                              <  2:5576108/93‑1 (MQ=255)

                              taCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGtt                                                              <  2:6390525/93‑1 (MQ=255)

                              taCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGtt                                                              <  1:5040522/93‑1 (MQ=255)

                                 tGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGaa                                                           <  2:3440073/93‑1 (MQ=255)

                                 tGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGaa                                                           <  1:4216445/93‑1 (MQ=255)

                                   ccAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAAtt                                                         <  2:2972840/93‑1 (MQ=255)

                                    cAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTACATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTa                                                        >  2:6566987/1‑93 (MQ=255)

                                     aggaggGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTAc                                                       <  2:4549607/93‑1 (MQ=255)

                                     aggaggGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTAc                                                       <  1:3852058/93‑1 (MQ=255)

                                      ggagggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTAcc                                                      <  1:7115271/93‑1 (MQ=255)

                                      ggagggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTAcc                                                      <  1:6435057/93‑1 (MQ=255)

                                      ggagggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTAcc                                                      >  1:1802026/1‑93 (MQ=255)

                                       gagggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCg                                                     <  2:3821968/93‑1 (MQ=255)

                                        agggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCgt                                                    <  1:1551139/93‑1 (MQ=255)

                                         gggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCgtg                                                   <  2:7401564/93‑1 (MQ=255)

                                           gagaAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTc                                                 <  1:1016000/93‑1 (MQ=255)

                                           gagaAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTc                                                 >  2:2380636/1‑93 (MQ=255)

                                              aaaTATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAAt                                              <  2:4207399/93‑1 (MQ=255)

                                                atatAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTa                                            >  2:7472606/1‑93 (MQ=255)

                                                atatAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTa                                            <  1:2634019/93‑1 (MQ=255)

                                                 tataGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTAc                                           <  1:5248079/93‑1 (MQ=255)

                                                 tataGTAGATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTAc                                           <  1:4782511/93‑1 (MQ=255)

                                                   tagtagATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTa                                         <  2:1959952/93‑1 (MQ=255)

                                                   tagtagATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTa                                         <  1:7014223/93‑1 (MQ=255)

                                                     gtagATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATc                                       <  2:3259800/93‑1 (MQ=255)

                                                      tagATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATcc                                      <  2:3304712/93‑1 (MQ=255)

                                                        gATTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCAc                                    <  2:7427726/93‑1 (MQ=255)

                                                         aTTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACt                                   <  1:437626/93‑1 (MQ=255)

                                                         aTTGATACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACt                                   <  2:4932851/93‑1 (MQ=255)

                                                             aTACCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCAt                               <  2:1620759/93‑1 (MQ=255)

                                                               aCCAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCAttt                             <  2:2842822/93‑1 (MQ=255)

                                                                ccAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTg                            >  1:559809/1‑93 (MQ=255)

                                                                 cAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCGTTTGt                           <  2:1285929/93‑1 (MQ=255)

                                                                 cAGATGATCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGt                           <  2:3334358/93‑1 (MQ=255)

                                                                   gatgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTa                         <  2:3516805/93‑1 (MQ=255)

                                                                   gatgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTa                         <  2:3133430/93‑1 (MQ=255)

                                                                   gatgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTa                         <  2:4494529/93‑1 (MQ=255)

                                                                   gatgatCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTCGTTa                         <  2:1673743/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTAt                        <  2:6568730/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTAt                        <  1:4505626/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTAt                        <  2:3139838/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTAt                        <  2:3758988/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATCTGTTAt                        <  2:6988427/93‑1 (MQ=255)

                                                                    atgatCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTAt                        <  1:2288268/93‑1 (MQ=255)

                                                                     tgatCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATg                       <  2:6962721/93‑1 (MQ=255)

                                                                       atCATTTTCATATTACCCCATGTGGCTTAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTc                     >  2:4854347/1‑93 (MQ=255)

                                                                       atCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTc                     >  1:5624573/1‑93 (MQ=255)

                                                                       atCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTc                     <  1:6482735/93‑1 (MQ=255)

                                                                       atCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTc                     <  1:6629971/93‑1 (MQ=255)

                                                                        tCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCt                    <  2:3793300/93‑1 (MQ=255)

                                                                        tCATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCt                    <  2:1167045/93‑1 (MQ=255)

                                                                        tCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCt                    <  2:2792862/93‑1 (MQ=255)

                                                                        tCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCt                    <  2:14608/93‑1 (MQ=255)

                                                                         cATTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCtt                   >  1:2348388/1‑93 (MQ=255)

                                                                          aTTTTCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATATCTTa                  <  1:1003213/93‑1 (MQ=255)

                                                                           ttttCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTAt                 <  2:6430763/93‑1 (MQ=255)

                                                                            tttCATATTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATc                <  2:3172193/93‑1 (MQ=255)

                                                                                 tatTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACg           <  2:399639/93‑1 (MQ=255)

                                                                                  atTACCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACgg          <  1:1505399/93‑1 (MQ=255)

                                                                                     aCCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTAt       <  2:3602065/93‑1 (MQ=255)

                                                                                     aCCCCATGTGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTAt       <  1:616732/93‑1 (MQ=255)

                                                                                         cATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTaa   >  2:5699256/1‑93 (MQ=255)

GGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAAT  >  NC_000913/1211073‑1211255

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