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%	551.4 / ‑4.7	170	Y344H (TAT→CAT)	mcrB	McrB
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (58/110); total (58/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 175 total aligned reads displayed.

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTT  >  NC_000913/4575788‑4575967

cTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGt                                                                                         <  2:1907881/93‑1 (MQ=255)

cTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGt                                                                                         <  2:1596660/93‑1 (MQ=255)

cTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGt                                                                                         >  1:2692223/1‑93 (MQ=255)

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

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

 tGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTc                                                                                        <  1:2123794/93‑1 (MQ=255)

  gCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTca                                                                                       <  1:2543584/93‑1 (MQ=255)

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

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

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

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

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

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

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

       tttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCaacaac                                                                                  >  2:842729/1‑93 (MQ=255)

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

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

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

            tCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCa                                                                             >  1:3114826/1‑93 (MQ=255)

             cAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCag                                                                            >  1:652002/1‑93 (MQ=255)

             cAGAAAAAAATTCCGGAACTGCGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAATGGCCag                                                                            >  1:652004/1‑93 (MQ=255)

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

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                           >  1:4043393/1‑93 (MQ=255)

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

                tAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagaga                                                                         <  1:1357931/93‑1 (MQ=255)

                tAAAAAATGCCGGAACTGTGGTGTATCAAATCCTGGCTCAATATCTATGAAAGAAAATCGTCTTCGTAGGGCATGGTCAACAACGGCCagaga                                                                         <  1:1357929/93‑1 (MQ=255)

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

                  aaaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGc                                                                       >  1:2925439/1‑93 (MQ=255)

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

                   aaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCg                                                                      <  2:1760444/93‑1 (MQ=255)

                   aaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCg                                                                      <  1:3332008/93‑1 (MQ=255)

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

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

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

                       ttCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCg                                                                  <  2:3784907/93‑1 (MQ=255)

                        tCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCgg                                                                 >  2:1533022/1‑93 (MQ=255)

                          cGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCa                                                               >  1:2144563/1‑93 (MQ=255)

                           ggAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAg                                                              >  1:4161597/1‑93 (MQ=255)

                           ggAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAg                                                              >  2:3165064/1‑93 (MQ=255)

                           ggAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGCCAg                                                              >  1:1886416/1‑93 (MQ=255)

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

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

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

                             aaCTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTa                                                            <  2:3248200/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                    gtgtATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTa                                                     >  1:4648740/1‑93 (MQ=255)

                                     tgtATCAAAACCTGGCTCAATATCTATGAAATAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCACTaa                                                    <  1:3715801/93‑1 (MQ=255)

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

                                        accAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAAcc                                                 <  2:328983/91‑1 (MQ=255)

                                           aaaaCCTGGCTCAATATCTATGAAAGAACATCGTCTGCGTAGGGCCTGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCgag                                              >  2:1362793/1‑92 (MQ=255)

                                           aaaaCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCgat                                              >  1:3128046/1‑93 (MQ=255)

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

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

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

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

                                                  ggCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATaa                                       >  2:309634/1‑93 (MQ=255)

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

                                                     tCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCAGGGTCAACAAGGGCCAGAGAGCGATCGGCAGTATTCTTTAAACCGATGATAAAAAct                                    >  2:13949/1‑92 (MQ=255)

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

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

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

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

                                                           tCTGTGAAAGAAACTCGTCTGCGTAGGGCATTGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTcc                              <  1:984680/93‑1 (MQ=255)

                                                           tCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTcc                              <  1:984681/93‑1 (MQ=255)

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

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

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCg                             >  1:3067009/1‑93 (MQ=255)

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

                                                             tATGAAAGAAAAGCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGCATTCATTAAACCGATGATATAAACATGCTCCgg                            <  2:1902715/93‑1 (MQ=255)

                                                               tGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGa                          <  1:4392485/93‑1 (MQ=255)

                                                                gaaagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGAc                         <  1:1710120/93‑1 (MQ=255)

                                                                  aagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACAt                       <  2:2471264/93‑1 (MQ=255)

                                                                  aagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACAt                       >  2:277376/1‑93 (MQ=255)

                                                                   agaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATa                      <  1:248506/93‑1 (MQ=255)

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

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

                                                                      aaaTCGTCTGCGTAGGGCATGGTAAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGaa                   <  2:2603881/93‑1 (MQ=255)

                                                                          cGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGt               >  1:3517860/1‑93 (MQ=255)

                                                                           gTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGtt              <  1:1300088/93‑1 (MQ=255)

                                                                           gTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGtt              >  2:2079525/1‑93 (MQ=255)

                                                                           gTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGtt              >  1:3263652/1‑93 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTTTACGGGACATATAATCGttc             >  2:823260/1‑93 (MQ=255)

                                                                             cTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttct            <  1:2749900/93‑1 (MQ=255)

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

                                                                                 gTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCAt        <  2:3926826/93‑1 (MQ=255)

                                                                                 gTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCAt        >  2:1326490/1‑93 (MQ=255)

                                                                                   aGGGTATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATct      >  1:1258432/1‑92 (MQ=255)

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

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

                                                                                     ggCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtt    <  2:3662570/93‑1 (MQ=255)

                                                                                     ggCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtt    <  2:2941680/93‑1 (MQ=255)

                                                                                      gCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGttt   <  1:937361/93‑1 (MQ=255)

                                                                                      gCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGttt   <  2:1902628/93‑1 (MQ=255)

                                                                                      gCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGttt   <  2:196079/93‑1 (MQ=255)

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTT  >  NC_000913/4575788‑4575967

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