BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	236,541	C→T	100%	R159* (CGA→TGA)	dnaQ →	DNA polymerase III epsilon subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	236,541	0	C	T	100.0%	338.6 / NA	103	R159* (CGA→TGA)	dnaQ	DNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (49/54); total (49/54)

Only 100 of 108 total aligned reads displayed.

GGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236444‑236637

ggTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatga                                                                                                >  1:2356121/1‑100 (MQ=255)

 gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCTATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               <  2:2585711/100‑1 (MQ=255)

  tCACCGATAGCATTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaac                                                                                              >  1:2578568/1‑100 (MQ=255)

   cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGAGAACAGTAaatgaacg                                                                                             >  2:3036268/1‑100 (MQ=255)

     ccGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCt                                                                                           <  1:2972588/100‑1 (MQ=255)

     ccGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTACTCGCTACGAAATAGATAACAATAAATGAACGc                                                                                            >  1:598859/1‑99 (MQ=255)

       gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         >  2:2541227/1‑100 (MQ=255)

           gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgg                                                                                     >  2:230429/1‑100 (MQ=255)

            ccTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACggg                                                                                    >  2:3002310/1‑100 (MQ=255)

                gCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAATAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCAt                                                                                <  2:1464317/100‑1 (MQ=255)

                   gTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTAc                                                                             <  1:2728398/100‑1 (MQ=255)

                     ggCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTc                                                                           <  1:2279049/100‑1 (MQ=255)

                     ggCGAGGAAAATGATTGCAGGTAAGAGCAACAGCCTCGATGCGTTATGTGCTCGTTACGAAATAGATAACAGTAAATGAACGCAGCACgggg                                                                                   <  2:2461464/92‑1 (MQ=255)

                      gCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTc                                                                           <  2:1313668/99‑1 (MQ=255)

                      gCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCg                                                                          <  2:660752/100‑1 (MQ=255)

                      gCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCg                                                                          <  2:1853408/100‑1 (MQ=255)

                       cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  2:1315483/1‑100 (MQ=255)

                       cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  2:2660006/1‑100 (MQ=255)

                       cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  1:2671740/1‑100 (MQ=255)

                         aGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATg                                                                       <  1:313166/100‑1 (MQ=255)

                         aGGAAAATGTTTCCCGGTAAGCGAAACAGCCTCGATGCGTTATGTGCTCGCTACGACATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCg                                                                          <  2:2580292/97‑1 (MQ=255)

                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  1:2386488/100‑1 (MQ=255)

                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATAAACGCTGCACGGGGCATTACTCGATGc                                                                      >  1:2090373/1‑100 (MQ=255)

                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGAAAACAGTAAATGAAAGCTGCACGGGGCATTACTCGAtcc                                                                      >  1:46900/1‑98 (MQ=255)

                           gAAAATGTTTTCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATggc                                                                     >  2:605247/1‑98 (MQ=255)

                           gAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  1:205634/99‑1 (MQ=255)

                            aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGAGGCATTACTCGATGccc                                                                    >  1:1204026/1‑100 (MQ=255)

                             aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                   <  1:2207682/100‑1 (MQ=255)

                             aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAAAAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                   >  1:850080/1‑100 (MQ=255)

                              aaTGTTTACCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAg                                                                  >  1:1582346/1‑100 (MQ=255)

                                 gTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGAt                                                                >  2:3010922/1‑99 (MQ=255)

                                  tttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATcc                                                              <  2:806667/100‑1 (MQ=255)

                                     cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           <  1:1075417/100‑1 (MQ=255)

                                      ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  1:2747203/100‑1 (MQ=255)

                                      ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  1:982573/100‑1 (MQ=255)

                                      ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  1:1691439/100‑1 (MQ=255)

                                       cGGTAAGCGCAACAGCCTCGATGTGTTATTTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         <  1:642572/100‑1 (MQ=255)

                                       cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         <  1:775586/100‑1 (MQ=255)

                                        ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         <  1:865403/99‑1 (MQ=255)

                                         gTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGa                                                       <  1:1517409/100‑1 (MQ=255)

                                           aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  2:1264024/100‑1 (MQ=255)

                                           aaGCGCAACAGCCGCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  2:2962538/100‑1 (MQ=255)

                                            aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    <  2:2373410/100‑1 (MQ=255)

                                             gcgcAGCAGCCTCGATGCGTTATGTGCTCGCTGCGAAATAGATAGCAGGAAATGAACGCTGCACGGGGCATTTCTCGATGCCCAGATCCTTGCGGAAGtg                                                   <  2:1569202/100‑2 (MQ=255)

                                             gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  2:2314562/100‑1 (MQ=255)

                                             gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   >  2:131260/1‑100 (MQ=255)

                                             gcgcAACAGCCTCGATGCGTTATGTGATCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCAAGATACTTGCGGAAg                                                     >  2:2719453/1‑98 (MQ=255)

                                             cttccgatctcTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  1:49557/90‑1 (MQ=255)

                                               gcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTa                                                 >  2:141484/1‑100 (MQ=255)

                                                cAACAGTCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTAt                                                >  1:1603998/1‑100 (MQ=255)

                                                cAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTAt                                                >  1:1603897/1‑100 (MQ=255)

                                                 aaCAGCCTCGATGCGTTATGTGCTCGCTACTAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATc                                               <  1:1576257/100‑1 (MQ=255)

                                                 aaCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATc                                               <  1:1537214/100‑1 (MQ=255)

                                                  aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGAAAGTTTATCt                                              <  2:2121966/100‑1 (MQ=255)

                                                  aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCACATCCTTGCGGAAGTTTATCt                                              <  2:1848507/100‑1 (MQ=255)

                                                   cAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGAGGAAGTTTATCTg                                             <  1:2038674/100‑1 (MQ=255)

                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:865916/100‑1 (MQ=255)

                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:1705388/100‑1 (MQ=255)

                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:1657114/100‑1 (MQ=255)

                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  1:2584528/1‑100 (MQ=255)

                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTAGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  2:2818719/1‑100 (MQ=255)

                                                       cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttagatcggaag                                         >  2:49557/1‑90 (MQ=255)

                                                       cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                         >  1:1969859/1‑100 (MQ=255)

                                                          gATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGa                                      <  2:2233851/100‑1 (MQ=255)

                                                            tGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGag                                     >  1:1558424/1‑98 (MQ=255)

                                                            tGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGAcc                                    >  2:1074546/1‑100 (MQ=255)

                                                            tGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGAcc                                    <  1:2224958/100‑1 (MQ=255)

                                                             gCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCg                                   >  2:730387/1‑100 (MQ=255)

                                                              cGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                  <  2:454112/100‑1 (MQ=255)

                                                               gatATGTGCTCGCGACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                 >  1:2045690/3‑100 (MQ=255)

                                                               gTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                 >  2:2261810/1‑100 (MQ=255)

                                                                   tgtgCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCATTGCGGAAGTTTATCTGGCGATGACCGGTGATc                             >  2:2891329/1‑100 (MQ=255)

                                                                    gtgCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCa                            >  2:2901559/1‑100 (MQ=255)

                                                                      gCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                          >  2:118325/1‑100 (MQ=255)

                                                                        tCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGTAAGTTTATCTGGCGATGACCGGTGGTCAAACg                        <  1:2891361/100‑1 (MQ=255)

                                                                         cGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                       >  2:128320/1‑100 (MQ=255)

                                                                         cGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGAAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                       >  1:1893989/1‑100 (MQ=255)

                                                                           cTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCg                     >  2:2086212/1‑100 (MQ=255)

                                                                           cTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCg                     <  1:1704225/100‑1 (MQ=255)

                                                                             aCGAAATAGATAACAGTAAATGAACGTTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGAt                   <  1:495843/100‑1 (MQ=255)

                                                                             aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGAt                   >  2:2782827/1‑100 (MQ=255)

                                                                             aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGAt                   >  2:1317936/1‑100 (MQ=255)

                                                                              cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  2:738465/1‑100 (MQ=255)

                                                                              cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  <  1:2824982/100‑1 (MQ=255)

                                                                               gAAATAGCTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTGTCTGGCGATGACCGGTGGTCACACGTCGATg                  >  2:1051223/1‑99 (MQ=255)

                                                                               gAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg                 >  1:1325066/1‑100 (MQ=255)

                                                                                    aGATAACTGGAATTGAACGCTGCACGGGGTATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCGGGCGATGACCGGTGGTCAAACGTCGATGGCtttt            <  1:1804866/100‑1 (MQ=255)

                                                                                      atctcAGTAAATGAACGCTGCACGGGGCATTACTCGATGCGCAGATTCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          <  1:1728507/96‑1 (MQ=255)

                                                                                      aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          >  1:2231086/1‑100 (MQ=255)

                                                                                      aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGCAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          <  2:1675401/100‑1 (MQ=255)

                                                                                       tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg         <  2:1527373/100‑1 (MQ=255)

                                                                                        aaCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGa        >  2:2405180/1‑100 (MQ=255)

                                                                                         aCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt       <  2:2551906/100‑1 (MQ=255)

                                                                                         aCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt       >  2:173238/1‑100 (MQ=255)

                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGcagat      >  2:1728487/1‑96 (MQ=255)

                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:2517813/1‑100 (MQ=255)

                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACGCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      <  2:1539667/100‑1 (MQ=255)

                                                                                           aGTAAATGAACGCTGCACGTTGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     <  1:1810138/100‑1 (MQ=255)

                                                                                           aGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     >  1:2987780/1‑100 (MQ=255)

                                                                                              aaaTGAACGCTGCCCGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAg  <  2:2323168/100‑1 (MQ=255)

GGTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236444‑236637

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