BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,596,295	A→G	V332A (GTG→GCG)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,596,295	0	A	G	99.7%	932.1 / ‑5.0	287	V332A (GTG→GCG)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (140/146); total (141/146)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 298 total aligned reads displayed.

GTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAGGTG  >  NC_000913/3596205‑3596384

 tgccgtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCgg                                                                                        <  1:5478316/93‑1 (MQ=255)

   ccgtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAg                                                                                      >  1:3070460/1‑93 (MQ=255)

   ccgtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAg                                                                                      <  2:4401304/93‑1 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    <  1:4572044/93‑1 (MQ=255)

      ttGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTc                                                                                   >  2:3643088/1‑93 (MQ=255)

       tGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCg                                                                                  >  2:355582/1‑93 (MQ=255)

        ggCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGc                                                                                 >  2:2950907/1‑93 (MQ=255)

        ggCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGc                                                                                 <  2:4310708/93‑1 (MQ=255)

        ggCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGc                                                                                 <  2:1742911/93‑1 (MQ=255)

          cGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCtt                                                                               <  1:85000/93‑1 (MQ=255)

           gTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCttt                                                                              >  2:6108541/1‑93 (MQ=255)

              ccAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAg                                                                           <  2:2618270/93‑1 (MQ=255)

              ccAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAg                                                                           <  2:3284436/93‑1 (MQ=255)

               cAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAgg                                                                          <  1:160092/93‑1 (MQ=255)

                aGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGt                                                                         <  2:917600/93‑1 (MQ=255)

                 gTCAAATACGCCGAACTCAAAGCCTTGCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTa                                                                        <  2:3360679/93‑1 (MQ=255)

                   cAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAtt                                                                      <  2:4242478/93‑1 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     <  1:5363776/93‑1 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     >  1:5285778/1‑93 (MQ=255)

                     aaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTg                                                                    <  2:843947/93‑1 (MQ=255)

                     aaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTg                                                                    <  1:1189416/93‑1 (MQ=255)

                      aTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTgg                                                                   >  1:1046638/1‑93 (MQ=255)

                           ccGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGAtt                                                              <  2:5404267/93‑1 (MQ=255)

                              aaCTCAAATCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCa                                                           >  1:6566931/1‑93 (MQ=255)

                                 tCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGcc                                                        >  2:1586080/1‑93 (MQ=255)

                                 ctaaaaCCTTTCAAATCCCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGcc                                                        <  2:2832492/91‑1 (MQ=255)

                                  cAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCg                                                       >  2:401850/1‑93 (MQ=255)

                                  cAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCg                                                       >  2:6296358/1‑93 (MQ=255)

                                   aaaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCgg                                                      >  1:1063048/1‑93 (MQ=255)

                                    aaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGa                                                     >  2:6556402/1‑93 (MQ=255)

                                    aaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGa                                                     <  2:6691444/93‑1 (MQ=255)

                                      gCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGCCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATc                                                   <  2:4874852/93‑1 (MQ=255)

                                         tttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTc                                                >  2:6177880/1‑93 (MQ=255)

                                         tttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTc                                                <  1:4789080/93‑1 (MQ=255)

                                         tttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTc                                                >  2:526575/1‑93 (MQ=255)

                                         tttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTc                                                >  1:2773420/1‑93 (MQ=255)

                                             aGATCGTCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGAc                                            >  2:5635790/1‑93 (MQ=255)

                                             aGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGAc                                            <  1:622809/93‑1 (MQ=255)

                                              gATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGGCGGATCGTCAGACt                                           >  1:4856351/1‑93 (MQ=255)

                                               aTCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTg                                          <  1:4762528/93‑1 (MQ=255)

                                                tCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGa                                         <  2:4469485/93‑1 (MQ=255)

                                                 cGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGTAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCTGACTGAt                                        <  1:6360520/93‑1 (MQ=255)

                                                   ccTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTc                                      >  1:3448246/1‑93 (MQ=255)

                                                    cTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCa                                     <  2:6590191/93‑1 (MQ=255)

                                                    cTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCa                                     <  1:6604729/93‑1 (MQ=255)

                                                      ttCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAgg                                   <  2:126876/93‑1 (MQ=255)

                                                        ctcATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGCGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGcc                                 >  1:2131248/1‑93 (MQ=255)

                                                        ctcATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGcc                                 <  1:1311836/93‑1 (MQ=255)

                                                        ctcATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGcc                                 <  1:757849/93‑1 (MQ=255)

                                                          cATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCg                               >  1:1721888/1‑93 (MQ=255)

                                                           aTCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGc                              >  1:2508080/1‑93 (MQ=255)

                                                           aTCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGc                              >  1:1930589/1‑93 (MQ=255)

                                                           aTCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGc                              >  1:5108823/1‑93 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             <  1:5889696/93‑1 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             <  1:601319/93‑1 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             <  2:2449019/93‑1 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             <  2:4924976/93‑1 (MQ=255)

                                                             cccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCt                            <  1:3121775/93‑1 (MQ=255)

                                                             cccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCt                            <  2:6454687/93‑1 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  2:4849122/1‑93 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  2:2739206/1‑93 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  1:3665452/1‑93 (MQ=255)

                                                               cAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGc                          >  2:220366/1‑93 (MQ=255)

                                                               cAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGc                          >  1:2021787/1‑93 (MQ=255)

                                                               cAGGTCAGCGGTCCCATTACGGTATCCGCGGAGGTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGc                          >  2:1992/1‑93 (MQ=255)

                                                                   tCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAg                      >  1:1080406/1‑93 (MQ=255)

                                                                   tCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAg                      >  2:3805106/1‑93 (MQ=255)

                                                                     aGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAt                    <  2:4731792/93‑1 (MQ=255)

                                                                      gCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   >  1:5728283/1‑93 (MQ=255)

                                                                      gCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   >  2:1401371/1‑93 (MQ=255)

                                                                      gCGGTCCCATGACGGTATTCGCGGAGTTCGCTTTCAGGTATTTGACGATTTCAGCCTGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   <  1:2874840/93‑1 (MQ=255)

                                                                       cGGTCCCATTCCGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTg                  <  1:3863575/93‑1 (MQ=255)

                                                                       cGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTg                  <  2:544585/93‑1 (MQ=255)

                                                                         gTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCa                <  2:6333216/93‑1 (MQ=255)

                                                                         gTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCa                >  1:6688966/1‑93 (MQ=255)

                                                                         gTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCa                >  2:1136709/1‑93 (MQ=255)

                                                                          tCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAg               >  2:4128294/1‑93 (MQ=255)

                                                                          tCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAg               <  1:4406154/93‑1 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgc              >  2:1450834/1‑93 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            <  1:4065169/93‑1 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            <  1:2251590/93‑1 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            >  1:3753072/1‑93 (MQ=255)

                                                                              aTTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGcgcg           >  1:297450/1‑93 (MQ=255)

                                                                              aTTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTGCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGcgcg           >  2:3571125/1‑93 (MQ=255)

                                                                               ttACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcgg          >  2:3436466/1‑93 (MQ=255)

                                                                                tACGGTATCCGCGGAGTTCGCTTTCTGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcggc         <  2:6479285/93‑1 (MQ=255)

                                                                                tACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcggc         >  1:3316473/1‑93 (MQ=255)

                                                                                  cGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGt       <  2:1926690/93‑1 (MQ=255)

                                                                                  cGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGt       <  2:5965507/93‑1 (MQ=255)

                                                                                   ggTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTa      >  1:2398617/1‑93 (MQ=255)

                                                                                   ggTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTAAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTa      <  1:5003445/93‑1 (MQ=255)

                                                                                    gTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAg     <  2:3487367/93‑1 (MQ=255)

                                                                                    gTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAg     <  2:1516379/93‑1 (MQ=255)

                                                                                     tATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAgg    <  1:5331769/93‑1 (MQ=255)

                                                                                      aTCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggt   <  2:5700948/93‑1 (MQ=255)

                                                                                      aTCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggt   <  2:4588878/93‑1 (MQ=255)

                                                                                      aTCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggt   <  2:4359049/93‑1 (MQ=255)

                                                                                       tCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggtg  <  1:14371/93‑1 (MQ=255)

                                                                                       tCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggtg  <  2:5021374/93‑1 (MQ=255)

                                                                                       tCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggtg  <  2:3910403/93‑1 (MQ=255)

GTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAGGTG  >  NC_000913/3596205‑3596384

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