BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	393,476	G→A	G413D (GGC→GAC)	yaiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	393,476	0	G	A	99.3%	526.9 / NA	136	G413D (GGC→GAC)	yaiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (79/56); total (79/57)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

TCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCT  >  AC_000091/393389‑393562

tCACCTACAGCGAAGATCCGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAacgacg                                                                                      <  1:1017609/90‑1 (MQ=255)

 cACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAAcgacga                                                                                     <  1:576033/90‑1 (MQ=255)

  aCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGAt                                                                                    >  1:3119604/1‑90 (MQ=255)

  aCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGAt                                                                                    <  1:2132565/90‑1 (MQ=255)

  aCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGAt                                                                                    >  2:1309120/1‑90 (MQ=255)

    cTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTc                                                                                  >  1:3334813/1‑90 (MQ=255)

     tACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTcc                                                                                 <  1:3044507/90‑1 (MQ=255)

     tACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTcc                                                                                 >  2:399730/1‑90 (MQ=255)

       cAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCAt                                                                               >  2:604852/1‑90 (MQ=255)

       cAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCAt                                                                               >  1:778611/1‑90 (MQ=255)

        aGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATc                                                                              >  2:2471521/1‑90 (MQ=255)

        aGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATc                                                                              <  1:1808876/90‑1 (MQ=255)

        aGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATc                                                                              >  2:1038921/1‑90 (MQ=255)

           gAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGAt                                                                           >  2:2479832/1‑90 (MQ=255)

           gAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGAt                                                                           >  2:919886/1‑90 (MQ=255)

           gAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGAt                                                                           <  1:79770/90‑1 (MQ=255)

                tCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACgg                                                                      <  1:1706814/90‑1 (MQ=255)

                tCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACgg                                                                      <  1:640088/90‑1 (MQ=255)

                 cAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGt                                                                     >  2:1936795/1‑90 (MQ=255)

                   ggCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTAc                                                                   <  1:1767721/90‑1 (MQ=255)

                   ggCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTAc                                                                   >  1:3073712/1‑90 (MQ=255)

                    gCTCTGACTCCGTACTTCAATCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACt                                                                  <  1:1114955/90‑1 (MQ=255)

                    gCTCTGACTCCGTACTTCAAACTGGCTTGCGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACt                                                                  >  2:1179910/1‑90 (MQ=255)

                    gCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACt                                                                  <  2:731791/90‑1 (MQ=255)

                    gCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACt                                                                  <  1:243916/90‑1 (MQ=255)

                    gCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACt                                                                  <  2:3244803/90‑1 (MQ=255)

                      tctGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGa                                                                >  1:627308/1‑90 (MQ=255)

                      tctGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGa                                                                >  1:2702775/1‑90 (MQ=255)

                        tGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAg                                                              <  2:931349/90‑1 (MQ=255)

                          aCTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAggg                                                            <  2:2541290/90‑1 (MQ=255)

                          aCTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAggg                                                            >  1:254333/1‑90 (MQ=255)

                           cTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGt                                                           >  1:3482673/1‑90 (MQ=255)

                             ccGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCt                                                         <  1:301284/90‑1 (MQ=255)

                             ccGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCt                                                         <  1:1146326/90‑1 (MQ=255)

                              cGTACTTCAAACTGGCTTGCGTCTACGACAACTCTAACAACGATGACGATGTGTACGACGATTCCATCGATAACGGTACTGAAGGGTCTg                                                        >  2:877086/1‑90 (MQ=255)

                              cGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTg                                                        <  1:3421828/90‑1 (MQ=255)

                                tACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCg                                                      >  1:1308951/1‑90 (MQ=255)

                                tACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCg                                                      <  2:2767375/90‑1 (MQ=255)

                                  cTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGt                                                    <  1:1713049/90‑1 (MQ=255)

                                  cTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGt                                                    >  1:3441230/1‑90 (MQ=255)

                                   ttCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa                                                   <  1:2261817/90‑1 (MQ=255)

                                   ttCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa                                                   <  1:1298224/90‑1 (MQ=255)

                                    tCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTAc                                                  >  2:3355722/1‑90 (MQ=255)

                                    tCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTAc                                                  <  1:1550540/90‑1 (MQ=255)

                                        aCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGAGGTACgtgt                                              >  2:2819137/1‑90 (MQ=255)

                                         cTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGtt                                             >  2:2216555/1‑90 (MQ=255)

                                         cTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGtt                                             <  1:1859454/90‑1 (MQ=255)

                                         cTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGtt                                             >  1:3116710/1‑90 (MQ=255)

                                          tGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTg                                            >  1:1926992/1‑90 (MQ=255)

                                          tGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTg                                            <  1:1035205/90‑1 (MQ=255)

                                          tGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTg                                            >  2:1665836/1‑90 (MQ=255)

                                           ggCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTgg                                           >  2:2358885/1‑90 (MQ=255)

                                            gCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGt                                          >  2:3109987/1‑90 (MQ=255)

                                             cTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTc                                         <  1:935782/90‑1 (MQ=255)

                                               tACGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTg                                       <  1:703211/90‑1 (MQ=255)

                                                aCGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTgg                                      <  2:245221/90‑1 (MQ=255)

                                                aCGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTgg                                      >  1:3054721/1‑90 (MQ=255)

                                                 cGTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTggg                                     >  1:3235989/1‑90 (MQ=255)

                                                  gTCTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGt                                    >  2:157312/1‑90 (MQ=255)

                                                    cTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTAc                                  >  2:1365227/1‑90 (MQ=255)

                                                    cTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTAc                                  >  2:3401997/1‑90 (MQ=255)

                                                    cTACGACGACTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTAc                                  <  1:3111508/90‑1 (MQ=255)

                                                       cgacgaCTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCa                               >  2:2667839/1‑90 (MQ=255)

                                                        ggcgacTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAg                              >  2:1750569/3‑90 (MQ=255)

                                                         acgacTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGt                             >  2:664455/1‑90 (MQ=255)

                                                          cgacTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGtt                            >  1:1338223/1‑90 (MQ=255)

                                                           gacTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGttt                           >  2:3141262/1‑90 (MQ=255)

                                                            acTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTa                          <  2:2088028/90‑1 (MQ=255)

                                                            acTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTa                          >  1:2638079/1‑90 (MQ=255)

                                                            acTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTa                          <  2:1352689/90‑1 (MQ=255)

                                                            acTCTAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTa                          >  2:2469159/1‑90 (MQ=255)

                                                             ctctAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACCTGTTGGTCTGGGTACTCAGTTTAg                         >  2:3305260/1‑90 (MQ=255)

                                                               ctAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCt                       >  2:1175097/1‑90 (MQ=255)

                                                               ctAACAACGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCt                       <  2:2695281/90‑1 (MQ=255)

                                                                 agcaacGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTc                     >  2:2101281/3‑90 (MQ=255)

                                                                 aacaacGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTc                     >  2:2023754/1‑90 (MQ=255)

                                                                 aacaacGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTc                     >  1:1824445/1‑90 (MQ=255)

                                                                  acaacGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCa                    >  1:680564/1‑90 (MQ=255)

                                                                   caacGATAACGATGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCAc                   <  1:2174004/90‑1 (MQ=255)

                                                                     acgataacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCa                 >  2:434895/1‑90 (MQ=255)

                                                                     acgataacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCa                 <  2:2977318/90‑1 (MQ=255)

                                                                     acgataacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCa                 <  2:2664525/90‑1 (MQ=255)

                                                                       gataacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAg               >  1:401340/1‑90 (MQ=255)

                                                                        ataacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGa              >  1:1189145/1‑90 (MQ=255)

                                                                          aacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAAc            >  2:1230720/1‑90 (MQ=255)

                                                                          aacgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAAc            >  2:15633/1‑90 (MQ=255)

                                                                            cgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACtt          >  2:1658350/1‑90 (MQ=255)

                                                                            cgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACtt          <  2:1005389/90‑1 (MQ=255)

                                                                            cgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACtt          <  1:2380827/90‑1 (MQ=255)

                                                                            cgatGTGAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACtt          >  2:2169102/1‑90 (MQ=255)

                                                                               tgtgAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAg       >  1:1997150/1‑90 (MQ=255)

                                                                                gtgAACGACGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgc      <  2:2886231/90‑1 (MQ=255)

                                                                                gtgAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgc      <  1:1225302/90‑1 (MQ=255)

                                                                                 tgcacgacgATTACATCGATAACGGAACTGAAGGGTCTGCGGTACGTGTTGGTCAGGGTACTCAGTTTAGATTCACCAAGAACTTCAgcg     >  2:1242638/4‑90 (MQ=255)

                                                                                 tgAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgcg     <  2:2933935/90‑1 (MQ=255)

                                                                                 tgAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgcg     >  2:1414400/1‑90 (MQ=255)

                                                                                 tgAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgcg     <  1:1011394/90‑1 (MQ=255)

                                                                                 ggAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgcg     <  2:667348/89‑1 (MQ=255)

                                                                                  gAACGACGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAgcgc    >  1:3486279/1‑90 (MQ=255)

                                                                                    acgacgATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGCTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCt  >  1:2026839/1‑90 (MQ=255)

TCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCT  >  AC_000091/393389‑393562

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