BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,973,216	C→T	A153A (GCC→GCT)	wecE →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,973,216	0	C	T	98.3%	606.6 / ‑2.6	175	A153A (GCC→GCT)	wecE	putative regulator
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (81/91); total (83/92)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 186 total aligned reads displayed.

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

cGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGAtgctgc                                                                                         <  2:736698/93‑1 (MQ=255)

cGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGAtgctgc                                                                                         <  1:1645159/93‑1 (MQ=255)

 gTGCGGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATgctgct                                                                                        <  1:965316/93‑1 (MQ=255)

  tGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTc                                                                                       <  1:1523140/93‑1 (MQ=255)

  tGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTc                                                                                       >  1:2591910/1‑93 (MQ=255)

    ccGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAg                                                                                     >  1:2100583/1‑93 (MQ=255)

    ccGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAg                                                                                     <  2:2189285/93‑1 (MQ=255)

    ccGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAg                                                                                     <  1:2503775/93‑1 (MQ=255)

     ctgTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAgg                                                                                    <  1:1918049/91‑1 (MQ=255)

     cGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAgg                                                                                    <  2:3061926/93‑1 (MQ=255)

       gTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGc                                                                                  >  2:2119910/1‑93 (MQ=255)

          cATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTg                                                                               >  2:2116518/1‑93 (MQ=255)

          cATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTg                                                                               <  1:1253014/93‑1 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGCAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGCAGATGCTGCTCCGCGCGTGa                                                                              >  1:2219831/1‑93 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGa                                                                              >  1:2969763/1‑93 (MQ=255)

             tACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATg                                                                            <  1:3318959/93‑1 (MQ=255)

              aCGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGt                                                                           >  1:1012949/1‑93 (MQ=255)

                 cgGGTGTGGCCTGCGAAATGGACACCATTATTGCGTTGGCGAAAAAGCATAATTTGATTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTACa                                                                        >  2:1728520/1‑93 (MQ=255)

                  gggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCAc                                                                       >  2:54309/1‑93 (MQ=255)

                   ggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACt                                                                      >  1:3887625/1‑93 (MQ=255)

                   ggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACt                                                                      <  1:566614/93‑1 (MQ=255)

                    gtgtGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGGGATGTCCACtt                                                                     >  2:4034675/1‑93 (MQ=255)

                     tgtgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTa                                                                    >  1:1510177/1‑93 (MQ=255)

                       tgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACa                                                                  <  2:1308951/93‑1 (MQ=255)

                       tgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACa                                                                  >  1:2951509/1‑93 (MQ=255)

                        ggCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaa                                                                 <  1:2397017/93‑1 (MQ=255)

                         gCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaaa                                                                <  2:1006733/93‑1 (MQ=255)

                         gCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaaa                                                                <  2:2151348/93‑1 (MQ=255)

                          ccTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAg                                                               <  2:347682/93‑1 (MQ=255)

                           cTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTATAAGATGCTGCTCAGGGCGGGATGTCCACTTACAAAgg                                                              >  2:3743833/1‑93 (MQ=255)

                            tGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             <  1:4632610/93‑1 (MQ=255)

                            tGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             >  1:2049947/1‑93 (MQ=255)

                            tGCGAAATCGACTCTATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGTTAGATGATTCTGCTCAGGGCGTGATGTCCTCTTACAAAggc                                                             <  2:496600/93‑2 (MQ=255)

                             gCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGc                                                            >  2:3377564/1‑93 (MQ=255)

                              cGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCg                                                           <  1:476062/93‑1 (MQ=255)

                               gAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGt                                                          <  1:4679041/93‑1 (MQ=255)

                                aaaTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTg                                                         >  2:2776831/1‑93 (MQ=255)

                                aaaTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGCTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTg                                                         >  2:4421254/1‑93 (MQ=255)

                                  aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCa                                                       >  1:1187798/1‑93 (MQ=255)

                                  aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGAGTGCa                                                       >  2:2011143/1‑93 (MQ=255)

                                   tGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCAc                                                      >  2:2973162/1‑93 (MQ=255)

                                     gACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTg                                                    <  1:4025110/93‑1 (MQ=255)

                                      acacCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTgg                                                   <  1:340162/93‑1 (MQ=255)

                                       cacCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTggg                                                  >  1:2997968/1‑93 (MQ=255)

                                        acCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGa                                                 <  2:4119512/93‑1 (MQ=255)

                                        acCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGa                                                 >  1:1842141/1‑93 (MQ=255)

                                        acCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGa                                                 <  2:1599583/93‑1 (MQ=255)

                                         ccATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGac                                                >  1:2872823/1‑92 (MQ=255)

                                          cATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAc                                               >  2:1086061/1‑93 (MQ=255)

                                           aTTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAcc                                              >  2:4609015/1‑93 (MQ=255)

                                           aTTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGCCGTGATGTCCACTTACACAGGGCGTGCACTGGGAAcc                                              >  2:4700151/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             <  2:440668/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  1:3026441/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  2:127741/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  1:1246608/93‑1 (MQ=255)

                                              aTGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAtt                                           >  2:927801/1‑93 (MQ=255)

                                              aTGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAtt                                           <  1:2491917/93‑1 (MQ=255)

                                               tGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTg                                          >  2:779453/1‑93 (MQ=255)

                                                ggCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGTGAGTGCACTGGGAACCATTgg                                         >  1:2723224/1‑93 (MQ=255)

                                                ggCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTgg                                         <  2:3877747/93‑1 (MQ=255)

                                                 gCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGt                                        >  2:1526600/1‑93 (MQ=255)

                                                 gCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGt                                        <  1:3683413/93‑1 (MQ=255)

                                                  cGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTc                                       <  1:3242718/93‑1 (MQ=255)

                                                  cGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTc                                       <  2:2442463/93‑1 (MQ=255)

                                                   gTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCa                                      >  1:1529281/1‑93 (MQ=255)

                                                   gTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCa                                      >  1:2658584/1‑93 (MQ=255)

                                                        cGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTg                                 <  2:2294633/93‑1 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTgg                                >  1:1064612/1‑93 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTgg                                >  2:2569914/1‑93 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTgg                                >  1:695832/1‑93 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGGGGTAGAAGATGCTGCTCGGGGCGTGATGGCCACTTACAAAGGGCGTGCACTGGGAAGCATTGGTCATATTgg                                >  2:2357052/1‑93 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATAttcgg                               >  2:3105578/1‑90 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               <  2:312470/93‑1 (MQ=255)

                                                           aaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgct                              <  1:4263263/93‑1 (MQ=255)

                                                           aaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgct                              <  2:4712343/93‑1 (MQ=255)

                                                            aaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctg                             >  2:2635102/1‑93 (MQ=255)

                                                              aGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgct                           <  2:1816566/93‑1 (MQ=255)

                                                              aGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgct                           <  2:2468141/93‑1 (MQ=255)

                                                              aGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgct                           >  1:1977376/1‑93 (MQ=255)

                                                              aGCATAATTTGTGTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTAAAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgct                           <  2:2801071/93‑1 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          <  1:2630377/93‑1 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          <  2:358101/93‑1 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        >  2:2046413/1‑93 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        <  2:313725/93‑1 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        >  1:3831961/1‑93 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        <  1:3983986/93‑1 (MQ=255)

                                                                  tAATTTGTTTGTGGTAGAAGAGGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAg                       <  1:4413131/93‑1 (MQ=255)

                                                                   aaTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                      <  1:124132/93‑1 (MQ=255)

                                                                    aTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                     <  2:3156975/93‑1 (MQ=255)

                                                                       tgtttgtGGTAGAAGATGGTGCTCAGGGCGTGATGTCCACTTTCAAAGGGCGTGCACTGGGAACCATTGGTATTATTGGCTGCTTTAGCTTcc                  <  2:3696616/93‑1 (MQ=255)

                                                                       tgtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTcc                  >  1:4559804/1‑93 (MQ=255)

                                                                        gtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCa                 >  1:331116/1‑93 (MQ=255)

                                                                         tttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCAt                <  1:1149418/93‑1 (MQ=255)

                                                                            gtgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaa             >  2:4269284/1‑93 (MQ=255)

                                                                              ggTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAc           >  1:2741681/1‑93 (MQ=255)

                                                                                 agaagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaa        <  1:4741496/93‑1 (MQ=255)

                                                                                   aagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaa      >  1:4493581/1‑93 (MQ=255)

                                                                                    agaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaaa     <  2:1756734/93‑1 (MQ=255)

                                                                                     gaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    <  2:173193/93‑1 (MQ=255)

                                                                                     gaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    <  1:854228/93‑1 (MQ=255)

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

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