BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,064,328	A→G	N357S (AAT→AGT)	yccE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,328	0	A	G	84.9%	562.9 / 90.4	238	N357S (AAT→AGT)	yccE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (22/14); new base G (118/84); total (140/98)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 248 total aligned reads displayed.

ATACCGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTATGACG  >  NC_000913/1064238‑1064420

aTACCGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGtg                                                                                            >  1:12812/1‑93 (MQ=255)

 tgccGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAgtgt                                                                                           <  1:5198171/91‑1 (MQ=255)

 tACCGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAgtgt                                                                                           <  1:7439605/93‑1 (MQ=255)

 tACCGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAgtgt                                                                                           <  1:2435617/93‑1 (MQ=255)

      acacGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCAt                                                                                      >  1:1042953/1‑93 (MQ=255)

      acacGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCAt                                                                                      >  2:2902916/1‑93 (MQ=255)

        acGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCa                                                                                    >  2:6547698/1‑93 (MQ=255)

          gTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAAt                                                                                  <  2:2320941/93‑1 (MQ=255)

          gTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAAt                                                                                  <  1:329944/93‑1 (MQ=255)

           tCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAAtc                                                                                 >  2:3898065/1‑93 (MQ=255)

             ggATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAAtctc                                                                               >  2:4632449/1‑93 (MQ=255)

             ggATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAAtctc                                                                               >  2:6355585/1‑93 (MQ=255)

             ggATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAAtctc                                                                               >  2:535979/1‑93 (MQ=255)

                 aGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTg                                                                           <  2:1968930/93‑1 (MQ=255)

                  ggATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTgg                                                                          >  1:5822715/1‑93 (MQ=255)

                  ggATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTgg                                                                          >  2:790466/1‑93 (MQ=255)

                   gATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTggg                                                                         <  1:4750508/93‑1 (MQ=255)

                   gATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTggg                                                                         >  1:2165257/1‑93 (MQ=255)

                       aaaGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGgaat                                                                     >  1:1571639/1‑93 (MQ=255)

                         aGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGgaatga                                                                   >  1:544686/1‑93 (MQ=255)

                          gCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGgaatgaa                                                                  >  2:3632469/1‑93 (MQ=255)

                           cGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGgaatgaat                                                                 <  2:4417067/93‑1 (MQ=255)

                              tGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCt                                                              >  1:2709929/1‑93 (MQ=255)

                               gaagaaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCtt                                                             >  1:4322699/1‑93 (MQ=255)

                                aagaaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCttt                                                            <  2:1850932/93‑1 (MQ=255)

                                 agaaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCtttt                                                           >  2:2684794/1‑93 (MQ=255)

                                 agaaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCtttt                                                           >  2:2116291/1‑93 (MQ=255)

                                 agaaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCtttt                                                           >  1:3359728/1‑93 (MQ=255)

                                   aaTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTTTTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCa                                                         >  1:2319151/1‑93 (MQ=255)

                                    aTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCaa                                                        <  1:6851946/93‑1 (MQ=255)

                                    aTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCaa                                                        >  1:1160342/1‑93 (MQ=255)

                                    aTTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCaa                                                        <  2:6491228/93‑1 (MQ=255)

                                      tttCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACa                                                      <  1:7093411/93‑1 (MQ=255)

                                        tCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGa                                                    >  2:6988130/1‑93 (MQ=255)

                                         cTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGAt                                                   >  2:5244466/1‑93 (MQ=255)

                                          ttgttgTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATg                                                  >  2:3146877/1‑93 (MQ=255)

                                           tgttgtTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGa                                                 >  2:4164866/1‑93 (MQ=255)

                                            gttgttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGaa                                                <  2:2791329/93‑1 (MQ=255)

                                             ttgttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAg                                               >  1:964915/1‑93 (MQ=255)

                                              tgttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGa                                              <  2:3405750/93‑1 (MQ=255)

                                               gttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGAt                                             <  1:5699593/93‑1 (MQ=255)

                                               gttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGAt                                             >  2:193962/1‑93 (MQ=255)

                                                ttCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATc                                            >  1:7019038/1‑93 (MQ=255)

                                                 tCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCt                                           <  2:6707921/93‑1 (MQ=255)

                                                 tCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCt                                           <  2:846454/93‑1 (MQ=255)

                                                 tCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCt                                           <  1:6390894/93‑1 (MQ=255)

                                                  cTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCtt                                          >  2:5551768/1‑93 (MQ=255)

                                                  cTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCtt                                          <  2:3972999/93‑1 (MQ=255)

                                                  cTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCtt                                          >  2:1563513/1‑93 (MQ=255)

                                                  cTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCtt                                          <  2:625070/93‑1 (MQ=255)

                                                   tttCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTg                                         >  1:3495157/1‑93 (MQ=255)

                                                   tttCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTg                                         <  2:4598708/93‑1 (MQ=255)

                                                   tttCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTg                                         >  1:7496617/1‑93 (MQ=255)

                                                   tttCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTg                                         >  1:2526676/1‑93 (MQ=255)

                                                     tCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAg                                       <  1:6818949/93‑1 (MQ=255)

                                                       aGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTa                                     >  2:6033474/1‑93 (MQ=255)

                                                       aGATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTa                                     <  2:3274066/93‑1 (MQ=255)

                                                        gATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTAt                                    >  2:6618835/1‑93 (MQ=255)

                                                        gATAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTAt                                    <  1:5634551/93‑1 (MQ=255)

                                                         aTAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATg                                   >  2:4082520/1‑93 (MQ=255)

                                                         aTAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATg                                   >  2:7339834/1‑93 (MQ=255)

                                                         aTAGCATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATg                                   >  2:2415641/1‑93 (MQ=255)

                                                          tAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATgg                                  >  2:1680584/1‑93 (MQ=255)

                                                             cATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTa                               >  2:2629929/1‑93 (MQ=255)

                                                             cATTTGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTa                               >  2:585940/1‑93 (MQ=255)

                                                               tttGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTAtt                             >  1:4906841/1‑93 (MQ=255)

                                                               tttGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTAtt                             >  1:607114/1‑93 (MQ=255)

                                                                 tGCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTa                           <  1:2002368/93‑1 (MQ=255)

                                                                  gCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTaa                          <  2:2727054/93‑1 (MQ=255)

                                                                  gCTCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTaa                          >  1:4782927/1‑93 (MQ=255)

                                                                  gCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTaa                          <  1:4549677/93‑1 (MQ=255)

                                                                  gCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTaa                          <  1:6987339/93‑1 (MQ=255)

                                                                   cTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAAt                         >  1:2745116/1‑93 (MQ=255)

                                                                    tCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATg                        >  2:7184982/1‑93 (MQ=255)

                                                                    tCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATg                        <  2:577343/93‑1 (MQ=255)

                                                                    tCATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATg                        >  2:1030400/1‑93 (MQ=255)

                                                                     cATTATTTATTGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGa                       >  1:6736868/1‑93 (MQ=255)

                                                                       ttatttatTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAac                     >  1:1319641/1‑93 (MQ=255)

                                                                        tatttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAaca                    >  2:4387769/1‑93 (MQ=255)

                                                                        tatttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAaca                    >  2:6097395/1‑93 (MQ=255)

                                                                        tatttattGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAaca                    >  1:886890/1‑93 (MQ=255)

                                                                         atttattGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAacac                   >  1:3663373/1‑93 (MQ=255)

                                                                          tttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACt                  >  1:3390766/1‑93 (MQ=255)

                                                                          tttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACt                  <  1:272019/93‑1 (MQ=255)

                                                                           ttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTg                 <  2:3541572/93‑1 (MQ=255)

                                                                           ttattGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTg                 <  1:2892398/93‑1 (MQ=255)

                                                                            tattGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTgg                <  1:1345867/93‑1 (MQ=255)

                                                                               tGTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCaa             <  1:3613837/93‑1 (MQ=255)

                                                                                gTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAAt            <  1:4652534/93‑1 (MQ=255)

                                                                                gTTAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAAt            <  2:7210691/93‑1 (MQ=255)

                                                                                 ttAATAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATc           >  1:4023626/1‑93 (MQ=255)

                                                                                   aaTAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCAt         <  1:3156570/93‑1 (MQ=255)

                                                                                   aaTAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCAt         >  1:4074410/1‑93 (MQ=255)

                                                                                   aaTAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCAt         >  2:2029078/1‑93 (MQ=255)

                                                                                    aTAGTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCAtt        <  1:7139255/93‑1 (MQ=255)

                                                                                    aTAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCAtt        >  1:347043/1‑93 (MQ=255)

                                                                                      agtagtGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTAt      <  1:6538449/93‑1 (MQ=255)

                                                                                      agtagtGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTAt      <  2:5712551/93‑1 (MQ=255)

                                                                                      aGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTAt      <  1:6011863/93‑1 (MQ=255)

                                                                                       gTAGTGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTATg     >  1:6762423/1‑93 (MQ=255)

ATACCGACACGTCGGATAGGATTAAAGCGGTGAAGAATTTTCTTGTTGTTCTTTCAGATAGCATTTGCTCATTATTTATTGTTAATAGTAATGTGGCATCAATCTCCTTGGGGAATGAATCCTTTTCAACAGATGAAGATCTTGAGTATGGTTATTTAATGAACACTGGCAATCATTATGACG  >  NC_000913/1064238‑1064420

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