BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,519	T→C	T122A (ACG→GCG)	yddE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,519	0	T	C	98.7%	473.3 / ‑2.6	151	T122A (ACG→GCG)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (67/82); total (68/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 156 total aligned reads displayed.

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTG  >  NC_000913/1533429‑1533607

aGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAc                                                                                        >  2:671647/1‑93 (MQ=255)

 gCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAcc                                                                                       >  2:1532849/1‑93 (MQ=255)

 gCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCAcc                                                                                       <  1:2009390/93‑1 (MQ=255)

    tGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTg                                                                                    >  2:2983780/1‑93 (MQ=255)

     ggCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGAGCACCTTGt                                                                                   <  1:3394934/93‑1 (MQ=255)

       cAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTc                                                                                 <  1:4296402/93‑1 (MQ=255)

       cAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTc                                                                                 >  1:2557157/1‑93 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGc                                                                             <  1:3793228/93‑1 (MQ=255)

           atatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGc                                                                             <  2:441938/93‑1 (MQ=255)

            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                            >  1:1821181/1‑93 (MQ=255)

            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCg                                                                            >  2:625207/1‑93 (MQ=255)

            tatCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGAGCACCTTGTTCCAGCg                                                                            <  2:746227/93‑1 (MQ=255)

              tCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaa                                                                          >  2:11528/1‑93 (MQ=255)

               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                         <  1:3964997/93‑1 (MQ=255)

               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                         <  2:838593/93‑1 (MQ=255)

               cGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGaaa                                                                         <  1:4210071/93‑1 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                       >  2:3875357/1‑93 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                       <  1:3071042/93‑1 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                       >  1:3302375/1‑93 (MQ=255)

                 tCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATa                                                                       <  1:4746946/93‑1 (MQ=255)

                  ccTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATAc                                                                      >  1:3414422/1‑93 (MQ=255)

                   cTCAGTGAGATGAAGCGCGTTGATAATCGCCGAACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACg                                                                     >  1:368097/1‑93 (MQ=255)

                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                    >  1:4346288/1‑93 (MQ=255)

                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                    >  2:4036023/1‑93 (MQ=255)

                    tCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGa                                                                    <  2:408433/93‑1 (MQ=255)

                     cAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGAt                                                                   <  1:283751/93‑1 (MQ=255)

                       gTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATaa                                                                 >  2:1678103/1‑93 (MQ=255)

                         gagaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatc                                                               >  1:942956/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatca                                                              >  1:4635592/1‑93 (MQ=255)

                          agaTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTCCCAGCGAAATACGATAatca                                                              <  1:323668/93‑1 (MQ=255)

                            aTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCATCGAAATACGATAatcatc                                                            >  2:852920/1‑93 (MQ=255)

                            aTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcatc                                                            <  2:4747658/93‑1 (MQ=255)

                            aTGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAatcatc                                                            >  2:3109410/1‑93 (MQ=255)

                              gAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGt                                                          <  1:2670579/93‑1 (MQ=255)

                                aGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCTGCGAAATACGATAATCATCGTtg                                                        >  2:3431582/1‑93 (MQ=255)

                                 gcgcGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgt                                                       <  1:4128979/93‑1 (MQ=255)

                                 gcgcGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTtgt                                                       <  2:3566340/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGGACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                     <  1:2042607/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                     <  1:4684599/93‑1 (MQ=255)

                                   gcgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGc                                                     <  2:1457940/93‑1 (MQ=255)

                                    cgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCt                                                    <  2:263792/93‑1 (MQ=255)

                                    cgTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCt                                                    <  2:263791/93‑1 (MQ=255)

                                     gTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCtt                                                   <  1:31009/93‑1 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                  >  1:2004921/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGtttt                                                  >  2:2432764/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                  >  1:3563386/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                  >  1:3568717/1‑93 (MQ=255)

                                      ttGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCttt                                                  <  1:71071/93‑1 (MQ=255)

                                        gATAATCGCCGCACGTGTTTCGGCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTc                                                >  1:845803/1‑93 (MQ=255)

                                          tAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGa                                              >  2:1165157/1‑93 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                             <  1:2307978/93‑1 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                             <  2:3659532/93‑1 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                             >  1:2058926/1‑93 (MQ=255)

                                           aaTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGAt                                             >  2:1902209/1‑93 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                            <  1:4697472/93‑1 (MQ=255)

                                            aTCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATa                                            >  1:4001664/1‑93 (MQ=255)

                                             tCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAg                                           <  1:2924514/93‑1 (MQ=255)

                                              cgccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGt                                          <  1:1604357/93‑1 (MQ=255)

                                               gccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCCCCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                         <  2:2757555/93‑1 (MQ=255)

                                               gccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                         >  1:1266739/1‑93 (MQ=255)

                                               gccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                         >  2:4166307/1‑93 (MQ=255)

                                               gacgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                         <  1:1326286/91‑1 (MQ=255)

                                               cccgcACGTTTTTCACCTTACAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTc                                         <  2:1231652/92‑1 (MQ=255)

                                                ccgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCa                                        <  1:4455353/93‑1 (MQ=255)

                                                 cgcACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAc                                       >  2:1575001/1‑93 (MQ=255)

                                                   cACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCCCTTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgc                                     <  1:1193819/93‑1 (MQ=255)

                                                    aCGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCAcgcg                                    >  2:305940/1‑93 (MQ=255)

                                                     cGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGa                                   >  1:902393/1‑93 (MQ=255)

                                                      gtgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCTAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAt                                  <  1:3790256/93‑1 (MQ=255)

                                                      gtgtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGAt                                  >  1:3246456/1‑93 (MQ=255)

                                                        gtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGt                                >  2:2974467/1‑93 (MQ=255)

                                                        gtTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGt                                <  1:1622343/93‑1 (MQ=255)

                                                         tttCACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGtt                               >  1:1012683/1‑93 (MQ=255)

                                                            cACCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTc                            <  2:2283116/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           <  1:464359/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           <  2:2718518/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           <  2:1560157/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           <  1:3756665/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           <  1:2064212/93‑1 (MQ=255)

                                                             aCCTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTcc                           >  1:2124741/1‑93 (MQ=255)

                                                               cTTCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTg                         >  2:3962653/1‑93 (MQ=255)

                                                                 tCCAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGcc                       >  2:234013/1‑93 (MQ=255)

                                                                  ccAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTCCTTTTCGATAGTCACGCGATGTTTTCCTGCCa                      <  1:2321199/93‑1 (MQ=255)

                                                                   cAGCGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAg                     <  2:1707296/93‑1 (MQ=255)

                                                                      cGGTGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGa                  <  2:4517593/93‑1 (MQ=255)

                                                                         tGGCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGt               <  1:2858303/93‑1 (MQ=255)

                                                                          ggCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGtt              >  1:3491488/1‑93 (MQ=255)

                                                                          ggCTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGtt              <  2:4403745/93‑1 (MQ=255)

                                                                            cTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTg            >  2:4720656/1‑93 (MQ=255)

                                                                            cTCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTg            <  1:1367300/93‑1 (MQ=255)

                                                                             tCAAAGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGc           <  2:2916174/93‑1 (MQ=255)

                                                                               aaaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCa         <  1:1253086/93‑1 (MQ=255)

                                                                                aaGCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAg        >  1:3341804/1‑93 (MQ=255)

                                                                                  gCCCGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGAt      <  2:3434580/93‑1 (MQ=255)

                                                                                   cccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTACTGCCAGCGATGTTTGCCAGATc     >  1:2743398/1‑93 (MQ=255)

                                                                                    ccGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCg    <  2:2484285/93‑1 (MQ=255)

                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCTTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt   <  2:4026689/93‑1 (MQ=255)

                                                                                     cGGCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGt   <  1:1683323/93‑1 (MQ=255)

                                                                                      ggCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg  <  1:1836362/93‑1 (MQ=255)

                                                                                      ggCGCACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTg  >  2:379741/1‑93 (MQ=255)

AGCCTGGCAGAATATCGTCCTCAGTGAGATGAAGCGCGTTGATAATCGCCGCACGTGTTTCACCTTCCAGCGGTGGCTCAAAGCCCGGCGTACCTTGTTCCAGCGAAATACGATAATCATCGTTGTGCTTTTCGATAGTCACGCGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTG  >  NC_000913/1533429‑1533607

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