BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,575,878	A→G	Y344H (TAT→CAT)	mcrB ←	McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,575,878	0	A	G	99.4%	492.2 / ‑4.7	155	Y344H (TAT→CAT)	mcrB	McrB
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (58/96); total (58/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

Only 100 of 156 total aligned reads displayed.

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTT  >  NC_000913/4575788‑4575966

 tGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTc                                                                                       >  1:3893051/1‑93 (MQ=255)

  gCTTTTTTATTGAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTca                                                                                      >  1:4335098/1‑93 (MQ=255)

    ttttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaac                                                                                    <  1:4460496/93‑1 (MQ=255)

    ttttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaac                                                                                    >  2:458786/1‑93 (MQ=255)

    ttttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaaa                                                                                    >  1:1462667/1‑92 (MQ=255)

     tttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaaca                                                                                   >  1:4654326/1‑93 (MQ=255)

     tttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaaca                                                                                   >  1:3288443/1‑93 (MQ=255)

     tttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaaca                                                                                   <  1:2088150/93‑1 (MQ=255)

      ttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaacaa                                                                                  <  2:157792/93‑1 (MQ=255)

      ttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaacaa                                                                                  <  1:1164375/93‑1 (MQ=255)

       tttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCaacaac                                                                                 <  2:499872/93‑1 (MQ=255)

       tttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCaacaac                                                                                 <  1:91750/93‑1 (MQ=255)

       tttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCaacaac                                                                                 <  2:2627583/93‑1 (MQ=255)

        ttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACg                                                                                >  1:552544/1‑93 (MQ=255)

        ttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACg                                                                                <  2:255636/93‑1 (MQ=255)

          aTTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGc                                                                              >  2:3855232/1‑93 (MQ=255)

           ttCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGcc                                                                             >  1:2830485/1‑93 (MQ=255)

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                          <  2:1930212/93‑1 (MQ=255)

               gTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagag                                                                         >  1:676860/1‑93 (MQ=255)

               gTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagag                                                                         <  1:164656/93‑1 (MQ=255)

                 aaaaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagagag                                                                       <  2:2053453/93‑1 (MQ=255)

                  aaaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGc                                                                      <  1:4253272/93‑1 (MQ=255)

                   aaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCATAGGGCATGGTCAACAACGGCCAGAGAGCg                                                                     >  2:3119031/1‑93 (MQ=255)

                    aaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGa                                                                    >  1:4662573/1‑93 (MQ=255)

                    aaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGa                                                                    <  1:3717279/93‑1 (MQ=255)

                    aaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGa                                                                    <  2:330672/93‑1 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                   >  2:452767/1‑93 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                   <  2:4139421/93‑1 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                   <  2:2290737/93‑1 (MQ=255)

                      aTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATc                                                                  <  1:1186926/93‑1 (MQ=255)

                      aTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATc                                                                  >  1:3865588/1‑93 (MQ=255)

                       ttCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCg                                                                 >  2:2200591/1‑93 (MQ=255)

                         ccGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGc                                                               >  2:1619520/1‑93 (MQ=255)

                            gAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGt                                                            >  1:2263019/1‑93 (MQ=255)

                              aCTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTAt                                                          >  2:2154677/1‑93 (MQ=255)

                              aCTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTAt                                                          <  1:1272333/93‑1 (MQ=255)

                                tgtgGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTc                                                        <  2:2112141/93‑1 (MQ=255)

                                tgtgGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTc                                                        >  2:1934162/1‑93 (MQ=255)

                                 gtggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCa                                                       <  1:961901/93‑1 (MQ=255)

                                  tggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCAt                                                      <  2:1057933/93‑1 (MQ=255)

                                  tggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCAt                                                      >  1:2805935/1‑93 (MQ=255)

                                  tggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCAt                                                      >  1:1854026/1‑93 (MQ=255)

                                   ggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCAtt                                                     <  2:1776871/93‑1 (MQ=255)

                                   ggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCAtt                                                     >  1:3029002/1‑93 (MQ=255)

                                      gtATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTaaa                                                  <  1:1355671/93‑1 (MQ=255)

                                       tATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAAc                                                 <  1:472848/93‑1 (MQ=255)

                                       tATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAAc                                                 <  1:1822901/93‑1 (MQ=255)

                                        aTCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAAcc                                                <  2:3214119/93‑1 (MQ=255)

                                          cAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCga                                              <  2:3192874/93‑1 (MQ=255)

                                          cAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCga                                              >  1:28819/1‑93 (MQ=255)

                                          cAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCga                                              <  1:719997/93‑1 (MQ=255)

                                              aCCTGGCTCAATATCTATGAAAGAAAAGCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCgatgat                                          <  1:3357134/93‑1 (MQ=255)

                                               ccTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGata                                         <  2:4157283/93‑1 (MQ=255)

                                               ccTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGata                                         <  2:3111486/93‑1 (MQ=255)

                                               ccTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGata                                         <  2:3514261/93‑1 (MQ=255)

                                                cTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGatat                                        >  1:327817/1‑93 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGAtata                                       >  1:1846790/1‑93 (MQ=255)

                                                   gCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGACCGATCGGCAGTATTCATTAAACCCATGCTATaaa                                     <  2:3903334/93‑1 (MQ=255)

                                                    cTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATaaaa                                    >  2:4134652/1‑92 (MQ=255)

                                                    cTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAAc                                    <  2:3837032/93‑1 (MQ=255)

                                                    cTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAAc                                    >  2:645475/1‑93 (MQ=255)

                                                     tCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACa                                   <  2:2750782/93‑1 (MQ=255)

                                                     tCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACa                                   <  1:3268504/93‑1 (MQ=255)

                                                     tCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACa                                   >  1:352258/1‑93 (MQ=255)

                                                      cAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACAt                                  <  1:1998696/93‑1 (MQ=255)

                                                      cAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACAt                                  <  2:3975440/93‑1 (MQ=255)

                                                      cAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACAt                                  <  2:4276782/93‑1 (MQ=255)

                                                       aaTATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACAtt                                 >  1:3001888/1‑93 (MQ=255)

                                                       aaTATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACAtt                                 <  1:1592555/93‑1 (MQ=255)

                                                        atatCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATtc                                <  2:944383/93‑1 (MQ=255)

                                                         tatCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATtct                               <  2:2854502/93‑1 (MQ=255)

                                                          atCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATtctc                              <  1:4477556/93‑1 (MQ=255)

                                                          atCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATtctc                              <  1:1430322/93‑1 (MQ=255)

                                                           tCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTcc                             >  1:2141687/1‑93 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCg                            <  2:1746351/93‑1 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCg                            <  2:4213029/93‑1 (MQ=255)

                                                             tATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCgg                           >  1:2509239/1‑93 (MQ=255)

                                                             tATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCgg                           <  2:1489219/93‑1 (MQ=255)

                                                             tATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCgg                           >  2:1493281/1‑93 (MQ=255)

                                                             tATGAAAGAAAATCGTCTGCGTACGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCgg                           >  2:2466934/1‑93 (MQ=255)

                                                               tGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGCCAGTATTCATTAAACCGATGTAATAAACATTCTCCGGGa                         >  1:3766939/1‑93 (MQ=255)

                                                                 aaagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACa                       <  2:1226569/93‑1 (MQ=255)

                                                                    gaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAg                    >  1:988135/1‑93 (MQ=255)

                                                                    gaaaATCGTCTGCGTAGGGCATGGTCAACAAAGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAg                    >  2:3755477/1‑93 (MQ=255)

                                                                      aaaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGaa                  <  1:1568531/93‑1 (MQ=255)

                                                                       aaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAAt                 <  1:1793878/93‑1 (MQ=255)

                                                                       aaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAAt                 >  2:4045759/1‑93 (MQ=255)

                                                                       aaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAAt                 <  2:2785030/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg               <  2:1812786/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg               <  2:597438/93‑1 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc            <  1:3133333/93‑1 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc            <  2:2094351/93‑1 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc            >  2:4507990/1‑93 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc            <  2:2500044/93‑1 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc            >  1:2940150/1‑93 (MQ=255)

                                                                              tGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttctt          <  1:2122518/93‑1 (MQ=255)

                                                                               gCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttcttc         <  1:605181/93‑1 (MQ=255)

                                                                                cGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCa        <  1:2648823/93‑1 (MQ=255)

                                                                                    gggCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGt    <  2:2362473/93‑1 (MQ=255)

                                                                                      gCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGttt  >  2:3354139/1‑93 (MQ=255)

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTT  >  NC_000913/4575788‑4575966

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